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Abstract 



This research is concerned with the inter-particle potentials for few-particle bound state sys- 
tems in a scalar model with a Higgs-likc mediating held and QCD. The variational method, in 
a reformulated Hamiltonian formalism of QFT, is used to derive relativistic few particle wave 
equations and, consequently, the corresponding inter-particle potentials for stationary states of 
systems in question. The cubic and quartic contributions to the inter-particle potentials, in the 
non-relativistic limit, are expressed as multi-dimensional integrals. In the scalar Higgs-like model, 
when the mediator is massive, it is found that the non-linear terms modify the attractive Yukawa 
potential. In the case of QCD, the derived quark- ant iquark potential is significantly affected by 
the non-linear terms and provides a signal of confinement. 
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1 Introduction 



Physics - where the action is. 

Unknown 

1.1 The Present State of Affairs in Particle Physics 

The science of Particle Physics is concerned with the description of nature at the most fun- 
damental level. The foundation, on which it rests, was laid down with the invention of Quantum 
Mechanics (QM) and Quantum Field Theory (QFT). These two frameworks provide, as it is 
believed, the correct description of the physical laws which govern the microscopic world. 

QM mechanics is based on principles which are non-intuitive in the macroscopic world. A 
particle or a system of particles is described by state vectors; in coordinate space they are known 
as wavefunctions. If a state vector is known then the measurable observables, such as position, 
momentum and energy, can be determined as statistical averages. The concept of a trajectory of 
a particle in Classical Mechanics is replaced by probability density, derived from state vectors, of 
a particle to be at a specific point in space and time. QM can be extended to the relativistic realm 
for highly energetic systems in a natural way. A primary application of QM is in the treatment 
of atoms and molecules which has given chemistry its modern appearance. 

QFT is an extension of QM to situations where the number of particles, not necessarily el- 
ementary, may not remain conserved, i.e. particle creation and annihilation is permitted. Its 
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novelty is to treat the mediating fields on par with particle fields. It is motivated by the belief 
that the vacuum is filled with virtual particles which can become physical (i.e. go on-the-mass- 
shell) provided that all conservation laws remain stringent. The experimental study of elementary 
particles requires large amounts of energy, therefore particle physicists are appropriately inter- 
ested in relativistic QFTs. The structure of QFT comprises many branches of mathematics and 
to have a grasp of it requires many years of dedicated learning. 

The Standard Model (SM) of Particle Physics is a multi-component QFT which describes the 
elementary particles of nature and their interactions. The interactions among elementary particles 
arc mediated by the electromagnetic, weak, and strong mediating fields which are included in the 
SM. The SM has been established by conducting experiments which inter-played with theory to 
exalt the SM to its current form. The objective of the present day experiments, for instance 
the ones at the Large Hadron Collider, is to go beyond the domain of validity of the SM and 
possibly discover contradictions to the SM predictions. Discoveries violating the SM would signal 
the existence of the highly anticipated physics beyond the SM (BSM). More importantly such 
discoveries could lead to possible resolutions of the shortcomings of the SM. 

The first successful QFT was Quantum Electrodynamics (QED) where electrically charged 
particles interact via the mediating photon field. QED paved the way for the perturbative method 
of calculating quantum amplitudes for scattering cross sections and decay rates. The famous 
Feynman diagrams, historically stemming from perturbative QED, have become an indispensable 
tool in particle physics. In addition, the method of renormalization, the process of removing 
infinities in a systematic and meaningful manner, was also initially proposed in the context of 
QED. Basically, QED has become a template for all subsequent QFTs in particle physics. 

The early scattering experiments in the 1940s and 1950s discovered numerous, then thought to 
be elementary, particles known as hadrons. At that time, it was dubbed the "particle zoo" for the 
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lack of explanation of their large number. As more and more hadrons were being discovered, physi- 
cists began to suspect that they are not truly elementary (i.e. structureless) 1 . They organized 
these particles into various categories and introduced new quantum numbers for this purpose. 
Shortly, it was realized that symmetries played an essential role in categorizing hadrons. The 
Parton Model, in which constituent partons are treated as being entirely free inside of hadrons, 
and the advent of gauge symmetries led to the invention of Quantum Chromodynamics (QCD) 
where the fundamental massive quarks interact with each other via the massless gluons. 

Simultaneously, particle theory was readily extended beyond QED to account for the existence 
of neutrinos. More importantly, it unified the electromagnetic and weak interactions into a new 
Electroweak Theory (EWT). In this theory, neutrinos are electrically neutral elementary particles 
that interact via the weak interaction only. Presently, experiments indicate that there are three 
flavours of nearly massless neutrinos. The full incorporation of neutrinos into the SM remains 
problematic because of their non-zero mass. This predicament, along with many others, motivates 
to seek BSMs. In addition, the W and Z bosons, the quanta of the weak interaction fields, 
were discovered, as predicted by the EWT. There is much debate about the SM in the particle 
physics community but, despite some significant shortcomings and inconsistencies, it remains 
the most successful description of nature. For recent comprehensive summaries of the SM, see 
references [1,2]. 

The shortcomings of the SM may be resolved by the BSMs. A principal question is the lack 
of explanation as to why elementary particles have masses. It appears that to preserve the local 
gauge invariance of the theory all particles must be massless. The as yet undiscovered Higgs 
boson, and the spontaneous symmetry breaking mechanism associated with it, circumvents this 
requirement [3]. However, as the experimentalists sift through more data at the Large Hadron 

1 In fact, it is not certain whether the particles believed to be elementary today are really such. 
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Collider, the window on the Higgs boson mass is being gradually reduced. As an additional 
challenge, it appeals to theorists to explain the mass hierarchy of particles - the mass range from 
neutrinos to the top quark spans 11 orders of magnitude. Moreover, the experimental observation 
of the rotation of galaxies suggests that there are new heavy particles, the so-called dark matter, 
which do not interact via the electromagnetic or strong interactions and have not been produced 
in colliders to date. There is no dark matter candidate in the SM as it does not address the dark 
matter problem. Lastly, neutrinos are assumed to be massless in the SM, which is contrary to the 
experimental observation. 

A long-standing problem of theoretical physics is to reconcile General Relativity, the classical 
theory of gravitation, with QFT. The desire for a Quantum Theory of Gravity (QTG) originates 
from the ambitious endeavour to describe all physics in one unified theory. However, given the 
energy of the unification scale 2 Api anc k « 10 19 GeV and the collision energy of today's colliders 
Acoiiidcr ~ 10 3 GeV, it seems that in the near and intermediate future no experimental testing 
of QTGs is foreseen. Presently, the only possible arena of testing QTGs is in observational 
astronomy and cosmology. Now, it is plausible there exist new interactions and particles below 
Apiank thereby conceivably thwarting the present idea of unification in particle physics alone [4,5]. 
It can not be excluded altogether that the framework of current QFT itself becomes obsolete at 
some high energy scale, and consequently the idea of unification could become erroneous. In the 
author's view, it is rather premature to address the problem of unification given the present state 
of technology and knowledge. 

Returning back to the current state of affairs, the majority of particle physicists arc nonethe- 
less concerned with phenomenology and model-building, which are still very important fields. 



2 It is believed that the arrangement of the fundamental constants \ , which has the dimension of energy, 

V 

is actually related to the scale where the gravity effects on the microscopic world become large. 
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However, in the author's opinion, it is even more fundamental to be able to solve QFT fully. 
This might be the key to understanding QFT better. Perhaps, the next interaction among new 
particles to be discovered is non-perturbative as well as non-linear. One needs to know how to 
calculate cross-sections, decay rates and energy spectra reliably in such cases. These issues are 
already noticeable in QCD where the strong coupling and the non-linear mediating field have 
triggered new ideas and techniques of calculation. In the view of the author, the fundamental 
task is to develop universal methods of calculation, regardless of the coupling strength and the 
interaction types. If such tools are invented, it would be a major step forward. 

1.2 The Bound State Problem in QFTs 

The analytic solution of the Schrodinger equation for the hydrogen-like atom problem in QM 
served as one of the primary confirmations of quantum theory. The Dirac equation (and to a lesser 
extent the Klein-Gordon equation) being relativistic generalization of the Schrodinger equation, 
incorporates the ideas of anti-matter and spin (the intrinsic angular momentum of a particle) 
which have become an integral part of physics. This equation serves as a starting point in the 
development of realistic QFT. The Dirac equation has provided the correct 0(a 4 ) contributions 3 
to the energy spectrum of hydrogen [6] which can not be accounted for in QM alone. The Dirac 
and Klein-Gordon equations have become the standard pedagogy. 

The inter-particle interactions, entering into these bound state equations, arise from classical 
electromagnetism. The Coulombic potential describes such bound states adequately in the non- 
relativistic limit. The number of particles in these equations are fixed; yet there can be an arbitrary 
number of them. For instance, the Breit equation [7], being a many-particle Dirac-like equation, 
describes systems with more than one electron in the vicinity of a heavy nucleus. It includes a 

3 The perturbative expansion parameter a m Y37 ls cou P un S constant of QED. 
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Coulombic potential along with spin-orbit, spin-spin and the other 0(a 4 ) corrections. The main 
weakness of the Breit equation is that it is not invariant with respect to Lorentz transformations 
(i.e is not Lorentz covariant), and therefore it describes the relativistic effects perturbatively to 
order 0(a 4 ) (excluding held theoretical effects such as virtual annihilation). 

In QED, the mediators of the electromagnetic interaction are photons. Photons are massless, 
chargeless, spin one vector bosons which are classically described by the Maxwell equation. The 
Proca equation [8] is an extension of the Maxwell equation to include the case of massive vec- 
tor bosons. However, the principle of local gauge invariance (i.e. that is the QED Lagrangian 
must have no photon mass term in order to be gauge invariant) provides a sufficient theoretical 
explanation as to why photons should be massless. Thus, even though one can write the Proca 
equation for massive photons, it currently has little practical use. In addition, no bound states of 
photons are observed in nature. This is consistent with the absence of self-interaction terms for 
photons in the QED Lagrangian. In contrast, electrically charged massive particles can exist as 
free or in bound states. 

The search for a quantum field theoretic description of bound states of massive particles in 
QED has led to the Bethe-Salpeter (BS) equation [9,10]. It is a fully relativistic covariant equation 
which can be cast in integral or differential form. Being such, it is not surprising that the equation 
can not be solved analytically for any realistic problem. The usual perturbative expansion of the 
kernel of the BS equation (i.e. the interaction) and its truncation to the given order makes 
it only an approximate description. Furthermore, it contains negative energy solutions as well 
as relative time coordinates pertaining to the constituent particles. The BS equation produces 
accurate perturbative results for bound systems of two massive particles in QED. However, its 
applicability to three and more particles becomes increasingly difficult. This is all the more so for 
QCD [11]. 
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In the early days of particle physics, before the discovery of the strong interactions, there 
was no particular reason to study bound states in strongly interacting theories with the notable 
exception of heavy atoms. The establishment of QCD has profoundly changed this attitude 
because bound states are an important aspect in this theory. The interactions among the quarks 
(i.e. massive spin one-half quanta of the strong matter fields) and gluons (i.e. massless spin 
one quanta of the strong mediating field) originate from the non-Abelian gauge symmetry of the 
theory. It means that, in addition to the QED-like linear interaction between quarks and gluons, 
there are two non-linear gluon self-interaction terms of order three and four respectively. In light 
of these new terms, the effect on the running of the coupling constant in QCD is opposite to 
that in QED [12,13]. In QED, the running of the coupling constant is attributed to screening 
by virtual pairs hence the opposite behaviour in QCD is suggestive of anti-screening. At large 
energies the coupling strength decreases - a property known as asymptotic freedom. On the other 
hand, the coupling strength increases at lower energies and becomes non-perturbative at the scale 
Aqcd ~ 200 MeV - this is known as infrared slavery. 

The discovery of asymptotic freedom was deemed worthy of the Nobel Prize in 2004 (Gross, 
Politzcr, Wilczek) thus it deserves to be reviewed in some detail. The beyond-leading-perturbative- 
order corrections to the correlation functions 4 are ultravioletly divergent in most QFTs. Nonethe- 
less, for certain types of QFT, the ultraviolet infinities can be absorbed into the parameters of 
the theory (i.e. coupling constants and mass parameters) in a consistent and meaningful manner 
by introducing field renormalization. The behaviour of the correlation functions and the param- 
eters of the theory as the energy scale varies is described by the Callan-Symanzik differential 

4 Correlation functions are the expectation values of time-ordered products of field operators. 
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equation [14,15]. The Callan-Symanzik equation for massless QCD is 

G ( "< m) =0, (1.1) 



+ P{9)TT+niA{g)+m^{g) 



dlog/u, dg 

where n and m are the number of the quark and gluons fields in the correlation function, the 
parameter fi is the renormalization scale, the function /3(g) is simply entitled the "beta" function, 
whereas 7/1(5) and 7^,(3) are called the anomalous dimensions of the quark and gluon fields 
respectively. These functions can be calculated to a given order in perturbation theory. 

The massless QCD Callan-Symanzik equation has a solution provided the following subsidiary 
condition holds: 

5.9(A) 



d log A 



/?(<?), (1.2) 



where A = — — is the ratio of the renormalizcd A_r and bare Ao scales. This equation is known as 
the group renormalization equation and asserts that the beta function is just the rate of change of 
the coupling constant with respect to the change of scale. The beta function can be obtained by 
evaluating Feynman diagrams using a regularization scheme (such as dimensional regularization) 
to isolate the ultraviolet divergences. There are seven diagrams which must be calculated in order 
to extract the leading term of the beta function in QCD. The result of a detailed calculation 
yields [12,13]: 



7g 3 
(4tt) 2 



0(9) = -7J^ + O(g% (1.3) 



where the negative sign reflects the behaviour of the coupling constant as stated above. One can 
solve the differential equation (1.2) with the approximation (1.3) to obtain 

s 2 (A) = 2t4tT775 5T ' (i-^) 
l + 7^1og(A)/(8 7r) 

where go is the value of g measured at A = 1. According to this solution, the coupling constant 
is stronger for small A. Furthermore, it conforms with the experimental observation that the 



interaction is stronger at the lower energy scale. After all, the masses of the known hadrons 
formed by composite quarks and gluons seem to be all below about 10 GeV [2] . The running of 
the coupling is a theoretical hint that quarks and gluons should be confined. 

The phenomenon of confinement of quarks and gluons in bound states is an experimental 
fact. It should be emphasized that no free quarks and gluons have been found in nature. It is 
no exaggeration, then, to claim that the QCD Lagrangian is written in terms of the "wrong" 
degrees of freedom (i.e. quark and gluon fields instead of meson and baryon fields). Indeed, it is 
a big mystery why it should be like that. Moreover, the masses of the bound states of quarks are 
greater than the sum of its constituent quark masses. This is a characteristic of confinement and 
a defining property of QCD (but not QED). The theory indicates that it is important to develop 
non-perturbative techniques of calculation. A very ingenious method called Lattice QCD (LQCD) 
has been developed to study the hadron spectrum of the strong interaction precisely [16, 17]. 

In LQCD, continuous space-time is replaced by a finite grid of space-time points. Consequently, 
the action becomes discretized with the quark fields defined at each lattice point connected by 
the link variables representing the gauge bosons. The discretization requires a momentum cutoff 
on the order of the inverse lattice spacing a -1 to regularize the theory. The continuum limit can 
be approached as the spacing a is reduced to zero which retrieves the exact QCD. The LQCD 
formulation is in principle gauge invariant but the practical implementation requires an extensive 
computational effort. To calculate the mass spectrum of hadrons one numerically calculates the 
correlation functions for given quantum numbers and then fits them to exponential curves [18]. 
The results of LQCD are highly successful and are presently considered the best tool for describing 
the QCD bound spectrum. 

The confined bound states in QCD are a feature that still awaits its ultimate explanation. 
There is no rigorous and unambiguous explanation why quarks are confined. Identically, there 
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is no satisfactory explanation of the phenomenon of "string breaking" (i.e. the observation that 
hadrons do not separate into free quarks but rather create more hadrons if provided with sufficient 
energy). Non-relativistically, the existence of confinement must correspond to an appropriate 
shape of the quark-antiquark potential. The potential must rise at large separations to accord 
with the strong attraction and consequently explain the confining bound states. 

LQCD provides a way to calculate the potential for a static quark-antiquark pair by calculating 
the expectation value of the Wilson loop 5 [19]: 



where Z is the QCD generating functional, U is the link variable connecting neighbouring sites on 
a lattice, q and q are the quark and anti-quark fields respectively, S is the discretized Euclidean 
action, P denotes a particular enclosed path on the lattice and the trace is calculated over the 
gauge indices of Us. The static quark-antiquark potential can be then calculated via 



The calculation of the static quark-antiquark potential using equation (1.6) is a laborious numeri- 
cal task. Nonetheless, it has been done and the resulting curves can be found in references [17-19]. 

The form of the quark-antiquark potential has also been postulated in phenomenological mod- 
els. The original model [20], which has become known as the Cornell potential, serves as a good 
description for the low energy bound states of heavy quarks. A more general treatment of the 
mesonic bound states is presented in the work by Godfrey and Isgur [21]. They confirmed the 
existence of known and predicted hundreds, then-undiscovered, bound states of mesons by incor- 
porating relativistic effects in the lowest order. Another extensive discussion on the nature of 

5 In the early days, due to the complications induced by dealing with the quark fields, the quark fields were 
absent in the Wilson loop. 




(1.5) 



V(r) = - lim 1 Iog<W c >. 

T->oo I 



(1.6) 
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the quark-antiquark potential along with its spin, angular momentum and delta function correc- 
tions is provided by Lucha, Schoberl and Gromes in reference [11]. A more recent and realistic 
model [22] accounts for the string breaking effects by requiring the potential to soften at large 
distances (i.e. asymptotically reach a constant as the separation increases). Although such phe- 
nomenological models are effective, they all suffer from the fact that the underlying interactions 
arc inserted by hand. As already mentioned, it is an outstanding problem in particle physics to 
derive the quark-antiquark potential in QCD from first principles. It is a challenge that keeps 
theorists preoccupied. The numerous attempts and techniques of various complexity to derive the 
potential are discussed in the book by Greensite [23]. However, even Greensite's book is not a 
complete review of the subject; as the author says "to include every proposal would require an 
encyclopedia" ( [23] p. 207). There is a biannual conference devoted to the subject of "Quark 
Confinement and the Hadron Spectrum" . Proceedings of the conferences are published regularly, 
the most recent one being [24]; these present current updates of the status of the field. 

The research presented in this dissertation is concerned with the development of a technique 
for calculating the inter-particle potentials in QFTs containing non-linear mediating fields. The 
inter-particle potential provides a means of calculating the bound state spectrum and may indicate 
whether a QFT is confining or not. The existing method of determining the inter-particle potential 
based on LQCD, in the view of the author, is rather abstract and requires prolonged numerical 
calculations. It is then prudent to seek alternative methods to acquire the inter-particle potentials 
in QFTs. The basis of a such method, used in this dissertation, are discussed in the next section. 

1.3 The Method and Dissertation Outline 

The derivation of the inter-particle potential is implemented in the Hamiltonian formalism. 
The classical equations of motion and their formal solutions for the mediating fields are used 
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to reformulate the original Hamiltonian. The resultant Hamiltonian has a dependence on the 
mediating field via the Green functions only. There are two advantages to performing this refor- 
mulation. First, there is no requirement to quantize the mediating field since only bound states of 
particles are of concern. Second, the derived relativistic few-particle equations contain functions 
describing particle quanta only; without the reformulation the equations would be coupled in the 
fields that describe both particles and mediators. Thus, the complexity of the principal equations 
is reduced by the reformulated Hamiltonian. 

The primary subject of interest are theories in which the mediating fields are non-linear; in 
particular a scalar model with a Higgs-like mediating field as well as QCD. To study the effect of the 
non-linear interaction terms on the inter-particle potential the following scheme is implemented. 
The matrix elements of the Hamiltonian operator are calculated by using suitable trial states 
which describe few-particle bound systems or their elastic scattering (the latter is not addressed 
in this dissertation). It is possible to probe some or all interaction terms in the Hamiltonian 
by selecting different trial states. Once the matrix elements for the desired bound systems are 
found, the variational method is applied and relativistic equations, in momentum space, for the 
few-particle systems under consideration are obtained. These equations provide an approximate 
description of the bound state systems in question. In the Hamiltonian formalism, the bound state 
energy is bounded from below [25] and this provides the possibility for systematic improvement 
of the variational approximation. The interaction terms of these equations are examined as is 
described below. 

A scalar model containing a Higgs-like non-linear mediating field is studied in the second 
chapter. The results of reformulation and quantization are presented in some detail. The bound 
systems of a particle-antiparticle pair, three-particle, and four-particle systems are investigated 
at first. The trial states for these systems are mono-component in Fock-space and consequently 
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relativistic equations contain only one unknown function (i.e. the approximate wavefunction of 
the system being studied). The interaction terms of these equations are reduced to the non- 
rclativistic limit and Fourier-transformed to the coordinate representation. Hence, these terms 
are just the usual non-relativistic inter-particle potentials and those corresponding to the non- 
linear terms are in the form of multi-dimensional integrals. It is shown that the inter-particle 
potentials are dependent on the inter-particle distances only. Unfortunately, there are no simple 
closed-form expressions for arbitrary positions of the particles. However, considering particular 
cases with restrictions on the particle positions, it is possible to obtain a general picture of the 
full inter-particle potential. 

A trial state containing a one-pair (particle-antiparticle) and a two-pair components is in- 
vestigated thereupon. This trial state is capable, at least in principle, of describing the process 
of string breaking which has been mentioned in the context of QCD. The inclusion of the sec- 
ond component in the trial state leads to coupled relativistic equations for the functions sep- 
arately describing the one- and two-pair components of the bound state. The non-relativistic 
limit and Fourier-transform, along with an ansatz for the four-component wavefunction, yield an 
intricate expression for the inter-particle potential between the particle and antiparticle of the 
two-component of the bound state where the dependence on the four-component coordinates gets 
integrated out. This expression is in the form of a multi-dimensional integral. It is solved numer- 
ically to obtain an inter-particle potential which, for the case where the mediating Higgs field is 
massive, shows no evidence of confinement. 

The third chapter is devoted to QCD. It commences with a review of the Dirac equation 
and the quantization of QED followed by a succinct review of the QCD Lagrangian and its non- 
Abelian gauge symmetry. A similar scheme for quantization is followed for QCD as in QED, 
however, it is complicated by the addition of the various extra indices. The reformulation of QCD 
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is performed and the reformulated Hamiltonian is presented in a fixed gauge. Armed with the 
knowledge gained from the scalar Higgs-like model, a mono-component trial state describing a 
three quark system (i.e. a baryon) is considered first. This trial state has the right number of 
ladder operators to probe the cubic term of the QCD Hamiltonian. Surprisingly, it turns out that 
the summation over the colour indices (i.e. the colour factors) makes the cubic contribution of the 
matrix element vanish in the lowest order of iterative approximation. Next, a multi-component 
trial state describing a quark-antiquark system (i.e. a meson), and possibly string breaking, is 
considered. The summation over the colour indices indicates that only the non-linear cross terms 
of the interactions do not vanish. Again, the resulting potential energy terms in the equation 
that describes a non-relativistic meson are in the form of multi-dimensional integrals. Here, they 
are further complicated by the presence of momentum scalar vector products in the numerator. 
Unfortunately, it turns out that the only available technique to solve such multi-dimensional 
integral expressions, the Monte Carlo method, is incapable of producing stable numerical results. 
Thus, the best that can be done is to approximate these multi-dimensional integrals by upper 
bounds, for which the results are stable. 

It is worth mentioning that the variational method has been used effectively in the study of 
few-boson states in the scalar Yukawa theory [26,27], in the Higgs theory [28], in few-fermion 
bound states in relativistic quantum mechanics [29] and QED [30-32], and even in QCD [33,34]. 
Many results that appear in chapter two have been published in the preprint [35] . 

The long and cumbersome derivations and results have been relegated to the appendices to 
keep the contents more presentable. The summation convention on the repeated indices, where 
no summation sign is explicitly included, is implied throughout. 
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2 A Scalar Model with a Higgs-like Mediating Field 



We can't solve problems by using the 
same kind of thinking we used when we 
created them. 

Albert Einstein 

2.1 Reformulation 

Inter-particle potentials in few-particle systems described by a theory with the following non- 
linear scalar Lagrangian density 6 shall be addressed in the current chapter (units: h = c = 1): 

C = d v <j>* d v <f> — m 2 (j)* <f) — g \ 4>*4> — A (^» 2 + \ [d v X du X - f X 2 ) - \v X 3 ~ \° X 4 - (2.1) 

The quantities g, r\ are coupling constants with dimensions of mass, while cr, A > are dimension- 
less coupling constants. The parameters m and /i represent the bare masses of the scalar and 
mediating field quantum respectively. 

This Lagrangian density is Lorentz covariant and the mediating field \ mimics the gluon field 
of QCD if /j. = 0. In QCD, the cubic term contains a contraction of the gluon field and its 
derivative thus the entire term is covariant. There is no such contraction in the cubic term of the 
scalar Lagrangian density (2.1). Thus the coupling constant r\ is required to have dimensions of 
mass to maintain the correct dimensionality of the Lagrangian density. Undeniably, the \ ne ld 

6 Often in the literature a Lagrangian density is referred to as simply a Lagrangian. The correct terminology 
will be maintained throughout this dissertation. 
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is of the form of the Higgs field of the Standard Model. The A term in the Lagrangian density 
is required to keep the classical ground state energy bounded from below. It leads to a repulsive 
contact (delta function) inter-particle interaction for few-particle systems [36,37]. Thus, one can 
set A = since it has negligible effect on the results. 

The reformulation corresponds to a partial solution of the equations of motion. The equa- 
tions of motion corresponding to the Lagrangian density (2.1), derived from the action principle 
S J dx C = 0, are 

(d 2 + m 2 ) <j>(x) = -g<f>(x) X (x) (2.2) 
(d 2 +fi 2 ) X (x)= -g<f>*(x)<j>(x) - vx 2 (x) - a X \x). (2.3) 

Equation (2.3) has the integral representation (i.e. "formal solution") 

X (x) = X o(x) + J dx' D(x - x 1 ) p(x'), (2.4) 

where x — (t,x), dx — dtdx, p(x) = —gcj)*(x)(f)(x) — r]x 2 (x) — ax 3 (x) is the "source term" of the 
inhomogeneous equation, xo(x) satisfies the homogeneous (or free field) equation with the right 
hand side of (2.3) equal to zero, while D(x — x') is a covariant Green function (or the Feynman 
propagator) for the mediating field x, such that 

(d 2 + n 2 ) D(x - x') = 5 {i \x - x'). (2.5) 

Recall that the Green function can be expressed as 

/dk 1 
^-j^ tt, e - ifc < x -*>. (2.6) 
(2tt) 4 h 2 - k 2 + ie y J 

where e specifies the prescription for handling the singularities when integrating over k°. 

No processes involving free quanta of the mediating field will be studied, henceforth xo will 
be left out. Then, substituting (2.4) into (2.2), omitting terms containing xo, one obtains the 
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equation 

(d 2 + m 2 ) 4>{x) = g 2 cf>(x) J dx'<t>*(x')<f>(x')D(x,x') 

+ gr]<j>(x) J dx'x 2 {x')D{x,x')+ga<l>{x) J dx' x \x') D(x, x 1 ). (2.7) 

Equation (2.7) is derivable from the following Lagrangian density 

C = d tl (j)*{x)d^(j){x) - m 2 <j)*(x)(t>(x) + \g 2 I dx'4>*(x)<f>{x)D(x,x')(t)*(x')<p(x') 



+ 



grj (p* (x)<j>(x) J dx>D(x,x')x i {x')+g(T<p(x)<f>(x) J dx'D(x, a/)xV) (2.8) 



provided that the Green function is symmetric, i.e. D(x — x') — D(x' — x). 

For the case of a linear mediating field, i.e. when 77 = a = 0, the reformulated Lagrangian 
density (2.8) gives field equations that involve the particle fields only; the mediating field \ appears 
only through the propagator D(x, x'). Such a reformulated Lagrangian density, with 77 = a = 0, 
is convenient for the study of few-boson relativistic bound states in the scalar Yukawa theory as 
will be pointed out below. 

However, for the case of a non-linear mediating field, i.e. 77, a ^ 0, it is not possible to obtain a 
closed-form solution of (2.4) for the field \ m terms of D(x, x') and <j>, thus one requires to resort to 
approximations. An obvious one is an iterative sequence based on (2.4). The first-order iterative 
approximation corresponds to substituting the formal solution (2.4), with p = —gcj)*<j) (with xo 
left out), into the Lagrangian density (2.8). This yields the first-order approximate expression for 
the Lagrangian density 

1 



C = d„(f)* (x) d^4>(x) - m z cj)* (x)<t>(x) + 2 9 j ^ W^) 7 ^' %W {x)<t>(x') 
+ g7](t)*(x)<j)(x) J dx' dz'dz" D(x,x')D(x',z')D(x',z") p(z') p(z") 
+ ga<j)*(x)(f>(x) J dx'dz'dz"dz"' D{x,x')D(x',z')D(x',z")D(x',z"')p{z')p(z")p(z"'). (2.9) 
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The interaction terms in the Lagrangian density (2.9) do not contain the mediating field x explic- 
itly but rather implicitly through the mediating field propagators. The advantage of reformulating 
the Lagrangian density is that it enables one to use simple few-particle trial states free of the \ 
field quanta while still probing the effects of all interaction terms, including the non-linear terms 
(i.e. those with 77, a ^ 0). 

The Hamiltonian and momentum densities corresponding to (2.9) are obtained from the 
energy-momentum tensor in the usual way: 

dC 

where the index i = 1,2 stands for the fields </>i = </> and $2 = 4>* ■ The T 00 component of the 
energy-momentum tensor is the reformulated Hamiltonian density 

H = 0IT + 4>* n r -£ = H^ + H Il +Hi 2 +H l3 (2.11) 

where 

= n . n + (V0*) • (V$ + m 2 0*0, (2.12) 

U h = -^g 2 (f)*(x)<t>(x) J dx'<t>*(x')(t>{x')D(x-x'), (2.13) 

H l2 =-gr](i>*(x)<t>(x) J dx' dz' dz" D(x,x') D(x' , z') D(x' , z") p{z') p(z"), (2.14) 

H h =-g<T<p(x)<f>(x) J dx'dz'dz"dz"' D{x,x')D{x',z')D(x',z")D(x',z"')p(z')p{z")p{z"'). 

(2.15) 

The conjugate momenta of the fields are defined in the usual way 

Observe that every term of the reformulated Hamiltonian density has dimensions of M 4 (as it 
should). Notice that the Hamiltonian density is non-local as it contains integrals over space-time 
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variables. Equations (2.11)-(2.15) specify the Hamiltonian density that will be used to study 
few-particle bound states in this chapter. 

The T 0t components of the energy-momentum tensor are the momentum density components: 

v i = -n a t </>-n r a^*. (2.17) 

Note that the interaction terms do not contribute to the momentum density above. This expression 
will be used when the total momentum of the few-particle systems will be discussed. 

2.2 Formalism and Quantization 

In the Hamiltonian formalism of QFT, the basic equation to be solved is the 4-momentum 
eigenvalue equation 

P"|tf) = Q"|tf), (2.18) 

where P M is the energy-momentum operator of the quantized theory which follows from equations 
(2.11) and (2.17), = (E, Q) is the energy-momentum eigenvalue and |*) are the corresponding 
eigenfunctions. It is not possible to obtain exact solutions for the fi = component of equation 
(2.18). The exact solution would comprise all the possible basis states and would be expressed 
as an infinite sum over the components; such an approach is evidently intractable. Instead, the 
variational method will be used to obtain approximate results. For /i = 1,2,3 the solutions of 
(2.18) are free-field-like as is pointed out below. 

Variational approximations to the n = component of the eigenvalue equation (2.18), i.e. the 
Hamiltonian component, require the evaluation of 

8(* t \H-E\* t ) = Q, (2.19) 

where \$f t ) ar e suitably chosen trial states that contain adjustable features (functions, parameters). 
The accuracy of a variational approximation depends on the choice of the trial states. In general, 
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the more flexible features a trial state possesses the better are the approximate results that it 
yields. 

The canonical quantization of the theory with the Hamiltonian density (2.11) is performed 
next. First, one replaces the matter fields with their Fourier mode components: 

<t>{x) ^ j dp [(2tt) 3 2lj p ] ~ h (a (p) e"^* + fet ( p ) , (2.20) 
4>*(x)= J dp[(2^) 3 2^ p p (at(p)e^ + 6(p)c-^), (2.21) 

where Wp = m 2 + p 2 . Note that since the mediating field \ does not appear in the Hamiltonian 
(2.11) explicitly, one doesn't have to write out its Fourier mode representation. The next step 
is to promote the fields to the status of operators and impose the following non-vanishing, equal 
time commutation relations: 

[<f>{x),TL+{v)] = [4>*{x), IV {V)\ = i *(x - y). (2.22) 

All other commutators of the field and conjugate momentum operators vanish. The operators 
a(p) and b(p), in equations (2.20) and (2.21), are the free particle and antiparticle annihilation 
operators respectively whereas the operators a^(p) and tf(p) are the free particle and antiparticle 
creation operators. The vacuum state |0) is defined by a(p)|0) = 6(p) |0) = 0. Note that the 
canonical quantization is performed in the interaction picture. 

Expressing the creation and annihilation operators in terms of the fields and its derivatives, 
the commutation rules become 

[a(p),at(q)] - [&(p),& f (q)] = S(p - q). (2.23) 

The vacuum energy question is not essential to the bound state problem addressed in this dis- 
sertation, thus the creation and annihilation operators in the Hamiltonian are normal-ordered. To 
obtain the Hamiltonian operator the spatial coordinates are integrated out from the Hamiltonian 

20 



density (2.11): 

H(t) = J dx : W(t,x) :, (2.24) 

Next, the Hamiltonian is expressed in terms of the particle and antiparticle operators a, 6, a) and 
&t . This is straightforward for the free- field part of the Hamiltonian as can be seen in standard 
books on QFT, see for instance [38,39]. However, obtaining the analogous expressions for the 
interaction terms is tedious and this has been accomplished with the use of a computer code. 

The momentum operator is similarly obtained by normal-ordering and integrating out the 
spatial coordinates from the momentum density (2.17): 

P = I dx : V(x) : . (2.25) 



Notice that since the momentum density operator docs not involve the interaction terms, it is 
time independent. 

For the description of stationary (i.e time-independent) few-particle bound states, it is con- 
venient to switch to the Schrddinger picture in which the time-dependence of the Hamiltonian is 
removed. The two pictures are related by 

|tf/(t)) = e iH ° t \* s (t)) (2.26) 
where H is the free-field part of the Hamiltonian which in this theory is given by equation (2.12). 

2.3 Particle-Antiparticle State 

The simplest particle-antiparticle trial state can be expressed in terms of Fock-states as follows: 
|* 2 ) = f dpi,2 F(P1, 2 ) a f (Pi) & f (P2) |0>, (2.27) 
where i^pi^), is an adjustable coefficient function which is determined variationally 7 . 



7 The subscript notation means that dpi,2 = <2pi <ip2 and F(pi ; 2) = F(pi, P2). This notation is used throughout 
the dissertation. 
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(a) 

Figure 2.1: The one-chion exchange (a) and virtual annihilation (b) diagrams of the Yukawa 
kernel equation (2.29). 

To implement the variational principle (2.19) the matrix element (^l H — E \^2) is worked 
out and is varied with respect to the adjustable coefficient function F*. This leads to the following 
particlc-antiparticle relativistic wave equation in momentum space: 



F (pi, 2 ) (w Pl +uj P2 -E) = f dp' 12 J 7 2,2(pi, 2J Pi,2)i ;i (pi i 2) J 



(2.28) 



where 3^2,2 is the relativistic Yukawa particle-antiparticle interaction kernel (inter-particle "momentum- 
space" potential). In the interaction picture, it is given by 

i(« p ; +W p / -W P1 -Upj)t 



^2,2 (p'i, 2 ,Pl,2) = , To , ^— (5(Pi +P 2 - Pi -P 2 ) 

8(2tt) 3 ^/wp-wp^wp^p, 



+ 



1 



+ 



1 



+ ■ 



(2.29) 



V - (Pi - PI? ' - (P2 - P2? ' V 2 - (PI + P2? ' M 2 - (Pi + P 2 ) 2 - 

The first two terms in the square brackets in equation (2.29) correspond to one-mediating-field- 
quantum exchange (or one-"chion" exchange) and the last two to virtual annihilation Feynman 
diagrams. These diagrams are shown in Figure 2.1. 

The time dependence of the particle-antiparticle interaction kernel 3^2,2 can be "rotated away" 
by the use of equation (2.26). If the trial state \$> 2 ) = c lHat \^2) is used in the above, then 
the corresponding interaction kernel in the Schrodinger picture will be time independent, i.e. 
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one obtains equation (2.29) with t = 0. Therefore, the Schrodinger picture (time independent) 
interaction kernels can be found, in effect, by setting t — in all interaction picture kernels. 

Unfortunately, the quantized version of the non-linear terms (2.14) and (2.15) of the Hamilto- 
nian density are not probed by the particle-antiparticle trial state (2.27), i.e. the matrix elements 
(^2 1 Hi i |*2) vanish for i = 2, 3. The vanishing occurs because of a mismatch in the number of 
the ladder operators between the trial state and the interaction terms. More elaborate trial states 
than (2.27) must be used in order to examine the effects of non-linear terms. As it is, with the 
trial state (2.27), the problem reduces to the scalar Yukawa model [26]. 

For the particle-antiparticle trial state (2.27) to be an eigenstate of the momentum operator 
(2.25) we require that P \^f 2 ) = Q | ^2) - In the rest frame, where Q = 0, this requires that 
^(Pi,2) = <5(Pi + P2) /(Pi)- Consequently, an energy calculation yields the rest energy (i.e. the 
mass) of the particle-antiparticle system. The centre of mass motion separates and the relativistic 
momentum-space particle-antiparticle wave equation simplifies to 

1 



m 2 + (p'-p) 2 -K<-^p) 2 

1 \ 1 ( 1 



. (2.30) 



2 \fi 2 -4w2y 2 \ yu 2 — 4 cj 2 , 
In the non-relativistic limit (i.e. p 2 << m 2 ) the Fourier transform of equation (2.30) yields 
the expected Schrodinger equation in coordinate space for the relative motion of the particle- 
antiparticle system: 

--V 2 ?/>(x) + l/(x)V(x) =eV(x), (2.31) 
m 

where e = E — 2mis the non-relativistic energy, and the potential energy depends on the particle- 
antiparticle distance: 

^) = -^^-7^W), (2.32) 

g 2 

where a„ = — is a dimcnsionless coupling constant. The first term, due to the one-chion 
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exchange, is the usual Yukawa potential and is always attractive (i.e. gravity-like) in this scalar 
theory. The second term, due to the virtual annihilation, can be either attractive or repulsive 
depending on the values of the parameters \i and m. It is a correction to the Yukawa potential 
that is a feature of the quantum field theory and has no analog in relativistic QM. 

The relativistic equation (2.30) is not analytically solvable, so approximation methods must 
be used. Perturbative and variational approximations are presented in the paper by Ding and 
Darewych [26], along with a comparison to results obtained using the ladder Bethe-Salpeter 
equation and various quasi-potential approximations. 

The particle-antiparticle trial state has served as a preamble. In the following sections, trial 
states that probe the non-linear terms are examined. 

2.4 Three-Particle State 

To observe the effects of the non- linear terms of the Hamiltonian density (2.11) on the inter- 
particle potential, one must consider trial states with either more particle content or more Fock- 
space components. The easier task is to consider a three identical particle trial state given by 

l* 3 ) - J dpi,2,3 ^(Pi,2, 3 ) a f (pi) a f (p 2 ) « f (p 3 ) |0>, (2.33) 

where -F(pi,2,3) is a three-particle function to be determined variationally. Note that this trial 
state can be taken to be an eigenstate of the momentum operator (2.25), i.e. P IW3) = Q ^3) 
with the choice F(pi,2,3) = £(pi +P2 + P3 — Q) /(Pi,2), where Q is the constant total momentum 
of the three-particle system. Therefore, the wavefunction will be of the form where the centre 
of mass motion is completely separable, just as was the case for the particle-antiparticle system. 
If Q = 0, the the corresponding eigenenergy will represent the rest mass of the three-particle 
system. 
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The addition of an extra third particle operator to the trial state significantly complicates the 
evaluation of the matrix element. However, the implicit symmetry of equation (2.33) with respect 
to interchanges of momentum variables due to the identity of the particles makes the task easier. 
It permits one to carry out the evaluation of the matrix element for the three-particle trial state 
in terms of the symmetrized function F$: 

6 

F s (pi,2, 3 )= E F (Pn^s), (2-34) 

il .12 ,*3 

where the summation is on the six permutations of the indices 1,2 and 3. Making use of this 
symmetrization enables one to keep the arguments of the kernels non-permuted while storing all 
information about the symmetry under interchanges of the momentum variables in F$- 

The matrix element (^ 3 | H — E\^ 3 ) is calculated by commuting the ladder operators and 
turning them into momentum delta functions. Once the commutation is completed, the momen- 
tum integrals can be found in a straightforward, although involved, manner. Next, the variational 
derivative with respect to F* is carried out and set to zero (The complete expressions for the 
matrix element as well as other intermediate steps in the calculations are given in Appendix A 
section 5.1). This yields the following relativistic momentum space integral wave equation for 
stationary states of three identical particles: 

fs(Pl,2,3) (w pi + W P2 + W P3 -E) = J dp' 12 y3,3(pi,2,3) Pl,2,3) *s(pi,2,3) 

+ yVl,2,3C3,3(Pi,2,3,Pl,2,3)i^Pl,2,3)- (2-35) 

The relativistic interaction kernels (i.e. the relativistic momentum-space inter-particle potentials) 
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Figure 2.2: The three-chion propagator vertex corresponding to the cubic interaction kernel C3.3 
equation (2.37). The two propagators on the left should actually overlap (impossible to draw) 
such that they are perpendicular to the direction of time. 



are given by 

3^,3 (Pl,2,3iPl,2,3) 



8(2tt) 3 



U,i2,«3 

^3,3 (Pl,2,3)Pl>2,3) 



S (Pl +P2 Pn - P» 2 )^(P3 ~ Pis) 



M 2 - ip'i -Pti) 2 



(2.36) 



g 3 rj 



— T 



8(2tt) 6 

«1,«2,»3 

1 



<^(Pi + P2 + P3 - P»l - P»2 ~ PjJ 

V W P'i W Pi W P5 U Pn% W Pi3 



1 



1 



(2.37) 



_M 2 - (Pi + P2 -Pii -Pii? M 2 - ip'i - Pn) 2 V 2 - O2 -PiiY 
where the summation is on the six permutations of the three indices 1, 2 and 3. It is evident from 
the covariant factors of equation (2.36), that the kernel 3^3,3 corresponds to the three inter-particle 
one-chion exchange interactions. There are no virtual annihilation terms since only particle opera- 
tors (and no antiparticle operators) are present in the trial state (2.33). The kernel £3^, equation 
(2.37), corresponds to the non-linear interaction term Hi 2 (see section 5.1 of Appendix A for 
details). It is similarly evident, from the covariant factors of equation (2.37), that the Feynman 
diagram corresponding to the cubic kernel Cz,3 contains a three-chion propagator vertex. This 
diagram is shown in Figure 2.2. Note that the three-particle trial state (2.33) does not probe the 
quartic interaction term %/ 3 equation (2.15), i.e. Hi 3 ^3) = 0. 
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Equation (2.35), together with the kernels (2.36) and (2.37), is a relativistic wave equation 
for stationary states of a system of three identical scalar particles. It can describe purely bound 
states of the three particles (spinless bosons) or elastic scattering among them, but not processes 
involving the emission or absorption of chions. The relativistic kinematics (i.e. kinetic energy) 
of the system is described without approximation, but the relativistic dynamics (i.e. potential 
energy) is described only at the level of one-chion exchange between the particle pairs by the 
relativistic kernel 3^3,3, along with a three-particle iterative first-order correction C 3j3 due to the 
non-linear interaction terms. 

In principle, one would wish to solve the relativistic three-particle equation (2.35). However, 
this is a formidable task which cannot be done exactly. Even the determination of approximate so- 
lutions is a very challenging task which will not be undertaken in this dissertation. Instead, it is of 
interest to consider the non-relativistic limit of equation (2.35), which in the coordinate represen- 
tation is just a non-relativistic three-particle Schrodinger equation and the Fourier-transforms of 
3^3,3 and £3,3 are the inter-particle potentials. Among other things, these inter-particle potentials 
can be used to calculate the bound state energy for systems of three heavy spinless particles. 

The non-relativistic limit of equations (2.35)-(2.37), just as for the particle-antiparticle system, 
is obtained by assuming p 2 << m 2 and then Fourier-transformed to coordinate space. From 
equation (2.36), the non-relativistic inter-particle potential for the three-particle system due to 
the Yukawa interaction is, as expected, 

{g— Q — P.X13 0-/^23 I 
+ + }. (2.38) 
X12 Xl3 X 2 3 J 

Unlike in the particle-antiparticle case (2.32), no virtual annihilation delta function contributions 
appear in this expression since there are no antiparticles present in the three identical particle 
system. 
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The non-relativistic "cubic" potential term that follows from equation (2.37) is 



T/ , , 3 /\ n'e-' q -' X - £(qi+q 2 + q3) ... 

Vfc(xi A a, M) = * _/ dq^ ^ - qf) ^ - qg) ^ + qg) , (2.39) 

3 g 3 ?7 

where a„ = — 5 — 5- is a coupling constant with dimensions of mass (sec section 5.1 of Appendix 

2boiT°m i 

A for details). Vc can be simplified to a three-dimensional quadrature (see section 5.1 of Appendix 
A for details): 

/ e -/x|xi+x| g-/i|x 2 +x g— /^|x 3 +x| 
|xi+x| |x 2 + x| |x 3 +x| 

This is an overall well behaved convergent integral for any fi > 0. However, it cannot be evaluated 
analytically in general. The expression for Vc, equation (2.40), with the substitution v = x + xi, 
can be written as 

/Q~V\v\ g-(u|v+X 2 l| C -M|v+X 3 i| 
dv ~ TT - y— p -j— p , (2.41) 
|v| |v + X 2 l| |v + X 3 i| 

where x 2 i = x 2 — xi and x 3 i = x 3 — xi. The expression for Vc, as shown in Appendix A section 
5.1, can also be written Vc in the form 

e - M 2 (/3i+/32+/3 3 ) 



Vc(xjj,/i) = -a^-n I dp 



x exp ( -ML±^L±Ml ) , (2 .42) 



' ~4( U-h h »': + h h) 
which shows explicitly that Vc depends only on inter-particle distances Xij = Xji = |xj — Xj|. In 
light of equation (2.42), one can conclude that the total potential energy function V = Vy + Vc 
depends only on the inter-particle distances, and it is invariant under 3D rotations and translations 
of coordinates as is expected of a closed system. Therefore, to emphasize this fact the notation 
Vc{xij,fj) is used in what follows. In addition, the representation (2.42) is suitable for numerical 
evaluations of Vc for arbitrary values of Xij . 

Unfortunately, it is impossible to render multi-dimensional plots of inter-particle potentials 
of three independent variables X12, X23 and xi 3 , even if they can be worked out numerically for 
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arbitrary Xij. The best one can do is to plot one and two dimensional sections, corresponding to 
certain restrictions on the coordinates, to gain some idea of the shape of the cubic term Vc- 

Case 1: V(xi )2 ,2,A*) = V(ar 2 i,/i) 

The first case is for /x > and x 2 = x 3 (i.e. two particles overlap) so that there is only one 
inter-particle distance, X21 = :ri 2 = |x 2 — xi| = |x 3 — xi|, to be concerned with. The integral for 
Vc, equation (2.41), in this case is expressible in terms of the special function, the exponential 
integral, defined by 

roo -t 

E 1 (z) = J —dt, (2.43) 

namely 

Vc(xi, 2j2 ,/i) = 

Vc(x 21 , f i) = -^^\e-^^[\n(3)-E 1 (x 2W )}+e x ^E 1 (3x 2W ) I. (2.44) 

X2lH [ j 

The length is chosen to be expressed in units of the Bohr radius and, correspondingly, the 

ma g 

energy in units of ma g . Equation (2.44) can be written in terms of the dimensionless variables 

r = x 21 ma g and M = - . Suppressing the singularity at x 2 = x 3 in the Yukawa term (2.38), 
ma g 

the total inter-particle potential V = Vy + Vc for the three-particle trial state in the case when 
fi > and x 2 = x 3 is 



V r (x 1 , 2 , 2 , M >0)-^(r,M) 

{„-rM /p-rM -r M r M \ ~) 

- 2— - i (7^ M3) - W E l( rM) + _ El (3rJ»f)J j, (2.45) 

4 7r 4 a 12 n 

where m = tt- = is a dimensionless constant. Figures (2.3) and (2.4) are plots of 

ma 2 g g 

V(r, M)/ma q , equation (2.45), as a function of r for the different values of M, with m = and 
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Ki = 0.1 respectively. It is apparent that the contribution of Vc lowers the overall potential with 
decreasing values of M. 

For the case where fj, = 0, Vc, equation (2.39), becomes logarithmically divergent and requires 
regularization. There are various ways of regularizing Vc with /j, = 0. One way is by subtracting 
off an infinite constant and absorbing it in the bound state energy, namely: 

Vc(xij,fJ.) = V c {xij,/J.) - V c (aij,fi) (2.46) 
e -> e = e- Vc(ay,/i), (2.47) 

where a,j = a; — aj are reference (i.e. constant) inter-particle vectors and e is the eigenenergy 
corresponding to e in equation (2.31). In the limit ^ = 0, such a procedure renders Vc to be 
finite. 

The inter-particle potential for the case when x 2 = x 3 and \i = was originally derived using a 
different method by Darewych and Duviryak [40] . One can retrieve this expression from equation 
(2.45) by expanding in the limit M — > 0, namely: 

Vc(xi, 2> 2, n > 0) = V c (r, M) = 4 tt 4 a, ( In(3) - 1 + 7 + Mr M ) ) . ( 2 - 4 8) 

where 7 = 0.5772 is the Euler's constant. Consequently, implementing a regularization, as 
prescribed by equations (2.46) and (2.47), with the reference separation ai\ma g — 1 (where 
021 = |a 2 - ai|) yields: 

Regularized (Xl,2,2, A* = °) = Regularized (r, A* = 0) = ma 2 g j - ^ + K x lnr j. (2.49) 

2 

Notice that V(r,fj, — 0) is dominated by the Coulombic term for small r but by the logarithmic 

r 

term for larger r. 

However, the results for the \i — case are regulator-dependent. Within a perturbative 
calculation, it is known that there are other infrared singular effects which lead to a cancellation 
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Figure 2.3: The inter-particle potential V(r,M) / ma 2 g of equation (2.45) as a function of r 
X21 ma g for M = 10~ 10 with m = (i.e. Vy only) and K\ = 0.1 (i.e. Vy + Vc)- 




Figure 2.4: The inter-particle potential V(r, M) / ma 2 g of equation (2.45) as a function of r 
X21 ma g for k\ = 0.1 and various values of M. 
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of all regulator-dependent terms. It is expected that a more inclusive ansatz would do the same. 
Therefore, the results of equation (2.49) are most likely unphysical. Nevertheless, in the theory 
with r\ = (i.e. without the cubic term in the Lagrangian density (2.1)), infrared-divergences do 
not arise, as is pointed out in Section 2.5. 

Case 2: V(x 2 i = x 31 = .t 23 = A,/i) = V(A,/x) 

Another case where there is only one distance argument in the inter-particle potential is when 
the coordinates are at the vertices of an equilateral triangle. Unfortunately, no analytical solution 
is possible in this case and one has to carry out a numerical integration. The numerical integration 
was performed with the GNU Scientific Library [41]. With the previous choice for the units of 
length and energy, the expression for Vc is 



V(A, n) = V{r, M) = —mail 3 



c -rM 



r 

, M_ f°° ,n e -M^ 1+fe+fe ) / r 2 (/3 1+ fe + fe) \| 

+ 47ri dPl ^ 3 (A/32 + A/33 +/3 2 /? 3 ) 3 / 2 P { + hfo + W ) j 1 ™ ] 

where r = Ama g is the dimensionless inter-particle distance, A — |x 2 — xi| = |x 3 — xi| = |x 3 — x 2 | 

and M = — - — is the dimensionless mass parameter of the mediating field. A plot of equation 

ma g 

(2.50) is given in Figure (2.5). The overall character of the potential retains the same features as 
in the particular case x 23 = 0, equation (2.45). The curves show no hint of confinement as they 
all asymptotically approach zero for this M > case. 



Case 3: Arbitrary inter-particle distances 

Finally, the inter-particle potential for arbitrary inter-particle distances is discussed in this sec- 
tion. There are no analytical solutions available for Vc, equation (2.40), and the potential is calcu- 
lated numerically with the help of the GNU Scientific Library. Using the Gaussian parametrization 
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Figure 2.5: The inter-particle potential V(r,M) / mag of equation (2.50) for k\ = 0.1 and the 
indicated values of M = ^ , where r — Ama g and A = x 2 \ = £31 = x 2 3- 



given by equation (2.42) and with the previous choice of units, one obtains 



V(xij , fi) = V{r l:j ,M) = -ma 2 



e -r 12 M e -r 23 M e ~r 13 M 



9 



ri2 r-23 ri3 



«1 f°° ja e -M^( A+&+& ) ( r 2 2f3i+r 2 3f32+r 2 3/33 



+ 4tt J df)l ' 2 ' 3 (ftft + A/33 + /3 2 /3 3 ) 3 / 2 6XP { W2 + Pifo + ftft) ) j (2 ' 51) 
where = Xij ma g are the dimcnsionlcss inter-particle distances and Xij = |xj — xj\. Four 
surface plots of equation (2.51) for different values of ri3 are given in Figure 2.6. Note how the 
inter-particle potential rises with increasing separations among particles. However, these plots 
show that the inter-particle potential is of non-confining nature as the curves do not cross the 
zero plane. 
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Figure 2.6: The inter-particle potential V(rij,M) / ma 2 g of equation (2.51) for K\ — 0.1, M 

— — 0.1 and ri3 = 0.01(top left), 0.1(top right), 1.0(bottom left), 10(bottom right) where r i7 - 
ma g 

Ay ma 9 . 
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2.5 Four-Particle State 

It has been shown above that the cubic Hi 2 interaction term of the Hamiltonian affects the 
inter-particle interaction and consequently the energy spectrum for a three-particle system. How- 
ever, the three-particle trial state (2.33) does not probe the quartic T-Li 3 interaction term of the 
Hamiltonian. Thus, one must consider a system of four identical particles and examine how both 
non-linear terms of the Hamiltonian, Hi 2 and Hi 3 , affect the inter-particle potential. 

The four identical particle trail state analogous to equation (2.33) is given by 

|* 4 ) = J d Pl . A F(pi.. 4 ) a+(pi) a f (p2) a f (p 3 ) a f (p 4 ) |0>. (2.52) 

For this four-particle trial state to be an eigenvector of the momentum operator (2.25), one 
requires that P ^4) = Q ^4) which can be achieved with the choice F(pi, A ) = 5(p 1 + p 2 + 
P3 + P4 — Q) /(Pi,2,3)- Thereupon, as for the particle- antiparticle and the three-particle cases, 
the wavefunctions will be of the form where the centre of mass motion is completely separable for 
this identical four-particle system. It is convenient to write everything in terms of the completely 
symmetrized function F$ because the trial state (2.52) is completely symmetric under interchanges 
of the momentum variables. The symmetrized function is 

24 

F s (pi.a)= F (P«i^3,u) (2-53) 

»l,*2)»3i*4 

where the summation is on the 24 permutations of the indices 1,2,3 and 4. The matrix element 
(\I>4| H — E 1^4) is calculated and the variational derivative with respect to F* is found and set 
to zero (refer to section 5.2 of Appendix A for details). This leads to the following four identical 
particle relativistic equation for the function Fs (in momentum space): 

^MPl..4) (w pi +UJ p2 +W P3 +W P4 - E) = j dp' h2 3 ; 4,4(pi..4>Pl..4)i ? s(pi..4) 

+ J dp' 1 .. 4 C4,4(p 1 ..4,Pl..4)i ; s(pi..4) + J d Pl.A Q4,4(Pi.. 4 ,Pl..4)F5(pi..4)- (2-54) 

35 




Figure 2.7: The four-chion propagator vertex corresponding to the quartic interaction kernel Q4.4 
equation (2.55). The two propagators on both sides should actually overlap (impossible to draw) 
such that they are perpendicular to the direction of time. 

The Yukawa and cubic interaction kernels 34,4 and €4,4 are similar in structure to those of the 
three-particle trial state except there is dependence on an extra momentum coordinate (details in 
section 5.2 of Appendix A). The relativistic quartic interaction kernel is 



Q4,4(Pi.. 4 ,Pl..4) 



ffV 1 
16 (2tt) e 



24 

E 



<KPi + P2 + P3 + P4 Pi P2 P3 P4) 



UJp' 2 UJp' 3 UJp' i UJ pi U} p2 UJ p3 UJ pi 



1 



M 2 - (p'i + P2 + Vz - Pn - Pt 2 ~ Pt 3 ) 2 
1 1 



1 



(2.55) 



M 2 - (p[ ~Ph) V 2 - (p' 2 -Pt2) 2 V 2 - (Ps-Pis) 2 . 
where the summation is on the 24 permutations of the indices 1,2,3 and 4. This contribution 
to the potential corresponds to a four-chion propagator vertex shown in Figure 2.7. The non- 
relativistic limit of the interactions in the coordinate representation is examined below. 

In the non-relativistic limit, the Yukawa and cubic kernel reduce just as their three-particle 
trial state counter-parts. There are six Yukawa interaction terms which can be worked out in 
analytical form: 



VV(xi.. 4 ,m) = -Oig 



X12 



Xl3 



X14 



X23 



X24 



X34 



(2.56) 



where a a — 



9 16 7TTO 2 



is the dimensionless coupling constant and Xij = Xji = \xj — Xi\ are the 
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inter-particle distances as before. 

The cubic interaction kernel for the non-relativistic four-particle case reduces to a sum of four 
terms, i.e. one for every three-way interaction: 

Vc( x i..4,m) = ~ot v J dqi..4 

J- n 4 c - tqi .x I(5(qi+q2+q3 ^ (q4) nie- i q' x 'J(qi+q 2 +q 4 ) ( 5(q3) 
I (M 2 +qi)(^ 2 +qi)(M 2 + q^) (v 2 + q?) (m 2 + q§) (m 2 + q|) 



n 4 e - zq ,. Xl(5(qi+q3+q4)(5(q2) n 4 e _ iq ,. Xl|5(q2 + q3 + q4)(5(qi) 
(M 2 + q 2 )(M 2 + qi)(M 2 + q|) (^ 2 + q§) (m 2 + q§) (m 2 + q|) 



3 g 3 r? 

where a„ = — ^ — ? is a coupling constant with dimensions of mass (refer to section 5.2 of 

2567r t, m J 

Appendix A for details). The delta function with a single momentum variable in every term 
indicates that the particle carrying that momentum is a "spectator" of the three-way interaction. 
This expression is obtained by following similar steps as those leading to equation (2.39). 

The quartic interaction kernel for the four-particle system in the non-relativistic limit reduces 

to 

v ( \ I a n^" q '- x ^(qi + q2 + q3 + q4) 

3g 4 a 

where a a = — j is a dimensionless coupling constant. This expression is obtained by using 
10247r y m 4 

similar steps as those leading to equation (2.57). The cubic Vc and quartic Vq contributions to 
the inter-particle potential are both symmetric under 3D rotations and translations even though 
it is not readily evident from equations (2.57) and (2.58). In Appendix A, it is explicitly shown, 
using Gaussian parametrization, that these equations indeed depend only on the inter-particle 
distances, i.e. Vc = Vc(xtj) and Vq = Vg(xij). 
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The expressions for Vc and Vq can be simplified to three-dimensional quadratures: 



Vc(xi.. 4 ) = - a v 7T 3 J dx | 



g— (U|X1+X| C -^|X 2 +X| g— /^|X3+X| g-(U|xi+X c -^l|x 2 +X c -^l|x 4 +X 

|xi+x| |x 2 +x| |x 3 +x| + |xi+x| |x 2 + x| |x 4 +x| 



g-ju|xi+x| c -^|x 3 +x| c -^i|x 4 +x| g- J u|x 2 +x| g— /^|x 3 +x| c -^i|x 4 +x| j 

+ n n n r + n n n r ^ ( 2 - 59 ) 

|xi+x| |x 3 +x| |x 4 +x| |x 2 +x| |x 3 +x| |x 4 +x| j 

/„-/i|xi+x| e -p|x 2 +x| e -fl\x 3 +x\ c -^i|x 4 +x| 
J " |x 1+ x| |x, + x| |x, + x| l2m > 

Here, analogously to the three-particle case (2.40), the integral expression for Vc, equation (2.59), 
is convergent for /j, > 0. However, if yu = the integrand of Vc goes as |x| _1 for large |x| and 
the integral diverges. The integrand in Vq, equation (2.60), behaves as |x|~ 2 for large |x| and 
the integral remains finite for fi > 0. It is shown below that the contribution of Vq to the total 
inter-particle potential of the four identical particle system is a diminishing correction as the 
inter-particle separations increases. On the other hand, for small separations, depending on the 
values of the coupling constants, it can have a significant effect. 

The expressions for Vc and Vq, equations (2.59) and (2.60), analogously to equation (2.41), 
can be written as 

e ~i u l v l f t .-/ J l v + x 2l| g-Ji|v+x 3 i| e -M|v+x 2 i| c -Ai|v+x 4 i| 

Vc(xi..4,m) =-71-% ldv— r - r {- ■ ■ - + - ■ -- ■ 

" + X21 V + X31 V + X21 V + X41 



f c — I v I ( e~' 
(x!..4,m) =-tt a, Jdv-^-l — 



g-/^|v+x 3 i| g-(u|v+X4i| e -M|v+x 32 | 6 -m|v+x 42 | I 

+ n r + i ri rh ( 2 -61) 

|v + x 3 i| |v + x 4 i| |v + x 32 | |v + x 42 | J 

/ e — A* | v | g-(u|v+x 2 i| g-ju|v+x 3 i| e -/^|v+x 4 i 
dv ill, r> ( 2 - 62 ) 
|v| |v + x 2i | |v + x 3 i| |v + x 4 i| 

where Xjj = Xj — Xj . It is impossible to obtain analytical expressions for Vc and Vq in general 
and one must resort to numerical evaluation of these integrals. However, as for the three-particle 
system, there is a particular solvable case, namely when Xi = x 3 and x 2 = x 4 , which shows the 
general features of the total potential energy. With this restriction on the coordinates there is 
only one inter-particle distance x 2 i = |x 2 — xi| = |x 4 — x 3 | to be concerned with. 
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Special Case : V(xi i2 ,i,2, A 4 > 0) = V(x 2 i,^ > 0). 

For the case when xi = x 3 and x 2 = X4 with yu > the cubic potential term Vc of equation 
(2.61) leads to four copies of what has been previously obtained for the three-particle case, namely, 

Vc(xi, 2 ,i,2, fJ- > 0) = V c (x 2 i,v > 0) 

8 7r 4 a 



1) e -Z2lM 



[In (3) - Ei (.T 2lA t)] + c x ^ Ei (3ar 2 i/i) } ■ (2.63) 



The single difference is in the numerical factor of four in front. 

The quartic potential term Vq, equation (2.62), can be reduced to an evaluation of a single 
quadrature (refer to section 5.2 of Appendix A for details): 



Vq(xi i2j i i2 ,/x > 0) = V Q (x 2 i,n > 0) 
2ir 5 a, 



f X21 rlv ( 

J — c- 2 ^|e 1 (2 M (x 2 i - v)) - E 1 (2 M (x 2 i + v)) 
+ r *L e - 2 ""{Ei(2/i(t; - x 21 )) - Ei(2 M (x 21 + v))\ 

JX21 V L J 



(2.64) 



The integrals in the expression above are not expressible in terms of common analytic functions, 
thus one has to evaluate them numerically using, as before, the GNU Scientific Library. 

The total potential energy V = Vy + Vc + Vq in the units of ma 2 as a function of r = x 2 \ ma g 
for the case when /j, > 0, suppressing two Yukawa singularities, is 

^(xi, 2 ,i, 2 ,A< > 0) = V(xi2,fj, > 0) = ma 2 < 



-rM 



2 ki ln(3) 



rM 



„ Ei (rM) + — Ei (3 rM) 
rM rM ' rM v ' 



r 



J ~ ^ | El ( 2 ( rM - w ^ ~ E i( 2 ( rM + w )) j 

+ / — e ~ 2w \ Ex(2(w - rM)) - E : (2(rM + «;)) 
Am w I 



(2.65) 



where w = v a is a dimensionless integration variable, M = — — — is the dimensionless mass 

ma g 
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4 7T^ G! CK 

parameter K\ = ^ and k 2 — — ^7. Two plots of equation (2.65) are shown in Figure 2.8 

ma* ma* 

for the parameter values K\ = 0.01 and n 2 = 0.1,0.005. The choice k 2 = 0.1 yields a potential 
which is repulsive for small separations and possesses a trough. The trough's depth increases 
with increasing values of the mediating held mass parameter M. For large separations and both 
values of n 2 , the potential approaches zero without crossing the zero of energy. Both choices of 
the parameter k 2 do not exhibit any confining features in the inter-particle potential; although 
both can support bound and scattering states. 

2.6 Improved Particle- Antiparticle State 

It has been shown in the previous sections that the non-linear terms of the Hamiltonian (2.11) 
alter the inter-particle potential. The simple particle-antiparticle trial state (2.27) is incapable 
of probing these terms. Hence, to observe the effects of the non-linear terms on the particle- 
antiparticle system, one must consider a multi-component trial state, such as, 

|*t> = C F |tf 2 ) + % |*4), (2.66) 

where the quantities Cf and Cq are linear coefficients specifying the relative size of each basis 
component, and 

|* 2 ) = J dpi, 2 F(p li2 ) a t( Pl )6t(p 2 )| ), (2.67) 
1*4) = j d Pl . A G(pi... 4 ) a f (pi) fo f (p 2 ) a f (p 3 ) & f (p 4 )|0), (2.68) 

are the adjustable functions containing variational parameters which describe the one-pair and 
two-pair components of the system respectively. The factor of 1/ \J\ is inserted for convenience to 
account for the identical particles in the trial state. Note that no such numerical factor is necessary 
with mono-component trial state. The linear coefficients must be determined variationally, along 
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Figure 2.8: The inter-particle potential V(x 2 i,n > 0) of equation (2.65) in units of ma 2 g as a 
function of r = x 2 i ma g for the indicated values of M, K\ — 0.01, and k 2 = 0.1 (top) and 0.005 
(bottom). 
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with other parameters included in F and G. The linear coefficients are not independent because 
of the normalization of |\&t). They obey the condition 

\Cf\ 2 + ^f- = 1- (2-69) 

The improved particle-antiparticle trial state (2.66) is flexible enough to describe the transition 
from a one-pair to a two-pair state, i.e. the so-called string breaking effect. However, it is not the 
purpose of this research to investigate this phenomenon. 

The matrix element (^l H—E\^f t ) consists of many contributions (see section 5.3 of Appendix 
A for details). Be it as is, many of them are irrelevant to the calculation of the inter-particle 
potential and hence are not included. Varying the matrix element with respect to the functions 
F* and G* leads to the following system of coupled relativistic equations: 

F(P1, 2 ) (W P1 + U P2 - E) = j dpi >2 y 2l 2(pi,2, Pl,2) nPi,2) 

+ Rj d Pl ... 4 |y2,4(pi...4,Pl,2)+C 2l 4(pi...4,Pl,2)}G(pi... 4 ) (2.70) 

G(P1.. 4 ) =R ~ 1 j dP l. 2 {^,2(pi, 2 ,Pl...4)+C 4 ,2(p li2 ,Pl...4)jF(p / li2 ) 

+ f <?Pi...4 |^4,4(Pi...4,Pl...4)+C4,4(Pl...4,Pl...4) + Q4,4(Pi...4,Pl...4)|G( P / 1 ... 4 ) (2.71) 
C G 

where the quantity R = —— is the ratio that specifies the relative contribution of each basis state. 
Cf 

The quantity R is a complex number and, in general, its value is determined via a variational 
calculation of the energy. 

Equations (2.70) and (2.71) are relativistic coupled equations, in momentum space, that de- 
scribe the processes involving one pair and two pair bound state systems but not the processes 
that involve the emission or absorption of the mediating field quanta. Even though one can, in 
principle, calculate the energy spectrum from these coupled integral equations, it is more practical 

to solve them approximately using the variational method which is based on the evaluation of the 
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matrix element of the Hamiltonian. The relativistic kinematics in these equations are described 
without approximation. On the other hand, the relativistic interactions are approximated by the 
kernels 3^2,2, 3^2,4, 34,2 and 3 7 4,4 arising from the Yukawa (i.e. linear) term of the Hamiltonian 
density (2.13), £2,4, C4.2 and C 4j4 from the cubic term (2.14), and Q44 from the quartic term 
(2.15). As explained above, these kernels are adequate approximations when the interaction cou- 
plings are weak. The y type kernels correspond to one mediating field quanta exchange and 
virtual-annihilation interactions, diagrammatically illustrated in Figure 2.1. The C and Q type 
kernels are three- and four-point interactions, illustrated by Figures 2.2 and 2.7 respectively. 

The primary goal of this dissertation is to study the inter-particle potential, henceforth the 
emphasis is placed on equation (2.70). The non-relativistic limit of equation (2.70) is found and 
Fourier-transformed to coordinate space. The symmetric definitions of Fourier transform is used: 

f(pi,2) = 7^3 / <ki,2 c ipi ' Xl e"*-" F(x li2 ), (2.72) 

1 r 4 

G( P i.. 4 ) - J dx i-4 II e ipi Xi G(xi.. 4 ). (2.73) 
The non-relativistic equation in coordinate space takes on the following form: 



2m 2m 



*2, 2 (X1, 2 ) + .R J rfx 3 ,4{^4(x 1 ..4) + C 2 ,4(xi..4)}^^ 



F(x li2 ), (2.74) 



where e = E — 2 m is non-relativistic energy. The two terms containing the function G are 
multiplied and divided by the function F. In this manner, the coefficients of F in the second 
line of equation (2.74) can be identified as the correction to the Yukawa potential ^2,2- Working 
entirely in the non-relativistic limit implies keeping only the leading terms in the inverse power of 
the parameter m. The higher order terms correspond to relativistic corrections; which shall not 
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be considered in this non-relativistic approximation. The non-relativistic kernels are 

g-/iX21 

*2,2(xi, 2 )= a g —— , (2.75) 

(g-/iXi4 g-jUX24 \ 
+ <5(x 3 - x 4 ), (2.76) 

/g-iqi-(x 4 -xi) e -iq2-(x 4 -x 2 ) 
dqi,2 7 -s <5(x 3 - x 4 ), (2.77) 
(> + ( qi + q 2 ) 2 ) (j** + q?) ( M 2 + q 2) 

where a;jj = |xj — Xj | are the inter-particlc distances, a g — — - is the dimensionless coupling 

constant, and a' = —i — 5 — 5- is the coupling constant with dimensions of mass. Note the slightly 
' 256 7r° m 6 

altered definition, i.e. a n = 3a' v . The analytical expression for F 2 ,2 and ^2,4 can be determined 
using the standard technique of integration in the complex plane. The expression for C 2 ,4 can 
not be simplified in terms of the known common functions and is thus left in this form. It should 
be mentioned that these kernels could have also been obtained by employing only particles, i.e. 
replacing anti-particles by particles, in the trial state (2.66). 

In order to proceed with the derivation of the inter-particle potential the functions F and 
G must be specified. It is difficult to solve equations (2.70) and (2.71) for these functions even 
approximately. Instead, in the spirit of the variational method, one can specify ansatze for them. 

The trustworthiness of any variational calculation rests on the choice of ansatze. It is best to 
require for an ansatz to be an eigenstate of the momentum operator, since then, in the centre of 
mass reference frame, an energy calculation corresponds to the binding energy of the system. This 
condition can be fulfilled if there is dependence on the coordinate differences only. At sufficiently 
low energies, the variables x 3 and X4, describing virtual particles and as such, are integrated out. 
It is then prudent to consider functions where the centre of mass motion of the particles described 
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by the coordinates xi and x 2 factors out: 

F(x 1;2 ) * exp(-*Q. £l + *T) expf-^1. 



3 / 2 ^ 2 7 2a 

^/ \ 1 / (xi+x 2 ) 

G(xi.. 4 ) =- ; r 7? exp zQ 



a 9 / 2 "V 2 

x exp ^ I x i2| + |xi 3 | + |xi 4 | + |x 24 | + |x 24 | + |x 34 | ^ ^ 27g ^ 

where x^ = Xj — Xj are the inter-coordinate vectors, Q is the total momentum of particles at 

xi and x 2 , and a is a variational parameter. The factor of two and the fractional powers of a 

are inserted for convenience. By acting with the momentum operator in coordinate space, one 

can verify that both of these functions, and hence the trial state (2.66), are eigenstates of the 

momentum operator with the eigenvalue of Q. Note that the exponentials c~l Xii l/ 2a are suitable 

for the ground state of the system, and this is the state that is considered here. 

Substituting equations (2.78) and (2.79) into (2.74), and after some extensive manipulations 

(see section 5.3 of Appendix A for details), the total inter-particle potential in units of ma 2 as a 

function of length in units of the Bohr radius — - — becomes 

ma g 

V{x 12 ) = Y 2a (x 12 ) + V 2 Y 4 (x 12 ) + y 2 c 4 (a; 12 ), (2.80) 

where, the contributions to the inter-particle potential emerging from the terms containing the 
kernels Y 22 , Y 2i and C 2j4 are 

Q-rM 

Y 2>2 (xi. 2 ) = - ma 2 , (2.81) 



V 2 Y 4 (x 12 )=R J dx 3A Y 2A (xi. A ) 



G(x 



l.A 



F(xi, 2 ) 



o n -r(MA+l)/A 

= - " < */m 2 m + 2) 2 ( A " ^ + 2 r M ^ M - 2 ^ M + 2) ■ (2 ' 82) 
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dx 3 ,4 C 2 ,4(xi..4) 



G(xi.. 4 ) 
^(xi, 2 ) 



2 7r Rum a 2 




-M 2 (Pi+fa+Ps) 



A 3 



x exp 



Pi r 2 + w 2 + /3 3 (w + r) 2 _ |w| + |w + r| 

4/3i23 A 



) 



(2.83) 



where k = — and /3i23 = /?i/3 2 + /?i/?3 + /3 2 /?3- In the above, the values of i? and A must be 
determined variationally. To do such a calculation is a numerical endeavour which is not worth 
undertaking in this scalar model. Instead, one can reason out some acceptable values for them. 
Note that the parameter R appears only as a multiplicative factor in front of the contributions to 
the potential V^ 4 and V^ 4 . In the non-relativistic regime, and with the perturbative values for the 
coupling constants, a reasonable value of R should be on the order of less than 0.1, i.e. the virtual 
component contributes less than 10% to the bound state system of a particle-antiparticle. The 
value of the dimensionless scale parameter A is the physical size on the system. In the absence 
of the non-linear terms of the Hamiltonian, it equals to unity for the ground state. With the 
perturbative coupling constants one should not expect it to deviate away from unity much. A 
small range of values of the parameter A is considered in Figures 2.9 and 2.10, but it should be 
emphasized that the best value is believed to be in the neighbourhood of unity. 

From results as seen in Figures 2.9 and 2.10, one observes that the non-linear interaction 
terms of the Lagrangian density (2.9) provide a slight modification of the Yukawa potential which 
supports bound states only for e = E — 2 m < 0. However, for M > they do not modify the 
potentials to a confining or a quasi-confining form which permits bound states with energies e > 
such as observed in QCD. 



46 



Figure 2.9: The inter-particle potential V, equation (2.80), as a function of distance in units 
of the Bohr radius with R = M = 0.1 and the dimensionless constant n = 0.1. The different 
curves demonstrate how the various contributions alter the potential for the given choice of the 
parameters. 
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Figure 2.10: The inter-particle potential V, equation (2.80), as a function of distance in units of 
the Bohr radius for R = M = 0.1, k = 0.1, and the indicated values of A. 



3 Quantum Chromodynamics 



In physics, you don't have to go around 
making trouble for yourself - nature does 
it for you. 

Frank Wilczek 

This chapter commences with a review of the Dirac equation and the quantization of spinors. 
It is prudent to proceed towards the reformulation of QCD in increments rather than attempting 
to do everything at once. Once the reformulation is complete the inter-quark potential of QCD 
bound states is investigated. 

3.1 QED Quantization 

The starting point of any QFT is the Lagrangian density. The free spinor Lagrangian density 

is 

£(x)=-ip(x)0il>(x)-mi>(x)il)(x) (3.1) 

where m is the bare mass parameters, the barred notation stands for -tjj = ^ 70; and the Dirac 
slashed notation stands for = <9 7° + V • 7. The gamma matrices, in the representation used in 
this dissertation, arc defined as 



I 




<7i 






-/ 




-<Ti 
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where I is the 2x2 identify matrix and at are the Pauli spin matrices defined by 








1 




-i 




1 


















1 







i 




-1 



(3.3) 

Note that other representations of the Dirac matrices are connected by a similarity transforma- 
tion, [42], 

^ = S^S-\ (3.4) 

where S is a unitary matrix. 

The free Dirac equation, obtained by using the Euler-Lagrange equations, is 

(i$ - m) ip(x) = 0. (3.5) 

It has four independent four-component solutions denoted by u s and v s , where s can take on 
two possible values (+, — ) designating the up and down spin orientations respectively [43]. In 
momentum representation, the spinors satisfy the following equations: 



(f -m)u s = 0, 
(i/i + m)v s = 0, 
u s {f-m) = 0, 
v s (p + m) = 0, 



(3.6) 
(3.7) 
(3.8) 
(3.9) 



where the barred notation stands for u — 70 and v — v' 70 and the slashed notation stands for 
= E-f° — p • 7. The explicit form of the spinors is 



u(p, s) = N p 



1 

v ■ P 

w p + m 



Xs, v(p,s) = 7V p 



q- • p 

cj p + m 
1 



where cj p = \]m 2 + p 2 , the normalization constant 



is 



uj p + m 



2m 



(3.10) 



(3.11) 
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and, the two component spin bases are 



X+ = 



1 














-1 




> X- = 




, v+ = 




, V- = 









1 




1 








(3.12) 



the plus and minus subscripts designate the spin up and down projections of particles and anti- 
particles respectively. A way to check this assertion is to apply the spin angular momentum 
operator on the one-particle free states. 

The spinors obey the following orthogonality conditions: 



u t (p,s)u(p,a) = vT(p,s)u(p,cr) = -S-8 aa , 

m 

u\p, s) v(-p, a) = v\p, s) u(-p, a) = 0, 



(3.13) 
(3.14) 



which one can easily verify from the explicit equations (3.10). The completeness relations, with 
the normalization given by equation (3.11), is 



u a (p,s)u fj (p,s) 



2m 



^2 ^ ce (p,s)w /3 (p,s) 



2m 



(3.15) 



afi s =± 

which, again, can be verified from equations (3.10). 

The most general solution of the free Dirac equation is the superposition 

lM*) = £ /(^(^) 1/2 U(p)u a (p, S )e-^ + bi(p)v a (p, S )e^ 

V4(*) = E /(2^(^) 1/2 U(P)uUP,s)e lp - x + bs(p)vl(p, S )e-^ 



a/3 



(3.16) 
(3.17) 



where it is well known that the choices of the coefficients a s (p), a\(p), b s (p) and b\(p) in equations 
(3.16) and (3.17) lead to the correct spin-statistics and positive definite energy in the Dirac 
theory [38]. The spinor index a is explicitly included in these equations. 

The canonical quantization is performed by imposing the following equal time anti-commutation 
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relations on the Dirac fields: 

{Va(x),^(x')} = <5(x-x')<U (3.18) 
tya(x), V6(x')} = {^(x),^(x')} = o (3.19) 
The coefficients written in terms of the fields are 

1/2 

a *(P) = / 7^72" ( — ) e'^u^p,*)^^), (3.20) 



a t(P) = I 77^77 ( — I' 7 ' e-^(x)u a (p, S ), (3.21) 





/ m 


(2tt) 3 / 2 


V W P 


dx 


/ to 


(2tt) 3 / 2 


V W P 


rfx 


/ to 


(2^)3/2 


V W P 


<ix 


/ to 


(2tt) 3 / 2 





1/2 



(3.22) 



1/2 

e^t(p,s)^ a (z). (3.23) 

Again, the spinor index have been explicitly shown for clarity. From equations (3.20)-(3.23), it 
follows that the anti-commutation relations for the coefficients are: 

{a.(p),4(q)} = {6 8 (p),6+(q)}= J(p-q)<y 8<J) (3.24) 

and all other anti-commutators vanish. As operators the coefficients a\ and b\ are identified as 
the particle/anti-particle creation operators respectively, whereas a s and b s are the particle/anti- 
particlc annihilation operators obeying a s \0) — b s \0) — 0, where |0) is the vacuum state. 

The free classical Dirac Hamiltonian follows from the Lagrangian density (3.1). It turns out 
to be 

H = j dxH(x)= J dxij(x)(-iV -7 + m)ip{x). (3.25) 

The quantized version of this classical expression is obtained with the use of equations (3.16) and 
(3.17). After some extensive algebra one arrives at 

H ^ = J^ J dp Wp {^ (P) as(p) + ^ (P) 5s(p) }' (3 ' 26) 
where the Hamiltonian has been normal-ordered to remove the infinite energy of the vacuum. 
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3.2 QCD Reformulation 

The QCD Lagrangian density (units: h = c = 1) [2,38,39], suppressing the spinor and flavour 
indices, is 

C= — {F^f-^A^ + ^ilp^-m)^ (3.27) 

where, the shorthand definitions are 

F% = d, Al -d v Al + g f abc 4 At, (3.28) 
{D li ) ij =d li 6 ij -igAlT? j . (3.29) 

The quark field ipi is a Dirac spinor where the index i = 1, 2, 3 is a colour index. The vector boson 
field A^ represents gluons carrying a Lorentz index \x = 0, 3 and a colour index of the adjoint 
representation a = 1,...,8. The dimensionless coupling constant g characterizes the strength of 
the strong interaction. F£ v is the non-Abelian field strength tensor, while is the covariant 
derivative which is the first ingredient in building a local non-Abelian gauge symmetry. The eight 
SU(3) group generators T?- and the structure constants f abc obey the commutation relation: 

[T a ,T h ]=if abc T c , (3.30) 

where f abc are completely antisymmetric. The term containing £ is known as the gauge fixing 

term. Physical observables do not, in principle, depend on the value of £. However, the canonical 

quantization of the gauge field is impossible in its absence since the conjugate momentum is 
dC 

undefined, i.e. 7r° = — — = if £ = oo. 
dA 

In the perturbative S-matrix formalism one has to include ghost fields. The ghost fields are 

non-physical auxiliary fields required to ensure the unitarity of the S-matrix and appear in the 

beyond-leading-order Feynman diagrams. They are irrelevant for the calculations presented in 

this dissertation, therefore they have been omitted. 
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The QCD Lagrangian density is invariant under the simultaneous transformations 

^"►$ = tf/>* ( 3 - 31 ) 
A-Ti^A'-Ti^UiA-T^ 1 [U-X+ 1 - [d,U\> k [CT 1 ]* , (3.32) 

n» T L -+ ^ T « - ^ F % T « \ u ~ x \ u • ( 3 - 33 ) 

where all the indices are explicitly exhibited. The gauge transformation matrix is 

[/ = e -is|T i c« i W ] (3 3 4 ) 

where the colour indices have been suppressed. The cti(x) are the eight independent rotation 
"angles" . 

It is convenient to express the Lagrangian density (3.27) in the following form 

£ = £a + dp + C^a + £-3A + £iAi (3.35) 

where 

£ A = -\{d»Al - d v A a ^A va - d»An - -L(0" A1)\ (3.36) 

C,^=^}{i$-m)ii, (3.37) 

L^ >A = g ^ i A a T^ h (3.38) 

C 3A = -gf abc (d^)A^A" c , (3.39) 

Ua = -^g 2 {f abc A , lAl){f ade A tld A ,Je ). (3.40) 

In this form, each term of the Lagrangian density is identified by the type of the interaction it 
represents. In the above contributions to the QCD Lagrangian density, the spinor index has been 
suppressed. 

The first step in reformulating the Lagrangian density (3.35) is to obtain the solution to the 
classical equation of motion for the gauge field A^. Having the solution for the non- interacting (i.e. 
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with g — 0) theory, one can form the "formal solution" for the gauge field A 1 ^ in the interacting 
case (i.e. with g ^ 0) similar to how it was done in the Higgsdike scalar theory. Integrating La 
by parts and discarding the total derivative term, since its contributes nothing to the action, one 
obtains 

= \A% (d 2 g^ - d»d v ) A a v + ^A^ff'Al, (3.41) 

d 2 

where d 2 — — V 2 is the d'Alcmbcrt operator and g^ v is the Minkowski metric with the 
signature [1,-1,-1,-1]. Consequently, the least action principle, i.e. S fcdx = 0, leads to the 
following equation of motion 

^V" - (1 - j) d^d^j Al = 0. (3.42) 

This equation has a Green function solution. The Green function (or the free propagator) satisfies 
the identical equation but with a delta function source on the right hand side: 

(flV - (1 - |) A a v b p (x -y) = S ab 6(x - y)8% (3.43) 

or, in momentum space, 

-*V + (1 - j) —J Af p (k) = 6 ab 6». (3.44) 

It is customary to solve for the momentum space propagator and then relate it to coordinate 
space. The momentum space propagator is given by 

A» = -^(^-(l-0^), (3-45) 

where the inclusion of e instructs to use the Feynman prescription to handle singularities. One 
can easily check that 

A^(fc)] _1 A^ bc (fc)=^, (3.46) 
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where A 1 is the operator acting on A in equation (3.44). The momentum and coordinate space 
propagators are related by a Fourier transform: 

A%(x -y) = J -0^ e -*-(*-v) A«t(fe), A«t(fc) = J d 4 xe+ ik <^ A%(x - y). (3.47) 

Hence, the coordinate space Green function is 

A- (, " ») = / (0 { " ^5 (*- " (1 - 0^) } ■ (3.48) 

Returning to the full Lagrangian density (3.35), one finds the complete classical equation of 
motion for the gauge field A® to be 

(d 2 g" v - (1 - |) d^d^j A a v = p» a (x) (3.49) 

where the "source" of this inhomogeneous equation is 

(f a {x) = - g$i{x) ^ T* i>j(x) + gf abc d v (a^x) A^{x)) 

~ 9 f abc ((d»A ub (x) - d v A^{xfj Al{x) + g 2 f abc f cde A b v {x) A» d {x) A ve {x). (3.50) 

Equation (3.49) can be written in integral from 

A a v {x) = (Ao)^ + J dx' Af„{x - x 1 ) p vb {x'). (3.51) 

This expression is dubbed the "formal solution" to equation (3.49). In actuality, it is the integral 
form of equation (3.49). The homogeneous solution A is irrelevant since no systems with external 
gluons will be considered. Therefore, A n is simply omitted. Everything that has been performed 
so far in reformulating the QCD Lagrangian density is in analogy to that of the Higgs-like scalar 
model. 

To proceed further one must specify a gauge. It is convenient to choose the Feynman gauge 
where £ = 1 since the Green function takes on the following simple appearance: 
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Substituting the formal solution (3.51) into the Lagrangian density (3.35), and after some extensive 
algebra, one arrives at the reformulated Lagrangian density 

C = dp + C^ A + Cf A + Cf A , (3.53) 

where, the reformulated contributions are 

C A +C* A ^C$ A = -±g 2 $>fT a il> J dx'Af u {x-x')p v \x'), (3.54) 
C 3A -+C* A =-g f ahc J dx' d» A a 4(x - x 1 ) p°\x') 

x J dx"A^ abe (x-x")p e a (x") J dx"'A^ cf {x-x"')p f (x"'), (3.55) 
Ua -± C? A = -\g 2 (V c J dx'A%(x - x') (Ttf) J dx" A%{x - x") p Tj {x")^j 

x (^f ade J dz'A» adk {x-z')p k a {z') J dz"A^ el (x-z")p l (z")^j . (3.56) 

The arrow represents the process of reformulation, the superscript R stands for reformulated 
quantity and the summation of the colour indices is implied, i.e. -07 M T a ^ = ipij^T^ipj. The 
letters a,b... stand for the colour indices of the adjoint representation and the Greek letters are 
Lorentz indices. To obtain this result, one has to work directly with the action S = J dx C rather 
than with just the Lagrangian density. 

Expression (3.53) is a non-local Lagrangian density, i.e. it involves integration over more than 
one space-time coordinate. Note that the free quark term £^ does not get reformulated since 
it contains no gluon fields. Recall that the purpose of reformulation is to eliminate the gluon 
field from the Lagrangian density while preserving its effects through its propagator. However, 
the "source" term in (3.53) implicitly contains the gauge field A a . Therefore, infinitely many 
iterations of the substitution of the formal solution (3.51) must be performed in order to completely 
eliminate Realistically, this is impossible. Instead, one must truncate the "source" term p°^ 
and work to a given order of iteration. In lowest iterative order, the "source" term is truncated 
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to 

pl{x) = -g^)l tl T^ j {x). (3.57) 

where the spinor indices are suppressed. Using this lowest-order truncation in the terms of the 
Lagrangian density (3.53) leads to 

1 f rlr' dh p-ik-(x-x') 

= + 2 a 2 $(x) rT a i>(x) j - 4,(x>) ^T^(x'), (3.58) 

/rlhrlnrln p-ik-(x-x') -iq-(x-x") p -ip-(x~x"') 
dx > dx » *r fL__ 

x j^Or'h^W)} {fc^0O7^V0O} {^"Ol^W")}, (3-59) 



e -ip-(x-z') e —il-(x—z") 



(2tt) 16 fc 2 9 2 



p 2 I 2 

x \ ^j(z') ^T d 2jj(z') \ \ i;(z") YT e ^(z") \. (3.60) 



The reformulated Lagrangian density contains only quark fields; the interactions involving the 
gluon field are represented by the gluon propagator. Notice, also, that C R contains terms up-to 
order g 6 , hence an energy calculation is expected to be accurate to this order in the coupling 
constant. As in the Higgs-like scalar model, one can examine the effects of the interaction terms, 
including the non-linear terms, in the context of bound states of quarks by considering trial states 
with quark quanta only. 

Finally, the reformulated Hamiltonian density corresponding to the Lagrangian density (3.53) 
follows from the usual consideration 

Hb = U^ + ^U^~C r , (3.61) 

where, the conjugate momenta are defined in the usual way: 

II^ = ^=i^ 7 ° II^ = ^ = -z 7 V (3.62) 
dip r dtp 
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Expressing the fields in terms of the operators, integrating out the spatial coordinates and normal- 
ordering the ladder operators to remove the infinite vacuum energy, one obtains the Hamiltonian 

H R = J ' dx : U(x) : = - j dx : (c* A + £* A + Cf A ) : = H i> +H^ A + H* A + H* A (3.63) 

where, the free Hamiltonian equals 

ff v> = EE /dp"p{<i(p)«-.i(p)+<i(p)6-,i(p)}. ( 3 - 64 ) 

Notice that it is practically identical to the one given in equation (3.26); there is an extra sum- 
mation of the colour index. The quantization of the quark field follows the same procedure as for 
colourless Dirac spinors described in the previous section. The principal difference is exactly in 
this index. Due to the inclusion of this index, the anti-commutation relations have to be slightly 
modified by supplementing a Kronecker delta function in the colour index: 

{as,»(p).<4,j(q)} = {6 s ,i(p),6^ (q)} = 6(p - q) 5 sa 5ij. (3.65) 

This generalization is straightforward and there should be not obscureness if this discussion is 
skipped. 

Unfortunately, there arc no simple expressions for the interactions terms H^ A , H^ A and H^ A 
when they are written in terms of the creation and annihilation operators (3.20)-(3.23) (now 
containing an extra colour index). They are non-trivial and will not be written out explicitly. 
Instead, the matrix elements in the context of a definite trial state will be provided. 

As before, to obtain relativistic equations for stationary states (i.e. time- independent) of bound 
state systems of quark one must switch from the interaction picture to the Schrodinger picture by 
means of equation (2.26). Henceforth, all matrix elements are written in the Schrodinger picture. 
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3.3 Particle-Antiparticle State in QED 

As a preamble, the potential of a fermion-antifermion state in QED is investigated to get 
acquainted with the derivation method in the case of Dirac spinor matter fields. The amount 
of algebra in QCD increases drastically, therefore it is wise to demonstrate the method and 
intermediate steps on an almost "trivial" example 

The simplest particle- antiparticle trial state in QED is 

l* 2 >= E /rfPi,2^(pi,2)at( Pl )&t(p 2 )|o), (3.66) 

s,cr=± J 

where ^F s ,<t(Pi,2) are coefficient functions containing adjustable parameters and s and a are the 
spin indices of the particle and anti-particle respectively. There is no colour index here; it will 
be added in the next section. The QED interaction is given by the term £^ A (3.58) except that 
there is no colour related features in it. 

The relativistic coefficient functions ^^(Pi^), in principle, should be selected a priori to 
be the eigenstate of the quantum numbers J PC (J - total angular momentum, P - parity, C 
-charge conjugation) as one desires; see for example the works [44] and [45]. However, for the 
purpose of deriving the non-relativistic inter-particle potential, it is not essential to specify these 
quantum numbers right away. In the end, only spin (index) function must be specified in the non- 
relativistic limit which, as is illustrated below, decouples from the momentum (and hence angular 
momentum) dependence. In some sense, the approach adapted by Terekidi [46] is followed in this 
dissertation. 

To implement the variational method the matrix element corresponding to the trial state (3.66) 
is calculated and then varied with respect to the coefficient function F* . Carrying out these two 
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operations leads to the following relativistic equation in momentum space: 



J- 5liCTl ( Pl>2 )( Wpi+Wp2 -£;) = E / ^2^?(Pi,2,Pi, 2 )^ 2 , CT2 (Pi,2), (3-67) 

S2,C"2 

where, the QED interaction kernel is 

V*».°Vn i»' ^ V 2 ™ 2 ^(Pl + P2 - Pi - P 2 ) 

Y si cri lPl,2i Pi 2J — o/o x-i ; 

2(2tt) 3 v /ZJp7ZJp^ p 7ZJ p 7 
f «(P 2 > ^2)7^(P2, ^(Pi, si)7^(pj, 52) u(pi, si)7^it(pj, a 2 ) ^(p 2 , ^2)7^^2, Qi) 

\ (K-pi) 2 (P2-P2) 2 

_ u(Pi, s 1 )j f "v(p 2 , <ti) v(p' 2 , a 2 )l f j,u(p' 1 , s 2 ) _ v(p' 2 ,a 2 )l tJ 'u(p' 1 ,S2) u(p 1 ,s 1 )-f l _ l v(p 2 ,a 1 ) \ 

(p[+p' 2 ) 2 (P1+P2) 2 y 

(3.68) 

As in the Higgs-like scalar model, the first two terms correspond to a mediating field quanta 
(i.e. photon) exchange and the last two correspond to a virtual annihilation interaction. This 
expression is nearly identical to equation (2.29) with fi = 0; where the exception is the spinor 
products in the numerators. 

The derivation of the inter-particle potential is continued in the centre of mass frame such that 
the wavefunction takes on the form ^^(pi^) = i^, CT (pi) <5(pi + p 2 ). The equation simplifies to 

F Sl , CTl ( Pl ) (2w p -E) = £ / dp'Y s ^(p,p')F S2!a2 (p% (3.69) 

S2,<?2 

where, the interaction kernel becomes 

ys2,a 2 / n is _ am 2 1 J u(p, Sl )j^u(p', s 2 ) v(-p', ^h^vj-p, Qi) 
.unW'P) (27r) 3 WpWp , I {p l_ p) 2 

-u(p,s l y^v(-p,a 1 )v(-p',a 2 ) 1 ^u(p',s 2 ) (^~~2 4^2") |' ( 3J0 ) 

where like terms have been collected, Once again, this is basically a repetition of the Higgsdike 
scalar model with the exception of the spinor products. Notice that there are spinor arguments 
with negative momentum; they come as a result of the delta function. 
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Unfortunately, and contrary to the customary practice, the products of the Dirac spinors and 
gamma matrices can not be worked out using the completeness relations (3.15) and the trace 
theorems of the gamma matrices. Therefore, the only option to proceed with the derivation of 
the inter-particle potential from the kernel is by using the brute method of matrix multiplication. 

In the non-relativistic limit (see Appendix B section 6.1 for the identities involving spinors 
and gamma matrices) , the Dirac spinors and gamma matrices multiply out to yield the simplified 
kernel 

V s 2 ,<T 2 ( n n '\ _ 9 2 1 ^1 «2 ^o-l cr 2 _ Xsi&i r}cri'nl 2 (J i Xs 2 ] /„ 71 s 

W.lP.P)- (27r)3 | (p/_ p )2 4m 2 J' [6 - a) 

Notice how the virtual annihilation term still carries a non-trivial dependence on spin. In addi- 
tion, in the non-relativistic limit, the spin and momentum dependencies of a coefficient function 
factorize, i.e. F Sli(7l (p) = Sl , CTl /(p). This implies that there are four different equations for 
/(p) in each spin configuration. Recall that spin states transform irrcducibly under the sym- 
metry group SU{2) in the fundamental representation. The four spin index functions S i,<xi for 
a particle-antiparticle pair are obtained by working out the product of two fundamental SU{2) 
representation. This is widely known [47] and the result is 



osinglct _ J_ otriplot 



©1,1 = 1 m z = +1 

©1,2 = 02,1 = 4= m * = ° (3 - ?2) 

V2 



02,2 = 1 m z = -1 

where e Sl<Tl is the two index antisymmetric symbol, m z is the spin projection and all the unspecified 
values of Sl ,o-i in the singlet and triplet configurations are zero. 

Having the results of equations (3.71) and (3.72) in command, one can perform a Fourier 
transform of equation (3.69) to coordinate representation. Simultaneously, multiplying by S1 . CT1 
and carrying out all index multiplications, including the index i in equation (3.71), leads to 
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a Schrodinger-like equation for the particle-antiparticle system corresponding to the trial state 
(3.66) 8 : 

-— ip(x) + V(x)ip(x) = Eipix), (3.73) 
m 

where, the potential energy in the singlet and triplet configurations of spin are 

V sinslct (x) = -a g K (3.74) 
V^(s)=-a fl Q-^*(x)). (3.75) 

g 2 

Here a„ = — is the dimensionless coupling constant of QED and x is the particle-antiparticle 

47T 

separation. Notice that the delta function arise naturally in this formalism, and as seen, is only 
present when the spins arc in the triplet configuration. It is, actually, a relativistic correction as 
is evident from the — - factor. If, while taking the non-relativistic limit of equation (3.70), the 
terms in the leading powers of the inverse mass have been kept, then one could have obtained 
spin-orbit, spin-spin and other delta function lowest order relativistic corrections to the potential 
energy. 

3.4 Three-Quark (Baryon) Trial State 

It has been learned from the Higgs-like scalar model that in order to probe the effects of the 
non-linear terms of the Lagrangian density it is necessary for a trial state to have a component 
with at least three particle quanta. It is anticipated that the same should concur in QCD. As it 
is simpler to work with mono-component trial states, a three quark trial state which models a 
baryon with quarks of different masses, is considered next. 

The simplest three quark trial state with each quark of a different flavour (i.e. mass) is given 



'It is evident that a trial state containing more Fock-spacc components leads to coupled equations. 
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by 

l* 3 ) = fi ^E / d Pi--3^,r,t(pi..3)< i , a (pi)4 i . i6 (pi)aJ ;fe;C (p 3 )|0), (3.76) 

ij'.fe s,r,t 

where the indices k — 1..3 refer to colour, s, r, t = ± refer to spin, and a, b 7 c are the dif- 
ferent flavours. The colour index function ft holds the information about the combined colour 
configuration of the quarks. It has been brought in front to accent that the colour dependence 
of the coefficient function F does not couple to the spin and momentum dependencies even rela- 
tivistically (i.e. the colour dynamics is decoupled from angular momentum as can be seen from 
the Lagrangian density). Recall, that the colour index function of a quark transforms in the 
fundamental representation of SU(3). Therefore, the combined colour configuration Q has to be 
obtained by finding the products of three fundamental SU(3) representations. The result is widely 
known [47] and its derivation also appears in section 6.3 of Appendix B. 

The dimensions of the product of three fundamental SU(3) representations breaks down as 
follows: 

3®3®3 = 1®8®8®10, (3.77) 

where the singlet representation corresponds to the physically allowable colour configuration. The 
index function of the singlet is a three index, totally anti-symmetric object with the appropriate 
normalization: 

^i,j,k — ~~/Q €i i k ' (3.78) 

Note that if the trial state (3.76) contained quarks of the same flavour only, then the overall trial 
state would have to be anti-symmetric under interchanges of quark momenta and indices if it is to 
model the ground state. Since the colour index function is always anti-symmetric, consequently 
the remaining (spin and momentum) parts would be symmetric. However, this is not the case here. 
Recall that the coefficient function F contains adjustable parameters which must be calculated 
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variationally via an energy calculation. 

The matrix element of the Hamiltonian (3.63) with the trial state (3.76) is given in section 6.2 
of Appendix B. Variationally differentiating it with respect to the coefficient function F* leads to 
the following relativistic equation in momentum space: 

^,r,t(Pl,2,3) « + < + ^p 3 ~E) = 

E / ''Pi..;.-'. {^''*\pi,2,3 ) Pl,2,3)+<i' t '(p;,2,3 ) Pl,2,3)}^Sr'^(Pl,2,3) I (3-79) 
s',r',t' J > 



where the slightly generalized definition wj, = y mf + p 2 accounts for different flavours. 
The relativistic interaction kernel y in equation (3.79) is 

y.,r,t W^ Pl ^-2j2^{3) 

m A m B <$(Pi + p 2 - Pi - Pa) S(p 3 - P 3 ) <W 2 "(Pi> s hv u (Pi, s')u(p 2 , r)^u{p' 2 , r') 

«« 2 <) 1/2 (P^-P?) 2 
m A m B <?(pi + p 2 Pi - P 2 ) ^(P3 - Pj) Kt 2 u(P2,r)-/^u(p' 2 , r')u(p!, s^ujp^, s') 

(««) 1/2 (Pi 4 -^) 2 
m A m c S(p 1 + p 3 - p^ - p 3 ) (5(p 2 - p' 2 ) (5 rir2 u(pi, s)7 m m(p' 1; s')u(p3, ^7Mp 3 > 
«o£, W g,a£,)V2 (p 3 C -pf) 2 

m^mc (5(pi + p 3 - Pi - p 3 ) £(p 2 - p 2 ) ^rirg m(P3,Q7m m (p^ O^Pj^sjyMpjV) 

««^) 1/2 (p?-p'l A ) 2 

m B m c S(p 2 + p 3 p' 2 P 3 ) ^(Pi - Pi) ^i£2 "(P2, rh^ujp'z, r')u(p 3 , t)^u(p' 3 , t') 

««-£) 1/2 (P^-P'f) 2 
m B m c 6(p 2 + p 3 - p' 2 - P 3 ) <5(pi - Pi) S SlS2 m(P3,^i)7m m (P3: t'MP2,t)-f tl u(p 2 , t') 

(«"&«£) 1/a (pf-P?) 2 

(3.80) 

This kernel gives rise to pairwise Coulomb interactions emerging from one-gluon exchange; there 
are no virtual annihilation terms in the absence of anti-quark quanta in the trial state (3.76). 



+ 



+ 



+ 



+ 



+ 
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The cubic interaction kernel C in equation (3.79) is 



W,r',t' f , „ ,_ id 4 jp_ <KPi +P2 + P3 - Pi - P 2 - P3) 



C s ,r,t' (Pi,2,3.Pi,2,3) = -7^-^ F c m A m B mc 



PI P2 P3 p' 3 ' 



u(pi, s)7^(pi, s')«(P2, r)7"«(P2. r')u(P3, t)7"«(P3. 7ZA _ „m w J? _ Jb^c _ ^2 



(pf-pl A ) 2 (pf -^ B ) 2 (pf-pf) 2 
(pf-p'i 4 ) 



+ «(Pi, s)7^(p'i, s')m(P3, 07 AI "(P 3 ^')"(P2, '")7 I/ w(p 2 , r') 
+ w(P2, r)7„u(p 2 , r')w(pi, s) 7 Al w(pi, s')u(p 3 , f^w^, t') 
+ «(P2, r)7„u(p 2 , r')w(p3, t)^u(p' 3 , i')w(pi, s) 7 ,y u(pi, s') 

+ M(P3,t)7^(P 3 ^ / )^(Pi.g)7^(pi,g / )^(P2,r-)7^(P2^ / ) ^. 1 _ ,/ h2( p B - p 'B)2 fr ,r T^rVi 



(Pi 1 - P'l A ) 2 (pi - P2 B ) 2 (P3 ~ P' 3 C ) 2 

^-Pf)^ 

(pf-p^) 2 (pf-p 2 B ) 2 (p^-pf) 2 



+ m(P3, <)7^w(p 3 , i')w(P2, rhMP 2 , r 'MPi, •s)7 ly «(pi, s )- 



' (Pf- Pi A ) 2 (P 2 B - P 2 B ) 2 (p£ - Pf) 2 

(3.81) 

This is a three-way interaction which, in addition to spin and momentum dynamics, involves 
colour dynamics. The colour dynamics reside in the factor Fc which surprisingly turns out to to 
be vanishing. There are six separate colour factors in each term but all lead to the same result 
(refer to section 6.2 of Appendix B for details). For example, 

jp„ fabc q ,0 , T a T b T 10 — — f abc c c , r F a T 10 T c 

r C — J s '»i,Ji,fei lL i 2 ,3 2 M 1 i 1 i 2 1 hii ± k 1 k 2 ~ g J fc «ljlfcl t *232k 2 1 hi 2 1 jih 1 k x k 2 



— ^ f abc T? li2 T n j 2 T£ ik2 (Si 1 j 2 5j 1 k 2 Sk 1 i 2 + fii 1 k 2 8j 1 i 2 5k 1 j 2 ) 



g J ~lll 2 ]lj2 

1 

g' \~J 2 1 2 ~ K2J2 ~l 2 K 2 ' ~ K 2 l 2 ~ l 2 ] 2 ~ ] 2 K 2 J g 



« f abC (T? 2 i 2 ?U T ^ 2 + TU T» 2j2 T9 aka ) = - r»< Tr (T° T c T b + T a T b T c ) 



= I /«bc Tr ,_ T a T b T c + T a T b T c\ = Q) / 3 82 x 

6 

where the fact that the generating matrices T are traceless, the structure constants f abc are anti- 
symmetric under index interchanges and the following identity involving cpsilon tensors have been 
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used: 







(3.83) 



S kii 2 *fcij 2 S kik 2 



This vanishing contribution to the energy from the cubic term of the Hamiltonian H^ A seems 
to be somewhat mysterious. The quartic term H± A does not contribute anything either; it is a four 
way interaction hence a trial state, such as in equation (3.76), with only three particle/antiparticle 
operators does not to probe it. It is improbable that the non-linear terms do not participate at 
all, especially when the effects of confinement and string breaking are ascribed to them. It might 
be that this surprising and counter-intuitive result only occurs in this leading order and for the 
simplified trial state (3.76). Indeed, reference [11] illustrates the same vanishing result for a system 
of three quarks, yet using a different method. In the next section though, it is shown that a multi- 
component trial state does produce non-vanishing non-linear contributions to the interaction and 
a reason for this is suggested. 

3.5 Multi-Component Quark- Antiquark (Meson) State in QCD 

The following multi-component trial state, consisting of a heavy quark-antiquark pair and a 
light quark-antiquark pair, is treated in this section: 



|v& 4 ) = ]T Kuki J d Pi--^M,^(Pi..4)< i ,„(Pi)6^,„(P2)aJ iM (p 3 )6 I I iIi6 (p 4 )|0), (3.86) 

i,j,k,l K.,A,/i,i/ 
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\*t) = C F \* 2 ) + C G \* 4 ) 



(3.84) 



where the two Fock-space components are 




(3.85) 



where the unsummed indices a, b label the two flavours of quarks, k, A, [i, v are the spin indices, 
and i,j,k,l are the colours indices along with the colour wavefunctions f2 and A for the two and 
four components respectively. As previously noted, the colour index functions factor out in front 
even relativistically. The relativistic coefficient functions F and Q are, in principle, such that 
\%?t) is an eigenstate of the quantum numbers J PC one desires. The coefficients Cf and Cq are 
parameters which, modulo normalization, must be determined variationally together with all the 
parameters contained in T and Q. 

The colour index of a quark (anti-quark) transforms in the fundamental (conjugate) represen- 
tation of SU(3). To find the colour index functions fi and A one must consider group products 
of SU (3) representations. The dimensions of the group products of SU (3) representation can be 
determined using the Young tableaux method as is explained in Appendix B section 6.3. 

The dimensions of the two required decompositions are 



It is known experimentally that the physically allowed representations are colourless (i.e. zero in 
each of the eight colour charges; see below) corresponding to the one dimensional representations. 
The singlet in equation (3.87) provides the colour wavefunction for the two component (3.85); it is 
nothing but the delta function with the appropriate normalization. Similarly, the two singlets in 
(3.88) provide the colour index function for the four component (3.86). The index wavefunction 
is, again, comprised of delta functions and, in this case, has to be equally weighted between the 
two singlets. The properly normalized index functions are 



3 ® 3 = 1 8 8, 



(3.87) 



3<g>3<g>3<g>3=l©l©8©8©8©8©10©10©27. 



(3.88) 



'ij ~ /t; 0i i ' 




(3.89) 




f 



{Sij S M + 5u 5 kj ) . 



(3.90) 



2V6 
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Recall, that the total colour charge operator can be obtained from the colour charge current 
density as follows: 

Q a = J dx^(x) 7 °7-^(x) 

i,j=\..3s=± J ^ J 

One can verify that the trial state (3.84) with the colour index functions given by (3.89) and 
(3.90) is an eigenstatc of Q a = for a = 1, 8. 

The procedure of extracting the quark-antiquark potential is the same as discussed in the 
previous sections. The relevant contributions to the matrix element (ty t \H — E\$F t ) are given in 
section 6.4 of Appendix B. The variation of the matrix element with respect to both coefficient 
functions T* and Q* leads to the coupled equations that describe a bound state of a meson: 

J- Kl , Al ( Pl , 2 ) « + < -E)= ]T J dp 3A {y 2 ,2)7± (Pi--*) -^(P3,4) (3.92) 

+ r y, J rf P3-e { w^iT'" 2 (pi..e) + {C2aT:±^ V2 fa.*)} g«,m,^<j*..*)i 

Gm.a^.^Pi^) -E)= (3.93) 

+ ^ E / d P^ { {y^)Z'Z^ Vl (Pi.,6) + (Cd.aCS,^,^ (Pi..e)} ^«CP6,a) 

+ E / { (^Ca^;" (P1..8) + (Si4):;^;: (pi-8)} ^2,A2,,2,, 2 (P4..8), 

where the colour indices have been summed over and the resulting colour factors included in 
C G 

the kernels, and R — — — is the ratio that specifies relative contributions of each Fock-space 
Cf 

component. The kernel C have been deliberately omitted because, just as in (3.82), the colour 
summation, despite the presence of the "spectator" particle which just produces a delta function 
in the colour index, sets this kernel to zero. 
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The relativistic kernels pertaining to the equation for the function F are 



n , x«2,a 2 , v 2 g 2 (5(pi + p 2 -p 3 -p 4 ) 

(J 7 2,2j Kl . Al (Pi. .4) - 



'«i,Ai V»' 1 "^ 3(2?r) 3 /AAA i\V2 



3(2^)3 (o 



x <( u(pi,Ki)7m u (P3)«2) w(p 4 , A 2 )7 M w(p 2 , Ai) ( — -^ + 7 — ^ 1 [>, (3.94) 



2 g m A m B 

0>2A) KlM (P1..6) 3(27r)3 

x / ^(Pl ~ P3 ~ Ps - Pe) ^(P2 - P4) S Xl x 2 
I (oj a uj a uj b u b ) 1/2 

K V Pi P3 Ps Pe/ 



x u(pi,Ki)7^u(p 3 ,K 2 ) v(p 6 ,u 2 )^ v u(p 5 ,fj, 2 ) 

+ 



(P1-P3) 2 {P5+Pe) 2 

S(P2 - Pi- P5 - Pe) <5(Pi - P3) S K1K2 



(««) 1/2 

r(pi, A 2 )7 i ,^(p 2 ,Ai) v{p Q ,v 2 ) :" "'p:..^; ( — ~ ( p5 + p6 )2 ) f ' ( ' ! - !):,) 



(t-2,4j KliAl (Pl..6j - 

fabcfabc g 4 ^ ^ ^ + p 2 - p 3 - p 4 - p 5 - p 6 ) 1 



(2tt) 6 U w ^ w A u A u B u B \ V2 ( pi - p 3 )2 (pjj _ p 4 )2 fa + pg )2 

x j(P3 -pi) Al «(pi,Ki)7f«(P3,K2) w(p 4 , A 2 )7 M u(p 2 , Ai) u(p 6 ,^ 2 )7 l/u (P5,A*2) 

- (P3 - "(Pi, «l) 7!/ w(P3, «2) u(pe, ^2) 7^ U (P5, M2) u(p4, A 2 ) 7" "(P2, Ai) 
+ (P4 -p 2 ) M w(P4, A 2 )7i, v(p 2 , Ai) v(p 6 ,u 2 )-f^u(p 5 ^ 2 ) u(pi,Ki)7 I/ ii(p 3 ,«;2) 

- (jPi -P2) M w(P4, A 2 )7 ! ,w(p 2 , Ai) u(pi,Ki)7 A 'u(p3,K2) C(p 6 ,^2)7 ,yu (P5,A*2) 
+ (P5 +P6)m«(P6,^2)7^ m (P5,A*2) u(pi,Ki)7 Ai w(p 3 ,K 2 ) w(p 4 , A 2 )7 ,y u(p2,Ai) 

- (P5 +P6)m«(P6,^2)7^ m (P5,A*2) w(P4, A 2 )7 AI w(p 2 ,Ai) u(pi , Ki) 7^ u(p 3 , K 2 ) >■ (3.96) 



Equations (3.92) and (3.93) are coupled relativistic integral equations for the wavefunctions T 
and Q which are difficult, if not impossible, to solve. The function T describes a quark- ant iquark 
pair of mass tha, while Q describes a two quark- ant iquark pair state with each quark flavour of 
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mass uiA and tub respectively. These equations are, in principle, capable of describing the process 
of string breaking in QCD, i.e. the transition from a single meson state qa q~A to a two meson 
state qAQAlB q_B- The relativistic kinematics are described fully and without approximations. 
On the other hand, the dynamics are described approximately by the kernels 34,2, 34,4, 34,2, 
34,4, ^2,4, C4.2 and because of the approximate nature of the trial state (3.84) and the first- 
order iterative approximation in the reformulation procedure. The emphasis in this dissertation 
is on the quark- ant iquark bound state described by the function J 7 , therefore it is required to 
calculate the kernels 34,4, £4,2 and Q4.4 for this purpose. Henceforth, equation (3.93) will not 
be investigated. The colour indices have been summed in 34,2 and 34,4 and the results appear 
as multiplicative factors in front. In £2,4, the colour factor is expressed by the contraction of the 
structure constants f abc f abc . Note that the spin and momentum dependencies are still coupled 
in these relativistic equations whereas the colour dependence is separate and calculated fully. 

In the non-relativistic limit, the spin and momentum dependencies of the wavefunctions de- 
couple (i.e. separate into factors). That is, the functions T and Q can be written as the products 

^K,A(pi,2) = e K , A F( Pl , 2 ), (3.97) 
Gk,\,h,v{vi.a) =£k,\,h,vG(pi.a)- (3.98) 

The spin indices n, A, [i and v transform in the SU{2) representation. Therefore, the spin indices 
in and S transform according to a particular product of SU(2) representation. The dimensions 
of these products of representations have been worked out using the Young tableaux method in 
Appendix B section 6.3: 

2®2 = 1©3, (3.99) 
2®2®2®2 = lffilffi3ffi3ffi3ffi5. (3.100) 

The index functions corresponding to these representations, such as 6 and H in equations (3.97) 
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and (3.98), are eigenstates of the spin operator such that, schematically (spin indices must be 
included to be precise), 

S 2 A= s(s + l)A (3.101) 

S 3 A = m 3 A (3.102) 

where A = 9, S is either of the spin index functions (i.e. equation (3.97) and (3.98)), s is the 
total spin quantum number and m 3 is the projection of the spin (i.e. the eigenvalue of the Pauli 
matrix a 3 in the representation given by the equation (3.3)). 

According to equation (3.75), the triplet spin configuration produces a delta function contri- 
bution to the potential energy. Hence, it is expected that an intricate delta function dependence 
might arise if the spin index function is not in the singlet representation. The simplest spin 
configuration to consider, then, is the singlet where the spin eigenvalues are S 2 — S3 = 0. The 
corresponding properly normalized spin index functions are 

e K ,A = ^e K ,A, (3.103) 

S K ,A,/i,i/ = -^-^ (e«A Zy,v + £kv c m a) , (3.104) 

where e is the totally-antisymmetric tensor. One can verify, by explicitly writing out all spin 
indices, that spin index functions (3.103) and (3.104) satisfy the spin eigenvalue equations (3.101) 
and (3.102). 

The non-relativistic equation for the function F in the spin singlet configuration can be ob- 
tained by multiplying the relativistically reduced equation (3.92) by 6, of equation (3.103), and 
summing over all the spin indices (see section 6.4 of Appendix B for some details). Simultane- 
ously, the functions F and G are expressed in the centre of mass frame, such that the total energy 
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corresponds to the rest mass of the system: 

J r (Pi,2) = /(Pi)*(Pi+P2), (3.105) 

G(Pl.. 4 ) = 5(Pl,2,3) *(P1 + P2 + P3 + P 4 ). (3.106) 

Consequently, one can integrate out one momentum dependence using the overall momentum 
conserving delta functions and thus reduce the number of variables in the equations. Such an 
integration leads to slightly altered interaction kernels which exhibit a degree of skewness (loss of 
symmetry) in the dependence on the momentum variables. 

The non-relativistic equation for the quark-antiquark system in the singlet configuration of 
spin and in the centre of mass frame is 

/(Pi) - £ ) = J dP3 y 2 ,2(pi,3) /(P3) + Rj rfp 3 ,4,5 i p i ff(p 3 ,4, 5 ). (3.107) 

where £ = E — 2 is the non-relativistic energy and all spin indices have been summed. The 
non-relativistic kernel 5^,4 has been omitted since its leading contribution is of order 0(m~ 3 ). The 
remaining non-relativistic kernels 3-2,2 and C 2j 4 have the leading contribution in order 0(m~ 2 ): 

y 2 , 2 (p,, s) =i «L( >J ( 1+ <Ei_Pl)! + <Pi±P|)), (3 . 108) 

3 (2tt) 3 y(p 3 - pi) 2 y ^ Am\ 2m 2 A J 

fabcjabc g i x I 13 
^2,4(Pl,3,4,5) — * , fr> n 6 2" 7 727 1 72 ^ 2" Pi ' P3 x P4 

4(27r) 6 m|j (p 3 -pi) 2 (p 4 + Pi) 2 \ 2m A 



3 (mg + toa) 



(pi • P5 x p 3 + pi • p 5 x p 4 + p 5 'p 3 xp 4 jj, (3.109) 



where an extensive simplification has been performed to bring the kernels to the present form. 
The appearance of the factor i in C24 looks troublesome since the inter-particle potential can not 
be imaginary. However, the quantity R is not restricted to be real and one can make a choice 
such that it is purely imaginary to enforce the overall term to be real. One can see from the 
spin and spinor indices summation in Appendix B section 6.4 that at least there is a degree of 
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consistency in this approach since both cross terms and 62,4 have an i in front. Even so, there 
still remains the choice of the phase which leaves the sign in front undetermined. The variational 
method, if implemented, would select the sign corresponding to the lower energy (more on this 
matter is written below). Also, notice that the discernible momentum dependence implies that 
the Fourier transform of equation (3.107) to coordinate space would not separate the inter-particle 
potential and the wavefunctions into separate factors. 

g-ipi-x 

To perform a Fourier transform of equation (3.107) one multiplies by and integrates 

(2ny> z 

over pi. As already mentioned, such an operation does not lead to an equation where the wave- 
function and inter-particle potential stand as separate factors in all terms. In the terms where the 
decoupling does not occur, one is forced to multiply and divide by /(x). Subsequently, an ansatz 
for / must be provided, such as the following choice 

/(^\At cx p(--)' /(p) = -rft-rr^' (3 - 110) 

V 7rer V aJ 7r (p J or + iy 

which is the properly normalized wavefunction for the ground state of hydrogen (or, more precisely, 
positronium) , with a being the characteristic size. 

Similarly, an ansatz for g has to be provided. In the coordinate representation, the function 
g has to be brought to a suitable frame where the particles described by the coordinates xi and 
x 2 must be in the centre of mass frame. The dependence on the coordinates x 3 and X4 can be 
arbitrary since it is integrated out. There is no unique choice for such a choice of coordinates. 
However, the following transformation has been found useful: 



x 3 




1 


1/2 










X 


x 4 




1 


-1/2 










X34 


x 5 




1 


1/2 


-1 







X45 


x 6 




1 


1/2 





-1 




X 5 6 



(3.111) 
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One can carry out a schematic Fourier transform of G using the coordinate transformation (3.111) 
to obtain the following dependence in momentum space: 

■(P3+P4 + P5+P6) 



/ dx3 - 6 ft eWX * G ( x 3..e) = J dXdx 34 dx 45 dx 56 e lX ( 

x e -il/2x 3 4-(p3+P5+P6-p4) e -«4 5 -p5 X56 ' P6 G(x 34 , X 45 , X 56 ) 

= (2 7r) 3 (5(p 3 + p 4 + p 5 + pe) g( P4, Ps, Pe ) (3.112) 

The delta function refers to the centre of mass motion which is a non-normalizable factor. For 
the remaining function g, in equation (3.112), one can make a variational ansatz such as 

a 9/2 (2^23/2^ 1 l__ 

31P4, P5, P6j ~ fabc fabc (p 2 a2 + 1)2 (p 2 fl2 + 1)2 (p 2 fl2 + 1)4 ' 

where the factors in front are included for convenience and the normalization of this wavefunction 
is absorbed into the definition of R (in fact, all constants in front can be absorbed into the 
definition of R). Equation (3.113) reflects, as in equation (3.110), a factorized hydrogen-like 
dependence for each three-momentum variable and is only appropriate for the ground state since 
it contains no angular dependence. A similar ansatz for the four-component coefficient function 
G has been used in paper [27]. 

Upon calculating the Fourier transform of equation (3.107) using the ansatze (3.110) and 
(3.113), the equation in coordinate representation becomes: 

- ^ /(x) + (V^x) + V 2 (x)) /(x) = £/(x). (3.114) 

The V\ contribution to the inter-particle potential is obtained from 



M*)/(x) = 7^372 / dpi,3e- <Pl - x F2,2(pi,3)/(P3) 



4 ^ Ui.e-'H^H* + ' /(P8) 



3 27r 2 7 1 (p 3 -pi) 2 ±m\ 



2^m\ (P3-P1) 2 (P|« 2 + l) 2 
74 



g 2 

where a s = - — is the dimensionless coupling constant of the strong interaction. The last term in 
the second line of equation (3.115) has been multiplied and divided by f(x), and (3.110) used to 
approximate / in coordinate and momentum spaces. It is manifest that the contribution to the 
potential energy from the last term is non-local and the procedure of multiplying and dividing by 
/ makes a local approximation for it. It follows from (3.115) that 

4 I 1 7r „. , a 3 rx\ f , exp(— i pi • x)(p? + po) I 

V 1 (x) = --aA- + —5{yi)+— i — T exp(-) / rfpi, 3 , 2 / Tg , 3.116 

3 [ir m A 4ir 4 m z A W J (p 3 - pi) 2 (pg a 2 + l) 2 J 

where the — term has been obtained by the usual integration in spherical coordinates followed by 

x 

the radial integration in the complex plane. With little effort, one can reduce the integral in the 
last line to a double quadrature: 

Vi{x) = -\m A a 1 s \ - + a 2 s ir5(r) 
3 r 



2 2 / r \ [°° J / r \ f°° pf+pi , 

+ a '^A2- r e <A)J dpiSm [ Pl A J dp3P3 UTTT 



Pi + P3 



(3.117) 



Pi ~P3 

where r = |r| = |x| tua a s is the dimensionless quark-antiquark separation variable, A — am a ot s 
is the dimensionless variational parameter and the remaining integrals are expressed in terms 
of dimensionless momentum variables. In light of the fact that the second and third terms in 
equation (3.117) are of the same order in a s , it is conceivable that the result of the intricate 
double integral of the third term also leads to a delta function dependence as in the second term. 
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The V2 contribution to the inter-particle potential comes from 

-V 2 {x) /(x) = - / G?pi j3i 4, 5 e^ P1 - X C 2 ,4(P1,3,4,5)5(P3,4,5) 



(2tt) 3 / 2 



a 2 ~ 6 



?-e-/»/(x) / 



i?z^e^ a /(x) / rfpi,3,4, 5 e- 



«pi-x 



1 



X 



>3-Pi) 2 (P4 + Pi) 2 (p|a 2 + l) 2 (p|a2 + l)2 ( P 2 a 2 + i) 4 y 

13 3 (m B + m A ) ( \\ 

J" Pi ' P3 X P4 o Pi ' P5 X p 3 + Pi • p 5 X p 4 + p 5 • p 3 X p 4 



(3.118) 

where, again, it has been multiplied and divided by / and the ansatze (3.110) and (3.113) have 
been substituted. Then, it follows that the V2 contribution to the potential is 



y 2 (x)= T \R>\^ e r/ A Jd Pl , 3A , 5 e 



i pi -r/A 

1111 



(P3-P1) 2 (p 4 + Pi) 2 (p| + i) 2 (p! + i) 2 (pg + 1) 4 ; 

(pi • P5 x p 3 + pi • p 5 x p 4 + p 5 • p 3 x p 4 H . (3.119) 



/ 13 3(£+l) 
x hr Pi • P3 x p 4 - 



2 1 1 1 ■> 1 1 ^2 

where the notation raj = has been introduced and the necessary substitutions have been 

made as before to make the integration variables dimensionless. The redefinition R' — e J7r / 2 R = 
i R does not provide the overall sign as there are two choices of the phase which yield a real 
number. The correct sign will be determined below. 

In the non-relativistic limit, the heavier quark flavour is suited to be either the charm with 
the mass of 1.25 ± 0.09 GcV or the bottom quark with the mass of 4.70 ± 0.07 GcV [2]. On the 
other hand, to model a realistic string breaking effect, the lighter quark flavour could be the up 
or the down quark bearing approximately identical but ill-defined masses in the neighbourhood 
of 1.5 — 7.0 McV [2]. Hence, the value of £ = 0.001 is a decent benchmark and approximately 
corresponds to the mass difference of the proposed flavours. With £ = 0.001 the first term in 
equation (3.119) can be ignored. 
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The numerical value of the coupling constant a s must be obtained from equation (1.4). It 
should correspond to the bound state energy of, say, a bottom quark-antiquark system in the 
singlet spin configuration which is 9.86 GeV in energy [2]. This is the QCD bound state whose 
potential will be modelled. As an input to equation (1.4), one can use the experimentally measured 
value a s (mz) = 0.117 where mz = 91.19 GeV [2]. Upon substitution, one retrieves a s (9.86) = 
0.167. When expressing the potential in units of ma 2 s , the contribution Vi is of order a 5 s which 
justifies the usage of such a large value for the coupling constant. The bottom quark-antiquark 
system, one unit of distance in Bohr radii corresponds to approximately 2.51 x 10~ 16 m of actual 
length. 

Equation (3.119), as it stands, contains two, yet undetermined, variational parameters A and 
R (or R'). The numerical values of these parameters should be obtained from a variational 
calculation of the matrix element ($> t \ H — E \ ^>t) ■ Such a calculation requires great effort which 
is not undertaken in this work. Instead, the approach shall be to multiply equation (3.107) by 
/*(pi) and, thereafter, integrate out the free variable pi: 



This is an expression for the total non-relativistic energy of a meson in the singlet spin config- 
uration and it could be used to estimate one of the two variational parameters. A variational 
calculation of energy, according to a theorem [25], always yields a result which is greater than 
the true value. As QCD bound states are confining (i.e. £ > 0), one can obtain the following 
expression containing the parameters A and R (see Appendix B section 6.5): 




R / dpi,3,4,8 f (Pi) C 2 ,4(P1,3,4,5) <7(P3,4, 5 )- (3.120) 




a 



■g(l + 16 7T- 2 Cl) 
A 3 



) 



) 



±\R'\ 



C 2 



(3.121) 
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where the constant C\ w 1.85044 emerges from the calculation of the first term on the left hand 
side of equation (3.120) (this calculation appears in Appendix B section 6.5). The constant Ci'- 
A 5 



A b f 

C2 = T~- I d PL3,4,5 /(Pl) C 2 ,4(Pl,3,4,5) 0(P3,4,5) 

m A Q s ' J 

_8_ 



dpi, 



3,4,5 



(p? + l) 2 ( P § + l) 2 (pi + l) 2 (p 2 + l) 2 ( P3 - Pl) 2 (P4 + Pl) 2 



1 13 3(e + i) 

x I — Pi • p 3 x p 4 ^ 



(pi • P5 x p 3 + pi • p 5 x p 4 + p 5 • p 3 x p 4 ^ (3.122) 



2 ^ L ro r e 

is a multidimensional integral expression which one has to solve numerically using, in practice, 
the Monte Carlo method. 

Unfortunately, it turns out that Monte Carlo integration of equation (3.122) does not produce 
reliable results. The troublesome pieces of the integrand are the triple scalar vector products in 
the numerator. To circumvent this difficulty, one can place an upper bound on the integral by 
considering upper bounds on the scalar vector products: 

p • q x k = |p| |q| |k| sm9 cos<£ < |p| |q| |k|, (3.123) 

where 4> is the angle between the vectors q and k and 6 is the angle between the vectors p and 
q x k. 

When all scalar vector products are replaced by their upper bound estimates, equation (3.122) 
can be reduced to a triple quadrature in the radial coordinates of the momentum variables: 
24(^+1) 4tt 4 



Co < 



4 



4tt 4 f 

— / dp 1 dp 3 dp A p 3 p A (p 1 p 3 +P1P4+ p 3 Pi) 



1 \J<»+nr+<*\j<*+»r+s\ (3 . 124) 



where only the leading terms in the parameter £ have been kept. To obtain this result, angular 
integrations resembling those in the calculation of the constant C\ have been employed. The 

78 



10- 2 
io- 4 

10 6 
10 s 
IO" 10 

10- 12 
10- 14 

10" 16 
10" 18 
IO" 20 




































































































































































































10" 3 10" 2 10" 1 1 10 



A 

Figure 3.1: The allowed values for the parameter A and R' according to equation (3.121) with 
the plus sign. 

parameter a; is a regulator which has been inserted to ensure that a numerical evaluation of this 
triple quadrature converges. In general, one must perform a handful of numerical integrations 
where the value of u is reduced in each successive trial. If the integral converges then its numerical 
evaluation should approach a constant value as lo tends zero. In this case, the triple quadrature 
is found to converge to the value 

C 2 < 1.005 x 10 16 . (3.125) 

Turning attention back to equation (3.121) itself, the energy of the bottom quark- antiquark 
state £, expressed in units ofmOj, is 3.55. Accordingly, one can solve the equation to determine 
the values of A and R which could produce such energy. Doing so, one finds that only the plus 
sign in equation (3.121) (correspondingly the negative sign in equation (3.119)) is realizable for 
this energy. Figure (3.1) shows a plot of the allowed values of the parameters A and R' . Recall, 
the optimal value of A is believed not to deviate much from unity. 
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Having figured out all the values for the inputs, it is time to return to the expression for 
V2, equation (3.119). The multidimensional expression can be reduced to the following triple 
quadrature: 



where similar steps have been taking as those in the calculation of the constant C\ . The parameter 
a; is a regulator which is, ultimately, to be taken to zero and, as before, R' — i R is a redefinition 
with the now-determined sign. The correct sign, according to equation (3.126), must be negative. 
For an instructive purpose, inter-particle potential curves corresponding to both signs in equation 
(3.126) are plotted in Figure (3.2). The correct sign (-) exhibits a confining potential at larger 
separation distances, whereas the wrong sign (+) yields a potential which does not support stable 
bound states. 

It is useful to examine the behaviour of the potential for somewhat different values of the 
parameters A and R. To this end, Figure 3.3 illustrates the quark- ant iquark potential for the 
indicated values of the parameters A and R which are obtained from the curve in Figure 3.1. The 
quark-antiquark potential is substantially altered from its Coulombic behaviour by the H^ A term 
of the QCD Hamiltonian (3.63). The apparent linear segment is characteristic of all three curves 
(at least in the domain r < 7 beyond which the numerical results become unreliable). The higher 
values of A seems to render the linear segment longer. Beyond the separation distances shown on 
the graph, the points become increasingly scattered. The likely cause of this scattering is in the 
difficulty of evaluating accurately the product of the exponential factor e r / A and the numerical 
calculation of the triple quadrature in V2, equation (3.126). This scattering can be diminished 



Va(x) 



T\R'\ 



m A a 7 s 4tt 4 (£ + 1) e r / A 




(3.126) 
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by requiring greater numerical accuracy in the evaluation of the triple quadrature but, of course, 
at the expanse of computational time. In any case, the validity of the curves is questionable 
beyond the linear segment and this is likely due to the limited accuracy in evaluating the triple 
quadrature. This is further discussed in the final section. 

It is of interest to extract the string tension a (i.e. the slope of the linear segment) of the 
quark-antiquark potential and compare it with the value known from lattice gauge calculations. 
In Greensite's book [23], the value obtained from LQCD is quoted to be a w 0.18 GeV 2 . In Figure 
3.3, the approximate values of the slopes are 0.16, 0.13, 0.11 GeV 2 respectively. These derived 
values are in reasonable agreement with the LQCD calculations. Thus, the present results are 
gratifying given the approximate knowledge of the parameters A and R. 

Note that in the approach of this dissertation, the effect of the non-Abelian terms on the inter- 
particle interactions is to modify the shape of the potential but leave the coupling constant a s 
unchanged. Within a perturbative calculation, the inclusion of virtual quark-antiquark pairs would 
give rise to vacuum-polarization effects. Ultimately, these would lead to a distance-dependent 
coupling constant a s that would modify the strength of the Coulomb potential. 

This is quite different from the perturbative S-matrix formalism where the coupling constant 
a s becomes energy dependent but the Coulombic shape of the potential remains unchanged. 
Nevertheless, it is gratifying to see that the inclusion of a virtual quark-antiquark pair in the trial 
state (3.86) changes the potential between the heavy valence quark-antiquark pair from purely 
attractive Coulombic one to one which exhibits linear confinement. 
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Figure 3.2: The quark-antiquark potential V x + V 2 for A = 1.0 and \R\ = 3.0 x 10~ 12 . The solid 
line shows the Coulombic V\ contribution only. The two sets of plotted points corresponds to 
the two choices of the overall sign in V 2 as indicated. The lower curve (labelled "plus") is clearly 
unphysical. 




Figure 3.3: The quark-antiquark potential V\ + V 2 for three choices of the parameters (A, \R\) 
with the appropriate (-) sign in V 2 . The choices are I: (0.8,9.9 x 10~ 13 ), II: (1.0,3.1 x 10~ 12 ) and 
III: (1.2,7.7 x 10~ 12 ). 
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4 Concluding Remarks 



... until non-abelian gauge theories are 
solved analytically, there is likely to be 
disagreement about ... and the origin of 
the confining force. 

Jeff Greensite (An Introduction to the 
Confinement Problem ) 

The research presented in this dissertation is concerned with the derivation of the inter-particle 
potentials in few particle systems in QFTs with non-linear mediating fields. The derivation 
of the potential is implemented in the Hamiltonian formalism of QFT. The method employs 
a reformulation of the original Hamiltonian to enable the usage of simpler trial states. The 
reformulation follows from the classical equations of motion and the Greens function for the 
mediating field. The equation of motion for the mediating field is non-linear and has no known 
analytical solutions. In the absence of other readily available alternatives, an iterative technique is 
used to approximate the solution to the leading order in the coupling constants of the underlying 
theory. 

The iterative technique is only valid in the perturbative regime. Therefore, the developed 
method of deriving inter-particle potentials is not universal. To extend its applicability to the 
non-perturbative regime would require an ingenious method of solving equations (2.4) and (3.51) 
possibly using a variational technique that would reformulate the Lagrangian densities (2.1) and 
(3.27) but differently. At this point, it is only a speculation and suggests a possible direction of 
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continuing this research. 

To apply the described method, in its current form, requires the calculation of matrix ele- 
ments of the Hamiltonian in the context of few particle systems. With the reformulation, it is 
not necessary to include mediating field quanta in the trial state, unless one wishes to include 
processes that involve their physical (as opposed to virtual) emission and absorption. Application 
of the variational method leads to relativistic equations where the interactions among the parti- 
cles are represented by kernels which are accurate to the leading order in the coupling constants. 
However, in order to capture as much relevant physics as possible one is compelled to employ 
multi-component trial states such as in equations (2.66) and (3.84). The multi-component trial 
states lead to coupled relativistic equations which are more difficult to analyze. The amount of 
algebra increases drastically. This is unavoidable and a substantial complication of the method. 

Once the relativistic equations are established, they are then reduced to the non-relativistic 
limit and Fourier-transformed to the coordinate representation. In coordinate representation, the 
interaction kernels become inter-particle potentials which, due to the complexity of the non-linear 
interactions, are expressed in the form of multi-dimensional integrals. These integrals are solvable 
using numerical methods. For the systems described by multi-component trial states, one requires 
the complete knowledge (i.e. the exact wave function), in general, of all but one component. Thus, 
one is forced to postulate ansatze which contain variational parameters to provide a reasonable 
representation the system. This complication is another drawback of this method since it requires 
a tedious calculation to determine the optimal values of the variational parameters. 

The result of the chapter devoted to the scalar model with a Higgs-like mediating held is that 
there is no confinement for few particle systems described by the trial states (2.33), (2.52) and 
(2.66) in the case when /x > 0. The non-linear interaction terms (2.14) and (2.15) do not raise the 
inter-particle potential to the region where energy is positive. On the contrary, the inclusion of 
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these terms seems to lower the inter-particle potential as can be seen in Figures (2.3) and (2.9). 
For the \x = case, the inter-particle potential must be regularized to be finite. This regularization 
introduces an infinite constant in the shifted energy of equation (2.47) and hence is considered 
unphysical. It is uncertain whether the appearance of this infinity is due to the inadequate choices 
of ansatz or the possibility that the theory given by the Lagrangian density (2.1) is ill-defined 
when \i = 0. 

The reformulation of the QCD Lagrangian is readily generalizable from the scalar Higgs-like 
model. Actually, the reformulation is more appealing in QCD than in the scalar Higgs-like model 
since it does not involve an intermediate Lagrangian as that of equation (2.8). On the other hand, 
the reformulation of QCD requires a choice of gauge. Ideally, the choice of gauge should not alter 
the final result. Notwithstanding, a computational technique usually loses its attractiveness when 
the gauge invariance becomes broken. The gluon propagator in the Feynman gauge, equation 
(3.52), leads to the reformulated Lagrangian terms given by equations (3.58)-(3.60). It would be 
quite interesting to see what these terms look like with no gauge fixing and how the results would 
be affected by different gauge choices. 

The vanishing of the colour index in the calculation of the matrix element of the Hamiltonian 
for the three quark system shows that the mono-component trial state (3.76) is insufficient to 
describe this bound state properly. To observe the effects of the non-linear terms on a three 
quark system would require constructing a three plus five particle trial state. The reason why 
a simple three quark trial state fails is probably due to the fact that it does not contain any 
virtual pairs and cannot accommodate string breaking. It is the inclusion of a virtual pair, as 
seen in the analysis of the two-component trial state (3.84), which invokes the extra interaction 
that induces confinement. Recall that there are no gluon bound states observed in nature despite 
the fact that there are LQCD predictions of them [2]. Therefore, this suggests that it is primarily 
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the virtual quarks, and not gluons, which provide the means for the strong interaction to trap 
valence quarks in strongly bound states. In fact, without explicitly using any gluon quanta in the 
trial state (3.84) (though virtual gluons are still represented by their propagators), the derived 
quark-antiquark potential shown in Figure 3.3 shows evidence of linear confinement. 

Perhaps, the most important results of this dissertation are reflected in Figures 3.2 and 3.3. 
These figures imply the possibility of linear confinement in quark-antiquark systems. The devi- 
ation from a linear form of the potential for large separations, in addition to inaccuracy of the 
numerical evaluation of the triple quadrature in y 2 equation (3.126), is from the limitations of 
the variational ansatz given by equation (3.84) (although, this is suspected to be a small effect). 
To improve upon this, one might consider more flexible forms for the functions F and G contain- 
ing more non-linear parameters. This should translate into more accurate and realistic potential 
curves for larger separation distances (i.e. perhaps a longer linear segment). The drawback of 
more flexible forms for the functions F and G is that it would require more variational parameters 
which have to be thoroughly calculated, and not estimated, if an improvement on what has been 
done here is sought. 

The derivation of the quark-antiquark potential for a heavy meson produces an estimate of the 
string tension a. The string tension is frequently written as a = k Aqcd, where k is a constant to 
be determined and Aqcd is the renormalization-scheme dependent QCD scale. In perturbation 
theory, the determination of Aqcd requires QCD theory and experimental input. In LQCD, the 
determination of the string tension a docs not require knowledge of Aqcd- In this dissertation, 
the string tension is similarly obtained directly from the linear form of the confining potential 
with the use of the coupling constant a s and the current quark masses as inputs. 

From the usual perturbative treatment of QFT, one expects loop effects to dominate the 
mechanism of confinement. In the present formulation, the Feynman-like diagrams for the cubic 
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and quartic interactions, shown in Figures 2.2 and 2.7, and consequently the matrix elements do 
not posses these properties. However, the relevant loop effects are believed to be incorporated by 
the reformulation of the Lagrangian density. One should also note that in the context of LQCD, 
the exclusion of the fermion propagator in the calculation (the quenched approximation) does not 
alter the potential drastically [48]. 

In this dissertation, the light quarks in the two plus four QCD calculation have been treated 
non-relativistically. This is acceptable since in the heavy quark-antiquark equation (3.107) the 
kinetic energies of the light quarks do not appear and the non-relativistic approximation is suf- 
ficient in the potential energy. Recall that in the non-relativistic limit, the potential terms have 
been expanded to the leading order in inverse powers of the mass tub- 

A main feature of the reformulation of the Lagrangian density is the non-inclusion of the 
gluon field quanta explicitly in the trial state (though they are included implicitly through the 
propagators). It is then prudent to wonder whether the inclusion of the gluon field explicitly is 
required at all to describe bound states of QCD. Certainly, the non-linear gluon interaction terms 
in QCD are necessary to explain confinement. However, is it possible to devise new non-linear 
Lagrangian densities without the gluon field, based on new symmetries and principles, which 
would copy their effects and describe the phenomenon of quark confinement? 

One possible direction of research to achieve this might have to do with the gauge indices 
of QCD. It is simple to notice that the terms of the QCD Lagrangian (3.27) are completely 
factorisable into products of the colour index quadratic forms and the quark and gluon fields. 
Such factorization between the momentum and spin dependencies does not occur in the QED and 
QCD Lagrangian densities. One might then consider exploring the idea of uniting the colour, spin 
and momentum dependencies by modifying the Lagrangian density in unconventional ways. Of 
course, the guiding principle, in such an approach, should be that in the appropriate limit the 
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standard QCD Lagrangian is retrieved. 

In this dissertation, the results of the QCD calculations show evidence of linear confinement in 
quark- antiquark pairs and this is consistent with LQCD. Besides the potential itself, it is desirable 
to calculate the bound state energy of quark-antiquark pairs using the trial state (3.84). This 
would require the calculation of all remaining interaction kernels in equations (3.92) and (3.93) 
and an extensive numerical toil. Even more so, it is prudent to perform a relativistic calculation 
such that, more common, light quark-antiquark states could be examined. The research was 
conducted using several approximations and simplification which could be improved upon in the 
future, as discussed in this section. The author hopes that the work recounted in this dissertation 
can be useful to other researchers. 
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5 Appendix A: Higgs-like Scalar Model 



5.1 Three-Particle State 

Some intermediate steps of the derivations with the three-particle trial state (2.33) are pre- 
sented in this section. All matrix elements are calculated using Maple [49]. 

Derivation of the kernels 

The matrix element in the Schrodinger picture for the three identical particle trial state (2.33) 
is comprised of 



(* 3 | H-E |* 3 ) = (*s| H+-E |tf 3 ) + (* 3 | H h |*3> + (*s| H h |* 3 ) 



(5.1) 



where the contributions are 



(tf 3 |#*-£|*3> = J dpi,2,3(w Pl +Wp 2 +Wp3--B) ^*(pi,2,3)^s(pi,2,3), 

d Pl,2,3 d Pl,2,3 



(y 3 \H h |* 3 )=- 



(* 3 |ff/ 2 |* 3 ) 



x 5(pi + p 2 - pi - p 2 ) <5(p 3 - p 3 ) 
dpi, 2.3 dPl,2,3 



^(Pi,2,3)^(pi,2,3) 
1 



x <5(pi + P 2 + P3 - Pi - P2 - P3) 

1 1 



M 2 - ip'i -Pi) 2 _ 

^(P1,2,3)^(P1,2,3) 



M 2 - (p'i + P 2 ~ Pi ~ p 2 ) 2 A* 2 - (pi - Pi) 2 M 2 - (p' 2 - P 2 ) 



(5.2) 



(5.3) 



(5.4) 
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In working out the variational derivative the following identity is used: 

S sf(qi2 3 ) = ? ^(Pii -qi) 5 (Pi 2 -q2)^(pi 3 -qs), (5.5) 

1:2,3 ' i1.i2.i3 

where the summation is on the six permutation of the indices i\, ii and i 3 . 

The relativistic kernels for the three-particle trial state are given in equations (2.36) and (2.37) 
in the body of the dissertation. In the non-relativistic limit the kernels become 

g 2 r£(pi+P2-Pi-P2)<Kp(j-p 3 ) 



^3,3 (Pi,2,3> Pi, 2, 3 ) — 



.(2ir) 3 m 2 



M 2 + (Pi - Pi) 2 



5(p'i + P 3 ~ Pi ~ P3) S(p' 2 - p 2 ) S(p' 2 + p' 3 - p 2 p 3 ) S(p[ - pi) 
M 2 + (p 3 - P3 ) 2 + M 2 + (P 2 -P2) 2 

C3 3 (Pi,2,3,Pl,2, 3 )-- 4( J )6 ^ 3 

5(Pi + P 2 + P 3 - Pi - P2 - P3) 



(5-6) 



(5.7) 



(M 2 + (Pi - Pi) 2 )(M 2 + (P' 2 - P 2 ) 2 )(^ 2 + (P 3 - Ps) 2 ) 

where the symmetry property of Fs is used to simplify the expression. 

3 

To obtain equation (2.39), the kernel C 3 . 3 is multiplied by J^J c lpi ' x ' and the variable of 

i 

integration is shifted by p[ — pi = q^. Integrating over the vectors leads to equation (2.40). 
This integral is found using the standard technique where the radial integral is evaluated in the 
complex plane using Cauchy's integration formula. 

Gaussian parametrization 

An alternative method for evaluating the cubic contribution Vc for the three identical particle 
system is presented in this section. The invariance of Vc under rotations and translations of the 
coordinates becomes explicit in this method. 

The denominators of equation (2.39) can be written in a way which enables one to perform 
Gaussian integration over the momentum variables. The following identity and the Gaussian 
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integration formula are used: 



f 

Jo 



d/3e 



for A > 0, 



n dxi 

exp 







" K~ 


2 


L jX j - W 


= Dct 


.2tt. 


exp 



(5.8) 
(5.9) 



where K is an invertible and symmetric matrix, L is a vector and W is a constant. 

The three-particle trial state (2.33) yields the cubic interaction kernel (2.39). After integrating 
over the delta function in (2.39) there remain only two momentum integrations. Using the identity 
of equation (5.8) on each factor separately, the the following expression is obtained: 



Vc(xij,fj, > 0) = -a v 7r 3 J d/3i,2,3 J dqi,2 



-(^ 2 +q?) /3l-(M 2 +q2)/32-(Ai 2 + (qi+q2) 2 ) /33„-iqi-x 2 i p -iq2-x 3 i 



(5.10) 



where Xjj = Xj — Xj are the inter- particle vectors. Upon expanding the squares in the exponentials 
and defining q = qi +qi, a six dimensional Gaussian integral is identified. The matrix K and its 
inverse K~ x , the vectors Lj and W for the Gaussian integral are given in block diagonal form 



K 



2(/Ji+/3 3 ) 2/3 3 
2/3 3 2(/3 2 + /? 3 ) 

L = % 



J123 



2(/3 2 +/? 3 ) -2/3 3 
-2/3 3 2(^+/3 3 ) 



X21, X 3 i 



(5.11) 
(5.12) 



where /?i 23 = /?i/?2 + /?i/3 3 + (faP3- Applying equation (5.9), one arrives at 

V(*+*+/fe) / Ax 2 i + / 3 2X 2 i + /3sX 2 2 



Vc(xy,M > 0) = -o^tf 3 / 

Jo 



2,3 



,3/2 
Pl2 3 



exp - 



. (5.13) 



4/3i 23 

This is the expression given by equation (2.42). It shows explicitly that the cubic potential term 
Vc depends only on the inter-particle distances. Numerical integration over the parameters /3j is 
required to complete the calculation. 
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Derivation of the cubic potential energy Vc for x 2 = x 3 

The cubic inter- particle potential for the three-particle trial state (2.33) in the case when \i > 
and x 2 = x 3 is determined from the following expression: 

/Q~^ v g- 2 Ml v + x 2l| 
dv : ; .J", (5.14) 
v |v + x 2 i| 2 

where x 2 i = x 2 — xi and the non-vector notation means the magnitude of the vector. After the 
trivial azimuthal integration one obtains: 

roo p+1 e -^i|v+x 2 i 

Vc{x2ufj) = -2ir 4 a n / dvv 2 / dw t. ; r , (5.15) 

Jo J-i v |v + x 2 i| 

where w = cos 9 with 6 being the polar angle of the vector v. The polar integration is performed 

using the substitution g 2 — x\ x + v 2 + 2x 2 i v w, with @\ = (x2i+v) and g 2 — \x2i — v\ as the upper 

and the lower limits of integration. Thereupon, the integral of equation (5.15) can be written as 

' Q2 

To integrate over the variable g, one splits the interval of integration accordingly and ends up 
with the following radial integral: 

Vc(x 2 i,n) = - 2lT 

{r X21 dv r°° dv 1 

This integral is expressible in terms of the exponential integral defined by equation (2.43). Thus, 
the cubic potential term for the three-particle trial state in the case when \i > and x 2 = x 3 
evaluates to 

V C (X21,») ^ -^^\e- x ^[ln(3) - E 1 (x2i^)}+c^^E 1 ('Sx2i^ 1, (5.18) 
which is equation (2.44). 
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V C {X2UH) = - 2 -^ 1 [°°-e- 2 ^ [ ei dge-™. (5.16) 

X21 Jo V J g2 



Derivation of the cubic potential Vc for the equidistant case with \i = 

It is instructive to consider the inter-particle potential for the three-particle system described 
by the trial state (2.33) in the case when fi — and the coordinates are at the vertices of an 
equilateral triangle. Even though the ji = limit leads to an unphysical result, the technique of 
integrating out the angular coordinates is interesting in itself. 

The cubic contribution Vc to the potential is determined from the following expression: 

Vc(xi, 2 ,3,M - 0) = -a v 7T 3 f dv- — (5.19) 

J v |v + x 2 i| |v + x 3 i| 

A cut off on the upper limit in the radial integral is inserted in anticipation of an infinity. This 
yields 

Vc(xi2 3,M = 0) = - lim a r ,7r 3 / vdv [dfl- - ■■ - r (5.20) 

JO J l x 21 + V| |X 3 1 + v| 

where dv = v 2 dv dQ, is the volume element in spherical coordinates. The angular integration is 
performed using a number of steps. First, the |r! — r 2 | _1 expansion identity is applied on the 
remaining denominators 

i 00 i 

K^r^^ p ' M (5 ' 21) 

where r<(r>) is the smaller(larger) of the lengths of ri and r 2 , and Pi are the Legendre polyno- 
mials. This leads to 

oo oo R h s i 2 

Vc(xi, 2 ,3,M = 0) = - lim a„7r 3 ^ V / dvv / P (l (a 2 ) P ; > 3 ), (5.22) 

where r < / > is the smaller/larger of x 2 i and v, s < f > is the smaller/larger of Xsi and v, a 2 is the 
angle between the vectors x 2 i and v, and 0:3 is the angle between the vectors x 3i and v. 

Next, the Legendre polynomials are rewritten in terms of the spherical harmonics and in- 
tegrated using the orthogonality relation of the spherical harmonics. The relevant expressions 
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arc; 



A 

Pl{a)= 2lTT £ Yi m *(0iAi)Yr(02,<t>2), (5.23) 

m=— I 

where (#i,0i) and (#2,^2) specify the directions of two unit vectors which are separated by the 
cosine of the angle a, and 

<W 6u> = j dn Yr*(6, 4>) Y/(6, cf>). (5.24) 

Using these two expressions in equation (5.22) enables one to perform the angular integration. 
Then, resumming over the index m leads to the Legendre polynomials with the argument being 
x 2 i • X31 where x^ = Xjj/|xjj|. The resulting expression is 

Vc(xi,2,3,M = 0) = - lim ar,TT 3 V — — Pz(x 2 i • x 3 i) / dvv (5.25) 

One can set x 2 i < x 3 i without loss of generality and break the intervals of integration into three 
segments. The subsequent expression holds for arbitrary coordinates and fi = 0: 

p ; (x 2 i -x 31 ) J , V 21 



V ( frt V A 4 V* P '(X21-X 31 ) I 



1 ^ I, •"J u '21 x 31 

21 J- 31 ^2:31 J 

The free variables in this expression are the coordinate separations x 2 i, £31 and the quantity 
x 2 i • X31. This last quantity can be written in terms of X21, £31 and x 32 , thus showing that the 
inter-particle potential depends on the three distances Xij only. Specializing to the case where the 
coordinates are at the vertices of an equilateral triangle leads to a simpler expression to evaluate: 

V c (x 2 i = x 31 = x 23 ,^ = 0) =V c (A,fi = 0) = 

where A = |x 2 — xi| = |x 3 — xi| = |x 3 — x 2 | and x 2 i • x 3 i = 0.5 for the present equidistant case. 

The integrals in equation (5.27) are elementary. The I = term in the summation has to be dealt 
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with separately so that 

Vc(A,/i = 0) = -4^V, 



1 , , 4 ^ ft (0-5) 1 \ 



(5.28) 



2 ' V A/ "" "'' 2 W f + 1, 

The summation is expressible in analytical form using identities involving the Legendre polyno- 
mials [50] . The expression can be simplified to 



Vc(A,^ = 0) = -4 7rV, 



'- bii £)- in (! 



(5.29) 



5.2 Four-Particle State 

Some intermediate steps of the derivations with the four-particle trial state (2.52) are presented 
in this section. All matrix elements were calculated using Maple. 



Derivation of the kernels 

The matrix element in the Schrodinger picture for the four identical particle trial state (2.52) 
is comprised of 

(* 4 | H — E |tf 4 ) = (* 4 | H^-E |tf 4 ) + (* 4 | H h |* 4 > + (* 4 | Hi 2 |*4> + (* 4 | H h |* 4 ) (5.30) 
where the contributions are 

(* 4 |i?0 -E\^ A ) = j dp!. A F*( Pl . A )Fs (pi. a) (w Pl + cj P2 + w P3 + w P4 - E) , 

? f dp[ 4 dpi. A 

2tt)3 J ' 



(*i\H h |* 4 ) =- 



:^(pi..4)^(pi.. 4 ) 



(* 4 |^J 2 |*4) 



16(2tt) 3 J yw p; w p ^ Pl Wp 2 
x 6(p[ + p' 2 pi p 2 ) 5(P3 - p 3 ) S(p' 4 - p 4 ) 
dpi 4 dpi. .4 



M 2 - (p'i -Pi) 2 



, (5.31) 



g 3 r] 



8(2tt) 6 7 ^w p; Wp^p^ Pl Wp 2 Up 3 
x 5(pi + p' 2 + P3 - Pi - p 2 - P3) <5(P 4 - P4 



^(p' 1 .. 4 )F s (p 1 .. 4 ) 



A* 2 - (pi + P 2 - -Pi - -P2) 2 M 2 - (Pi ~ Pi) 2 M 2 - (P2 - -P2) 



(5.32) 
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and 



(vP 4 |iJ 73 |v? 4 ) 



2^) 9 J 



dpi 4 rfpi..4 



16(2tt) 9 7 ^w pi CJ p / £^7 CJ p / w pi uj P2 uj P3 cj P4 

X <Kp'l + P2 + P3 + P4 - Pi - P2 - P3 - P4) 



^(p' 1 .. 4 )^ S (Pl..4) 



1 



A* 2 - (K + P*2 + P3 -P1-P2- P3) 2 
1 1 



- (Pi - Pi) M 2 - (P 2 - P2) 2 M 2 - (P3 - Ps) 2 
In working out the variational derivative the following identity is used: 

n a 

SF s (pi. A ) 



(5.33) 



^(qi.. 4 ) 



(5.34) 



where the summation is on the 24 permutation of the indices ii, i<i, 13 and i±. 

The relativistic Yukawa and cubic interaction kernels for the the four-particle trial state (2.52) 



are 



34,4 (pi.. 4 ,Pi.. 4 ) = - 



16(2tt) 3 



E 



<5(Pl +P 2 - P*i - P* 2 )<*(P3 - Pi 3 )^(P4 ~ Pu) 



V 2 - ip'l - Pil ) 2 



(5.35) 



r // \ a 3 v sr ^(pi + P2 + p'3 - Ph - Pi2 - Pt 3 )S(p' 4 - Pu) 

C 4,4 (Pl..4.Pl-4j - - 0/ - _^ 6 2 < 



8(2tt) 



1 1 



V - (Pi +P2 - Ph -Pt 2 ) 2 V 2 - (Pl-Pii) 2 M 2 - (P2 -J^) 2 . 

The relativistic quartic kernel for the four-particle trial state is given in equation (2.55). 



(5.36) 
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In the non-relativistic limit the interaction kernels for the four-particle trial state reduce to 
y (Pi..4.Pi..4) 



S (Pl + P2 ~ Pi ~ P2) <?(P3 - P3) S(p' 4 - p 4 ) 

" /' 2 • <:p' ; Pi) 2 ' ' 



4(27r) 3 m 2 

+ <*(pj + P3 ~ Pi - PjO S(p' 2 - p 2 ) S(p' 4 - p 4 ) + 5(p[ + p 4 - Pi - p 4 ) S(p' 2 - p 2 ) 6{p' 3 - p 3 ) 



M 2 + (p'i-Pi) 2 M 2 + (pi-Pi) 2 

<*(p 2 + P 3 ~ P2 ~ Pa) <5(p'i - Pi) Hp'i - Pi) , S (P2 + P 4 - P2 - P 4 ) 5(p[ - pi) 5(p' 3 - p 3 ) 
M 2 + (P 2 -P2) 2 

+ S(p' 3 + P 4 - P3 ~ P4) ^(Pi - Pi) S(p' 2 - p 2 ) 



,i 2 + (p 2 - P2 ) 2 



M 2 + (P 3 - P3) 2 

nl , x 3g 3 77 <5(pi +P 2 +P 3 - Pi - P 2 - P3)<5(P 4 - P 4 ) 

^(Pl..4.Pl..4) - 77^6^3 ~ _ " _ _ 



(5.37) 



4(27r) u ro 



(M 2 + (P'1 - Pi) 2 )(M 2 + (P' 2 - P2) W + (P 3 - Pa) 2 ) 



+ 



+ 



<HPi + P2 + Pi - Pi - P2 P4) <5(p 3 - Pa 



(M 2 


+ (pi 


Pi 


2 )(M 2 + (P 2 


-P2) 2 )(M 2 + (p 4 ^ P4 ) 2 ) 




«(pi 


+ Pa 


+ P 4 - Pi - 


P3~P4)<5(P 2 -P 2 ) 


(m 2 


+ (pi 


Pi 


2 )(m 2 + (p: 3 


-Pa) 2 )(M 2 + (p 4 -p 4 ) 2 ) 




<KP 2 


+ Pa 


+ P 4 P2 


P3 - P4)<5(Pi - Pi) 



<9(p'l..4)Pl..4) = - 



3 5 4 a 



2(27r) 9 m 4 

S(p[ + P2 + P3 + P4 Pi P2 ~ P3 - P4) 

(P 2 + (Pi " Pi)W + (P' 2 - P2) 2 )(M 2 + (P 3 ~ P 3 ) W + (Pi - P4) 2 ) 



(5.38) 



(5.39) 



Gaussian parametrization 

An alternative method for evaluating the cubic and the quartic potential contributions Vq and 
Vq for the four identical particle state is presented in this section. The invariance under rotations 
and translations of the coordinates becomes explicit in this method. 

The denominators of equations (2.57) and (2.58) can be written using a technique which 
enables one to perform Gaussian integration over the momentum variables. The relevant identity 
and the Gaussian integration formula are given in equations (5.8) and (5.9). 

The four-particle trial state (2.52) yields the cubic interaction kernel (2.57). The calculation 
of Vc for the four-particle trial state is identical to that of the thrcc-particle trial state. There are 
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basically four copies of the three-particle results with different inter-particle distances involved. 
The matrix K and its inverse K -1 are identical for all four terms and the same as for the three- 
particle trial state equation (5.11). The vectors Li pertaining to each term in equation (2.57) 
are 



L\ =i 


x 2 i, 


x 3 i 


, L 2 =i 


x 2 i, 


x 4 i 


L 3 =i 


x 3 i, 


X41 


, £4 = i 


X32, 


X42 



(5.40) 

where the subscript indicates the corresponding term in equation (2.57). The constant W is the 
same as in equation (5.12). Application of the Gaussian integration formula (5.9) leads to the 
result 



VcixijH > 0) 



Jo 



dpi, 



V(/3i+/32+/3 3 ) 



2.3 



, /3ixf 2 +/? 2 xf 3 +/?3xj 3 \ / /3 lX f 2 +/3 2 xf 4 + fex^ 4 \ 

/ Pj^+M±±^k \ + CXD A ftxi 3 +^x| 4 + /3 3 xl 4 \ 1 

From equation (5.41) it is evident that Vc depends on the inter-particle distances Xy = |xj — Xj| 
only. 

The four-particle trial state yields the quartic interaction kernel equation (2.58). The calcula- 
tion follows the same steps as for the cubic interaction kernel. Applying the identity of equation 
(5.8) leads to a 9 dimensional Gaussian integral. The matrix K and its inverse K~ x in block 
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diagonal form are 



K =2 



{Pi + Pi) Pi Pi 

Pi (fa + Pi) Pi 

Pi Pi (Ps + Pi) 



(5.42) 



K- 1 



1 



'20 



1234 



p23i -0304 -0204 
-0304 0134 -0104 > (••• 4: '»» 
—0204 —0104 0124 

where /3i 23 4 = 010203 + PiPzPi + PiP20i + 02030i and ijk = p i3 + p ik + jk . The vector L and 
the constant W of the Gaussian integration are 



L = i (x 4 i, x 42 , x 43 ) , W=(Pi+p 2 + Ps + Pi) At 2 . 



(5.44) 



Applying the Gaussian integration formula (5.9) and after some algebra one ends up with the 
expression 



V Q (x y ,^ > 



0)=a^ 2 / 
Jo 



(^1234) 3 / 2 



x exp — 



P 3 Pi^ 2 12 + P20iAz + P2p3^ii + PlPiAz + PlPlAi + P1P2A4 



4P 



1234 



(5.45) 



This expression shows explicitly that the quartic potential term Vq depends only on the inter- 
particle distances. Numerical integrations over the parameters Pi are required to complete the 
calculation of Vc and Vq . 



Derivation of the cubic potential energy Vq for xi = x 3 and x 2 = x 4 with \i > 

The quartic inter-particle potential term for the four-particle trial state (2.52) in the case 
when n > 0, xi — x 3 and x 2 = X4 is determined from equation (2.62), which upon imposing the 
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restrictions is the following expression: 

r 6 -2m|v| c -2^|v+x 2 i| 

^q(xi,2,i,2,M = 0) = Vq{x2uH = 0) = a a TT J d\r — — | V + X21 |2 • ( 5 - 46 ) 

Following similar steps as those leading to equation (5.16), one obtains the result 

V Q (x 1>2 ,i,2,n>0) = V Q (x 21 ,n>0) = — f ± e ~^v r^e" 2 ^, (5.47) 

X21 J v J e2 g 

where gi, g 2 and x 2 \ are as before. The result of the integration over the variable g can be 
expressed in terms of the exponential integral (2.43) and is given in equation (2.64). The remaining 
integrals have to be integrated numerically. 

5.3 Improved Particle- Antiparticle State 

Some intermediate steps of the derivations with the improved particle-antiparticle trial state 
(2.66) are presented in this section. All matrix elements are calculated using a Maple worksheet. 

Derivation of the kernels 

The matrix element in the Schrodinger picture for the improved particle-antiparticle trial state 
(2.66) is comprised of 

(* t | H-E |* t ) - (* 2 | H+ - E|tf 2 ) + (* 2 | H h |* 2 ) 

+ (* 2 \,H h 1*4) + (* 2 |i? /2 |* 4 ) + (^A\H h |* 2 ) + (tf 4 | Hi, |*2> 

+ <* 4 | % - E 1*4) + (* 4 | H h 1*4) + (* 4 | H l2 |tf 4 ) + (tf 4 | H l3 |* 4 ), (5.48) 

where the kinetic energy contributions are 

(* 2 | - E |*2) - \C F \ 2 J dpi, 2 ^*(Pi, 2 ) F(pi, 2 ) (w Pl + lo P2 - E) , (5.49) 
(* 4 | H^-E |*4) = \C G \ 2 J d Pl ..4 G*( Pl .. 4 ) G( Pl .. 4 ) (w pi + lu P2 + cu p3 + co pi - E) , (5.50) 
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and, the relevant interaction contributions are 

icviv 



8(2 
1 



J dp' li2 d Pl , 2 F*(p 12 )F( Pl , 2 ) 
1 



+ 



+ 



+P 2 - Pi -P2) 
V""pi u pi w pi w w 

1 1 

+ ■ 



M 2 -(p'i-Pi) 2 ' M 2 -(P 2 -P 2 ) 2 ' M 2 - (Pi + P2) 2 ' M 2 -(p'i+P 2 ) 2 



(5.51) 



(* 2 |i? 7l |*4 



C* F C G g 2 



4(2tt) 3 

5(p'i - Pi) <Kp 2 - P2 - P3 - P4) 



dp'1,2 d Pi..4 ^(p'1,2) G(pi.. 4 ) 

1 



7' 



a; p^P2 W P3 ti; P4 



1 



+ 



(tf 2 |ff/ 2 |*4> 



S (P2 - Pg) <5(pj Pi P3 ~ P4) 

y^p^pT^w^pT 



/ i2 - (P2 - V2f M 2 - (P3 + P4) 2 
1 1 



>, (5.52) 



C* F C G g 3 rj 
4(2 



G9 V f , / 



2 dpi.. 4 F*(pi i2 )G( Pl .. 4 ) 



A* 2 - (p'i - Pi) 2 M 2 - (P3 + P4) 2 _ 

8{p'l + P2 - Pi P2 P3 - P4) 



Wp'j Wp^W pi Wp 2 Wp 3 Wp 4 



+ 



+ 



^ 2 


-(p'l+p 2 


-Pi 


- P2) 2 M 2 


-(p'i+P 2 ) 2 M 2 


- (Pi +P2) 2 




1 






1 


1 


M 2 


-(P1+P2 


-pi 


-P2) 2 /i 2 


- (Pi - Pi) 2 M 2 


-(p 2 -p 2 ) 2 




1 






1 


1 


M 2 


- (pi - Pi 




- pi) 2 11 2 


- (Pi - Pi) 2 M 2 


- (P3 +P4) 2 




1 






1 


1 



[i 2 - (p 2 - p 2 - p 3 - Pi) 2 M 2 - (P 2 - P2) 2 M 2 - (P3 + P4) 2 



1 



(5.53) 



_/i 2 - (pi + p 2 +p 3 +p 4 ) 2 M 2 - (Pi +P 2 ) 2 M 2 - (P3 +P4) 2 . 

In working out the variational derivative expressions similar to those in equations (5.5) and 
(5.34) are used, except that here there is no permutation index to keep track of. The relevant 
interaction kernels are given below: 

v fW r, \ \ c f\ 2 9 2 f , ( , > <5(pj +P 2 - Pi - p 2 ) 

y2,2(Pl,2>Pl,2) = g , , 3 / dPi )2 F(p 12 ) = 

8(2^)^ J ' ' y/w^ovw^ 



+ 



+ 



pi^p 2 ^pi^p 2 
1 



/i 2 - ^ - pi) 2 ^i 2 - (p 2 - p 2 ) 2 ^i 2 - (pi + P2 ) 2 ^ 2 - ^ + p 2 ) 2 



(5.54) 
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y2, 4 (Pl.. 4 ,Pl,2) = - C ^Y j d Pl.AG(Pl. 



/ S(P[ - Pi) S(P2 + P3 + P4 - P2) 



V 1 



1 



1 



+ 



$(p' 2 - P2) S(p[ + P3 + P4 - Pi) 



a/ 



1 1 



+ 



/i 2 - (pi - pi) 2 p 2 - (p' 3 + p' 4 ) 



C2,4(P' 1 ..4,P1,2)=-^|^^ / dpi..4 G( P i.. 4 ) 



<*(Pl + P 2 + P3 + Pi Pi P2) 
1 1 



+ 



+ 



+ 



P 2 


- {p[ +P2-P1- 


-P2) 2 M 2 


-(p' 1+ ^) 2 M 2 


- (Pi +P2) 2 




1 




1 


1 


P 2 


- (K + p'i - pi 


-p 2 ) 2 M 2 


- (p'i-Pi) 2 V 2 


-(p' 2 -p 2 ) 2 




1 




1 


1 


P 2 


- (p[ +P3+P4 


-Pi) 2 M 2 


- (Pi - Pi) 2 M 2 


-(P3+P4) 2 




1 




1 


1 



A* 2 - ip'2 +P' 3 +P4- P2) 2 M 2 - (p'2 - P2) 2 M 2 - (P' 3 + P'i) 

1 1 1 

M 2 - (Pi + P' 2 + P'3 + Pi? V 2 - ¥1 +P2) 2 V 2 -(Ps+ Pa) 2 



(5.55) 



(5.56) 



The remaining interaction kernels are non- vanishing; they are just not included here. 

Derivation of the term V^i 

The derivation of the V^ 4 terms, as given by equation (2.82), begins with the substitution 
of equations (2.78) and (2.79) into the appropriate terms in equation (2.74). This substitution 
reduces to the following expression: 



/Gv 
d*3,i ^2.4(xi.. 4 )^7 



g(xi„ 4 ) 

Xl, 2 ) 



exp - 



|Xi 3 | + |X 23 



. (5.57) 



|Xl 3 | |X 23 | 

The dummy variable x 3 can be integrated out by shifting the integration variable x x — x 3 = z 
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and x 2 — X3 = z in the first and second terms respectively: 

jrY , s 1Ra 9 f , ( |Z| + |Z + X1 Z , ., 

V 2 4 {x 12 ) = ~ / dz —— cxp . (5.58) 



It is already evident that V 2 \ depends only on the inter-coordinate vectors. It is convenient to 
process in the spherical coordinates. The azimuthal integration is trivial and yields 2 n. The polar 
integral can be calculated with the substitution: 

z 2 + x 2 + 2zxa = /3 2 , (5.59) 

where z = |z| is the radial coordinate, x = |xi2| is the length of the separation vector and a is the 
cosine of the polar angle. Upon integration the new variable /?, one obtains the following radial 
integral: 



V^{x 13 ) = -^pL [dze-- **(-*) 
a 6 x J V a J 



x I - a {p 2 + a) cxp (-P^j+a (/?i + a) cxp (~^f) | (5-60) 

where (5\ = \z — x\ and j3 2 = z + x, are the familiar quantities from the previous calculations. To 

perform the radial integral in the parameter z, one has to split the interval of integration into 

two intervals in order to account for the absolute value in i.e. / = / + / The radial 

Jo Jo Jx 

integration is straightforward, although tedious, and can be calculated with Maple. The result 
comes out to be 



8 7T Rp~ Xl2 (Aii+l)/a 

V 2 \{x l2 )= „ r^2 (« ^12 V 2 e Xl2f * + 2 x 12 nc x »»-2c x »» + 2), (5.61) 

X12 a z /z z (fio + 2) z 



where, recall, the parameter a and R have to be determined through a variational calculation. 

Derivation of the term V 2 C 4 

The derivation of the V 2i terms, as given by equation (2.83), begins with the substitution 
of equations (2.78) and (2.79) into the appropriate term in equation (2.74). This substitution 
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reduces the term to the following expression: 

G(X!.. 4 ) 



V 2^4( x 12) = R J dx 3 .4 C 2 ,4(X1..4 I 

= J dx 3 dqi )2 



e -ix 3 -(qi+q 2 )gixrqi e !X2-q2 



exp 



X13 + |x 23 



. (5.62) 



(W(qi+q 2 ) 2 ) 

The momentum integrals in the expression above can be solved in the Gaussian parametrization. 
Leaving out the constants in front and using equation (5.8) in equation (5.62), one obtains 



(5.63) 



V 2 a ~ J dx 3 dq_i,2 d/3i,2,3 e" 

x e -(M 2 + (qi+q 2 ) 2 )ft c -(M 2 +q 2 )& c -(p 2 +q^)fe gxp ^ ]*13]+ jX 2 , 

where Xjj = Xj — xj are the inter-particle vectors. Upon expanding the squares in the exponential 
and defining q = qi + q 2 , a six dimensional Gaussian integral is identified with the following 
parameters in block diagonal notation: 



K 



2(01+ As) 2ft 
2ft 2(ft+ft) 

L= i 



K- 1 



2(ft+ft) -2ft 
-2ft 2(ft + ft) 



X 3 i, X 2 1 



4ft 23 

-zpi 

^ = (ft + ft + ft)^ 2 , 



(5.64) 
(5.65) 



where ft 2 3 = ftft + ftft + ftft- Applying equation (5.9) and making the change of variable 
x 3 = z + Xi leads to the following expression: 



P? 2 



J 123 



x exp 



ft X 2 2 + ft Z 2 + ft (Z + Xi 2 ) 2 |z| + |Z + Xi2 



->123 



(5.66) 



Already in its present form, it is evident that Vf 4 depends only on the inter-coordinate vector 
xi 2 and no angles. 

As before, length is expressed in units of the Bohr radius — - — , and, correspondingly, energy in 

mag 

units of m a 2 . Consequently, the contribution of V 2 C 4 can be written in terms of the dimensionless 
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variables r = x i2 ma g and w = zma g , and the dimensionless variational parameter A — ama g . 
Reinserting the constants in front and trivially integrating out the azimuthal coordinate of w, one 
obtains: 

^2,4^12 J - A 3 ft 3 / 2 

Pl23 



/ fh r 2 + /3 2 w 2 +/? 3 (w + r) 2 |w| + |w + r| ^ 
XCXP l ^ 1 J" (5 - 6?) 

There remains a five dimensional integral to perform which has to be numerically. 
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6 Appendix B: QED and QCD 



6.1 Particle- Antiparticle State in QED 



This section is reserved for the details involved in and related to calculations dealing with 
Dirac spinors. 

The following multiplication identities involving a single type of the Dirac spinor and the 
gamma matrices have been used: 



u(p,s)^iu(p',a) 
w(p',s)7o v(p,a) 



K 


+ m) (w p ' 


+ m) 


4 m 2 




+ to) (w p ' 


+ m) 


4 TO 2 




+ m) (w p / 


+ m) 


4 TO 2 


K 


+ m)(u;p 


+ m) 


4 TO 2 



<rp' 



cr • p 



(w p ' + to) (w p + m) 
CTi cr • p' cr • p Oi 



(ojp> + to) (w p + to) 
crp' crp 



(w P ' + m) (w p + m) 



<T • p' (Tj cr, , cr • p 



(w p / + m) (cjp + m) 



X<T, 



»7<7, 



(6.1) 
(6.2) 
(6.3) 
(6.4) 



The set of identities below involves both types of the Dirac spinors together with the gamma 
matrices: 



u{p,s)7°v(p',a) = 
u(p,s)Yv(p,a) = 



(w p + m)(wp> 


+ m) 


4 m 2 


(w p + m)(w p > 


+ m) 


4 m 2 



Xs 



Xs 



cr p' 



CT • p 



(w p + to) (w p + to) 



CT • p (T t CT • p 



(w p + to) (up + to) 
The complex conjugates of these equations are straightforward to obtain. 



(6.5) 
(6.6) 
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In the strict non-relativistic limit all of the above expressions reduce to 

u(p, s) 7° u(p', a)=S sa , (6.7) 

u(p,«)7 i «(p',a)=0, (6.8) 

v(p', s) 70 v(p, a) =S sa , (6.9) 

v(p',s) 7i v(p,a)=0, (6.10) 

u(p,s)7°«(p,<r)=0, (6.11) 

u(p, s) 7 4 u(p, a) = xl <?i V<r- (6-12) 

The fundamental and defining relation of the Pauli matrices is 

OiOj— ieijk&k- (6.13) 

The following, and frequently used, identity involves three Pauli matrices and can be obtained 
with a repeated use of equation (6.13): 

&j (Tj ak = Sij (Tfe — Sik a-j + 5jk <Ji + i €ijk- (6.14) 

6.2 Three Quark Trial State 

The contents of this section pertain to the derivations related to the three quark trial state 
(3.76). 

The matrix element in the Schrodinger picture of the trial state (3.76) with the reformulated 
QCD Hamiltonian (3.63) is comprised of 

(* 3 | H R -E |* 3 ) = (* 3 | - E |* 3 ) + (* 3 | H$ A 1*3) + (* 3 | H* A 1*3), (6.15) 

where, upon the integration over the momentum variables and summation over the spin and colour 
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indices, the contributions are 



(* 3 |ff*-£|*3> = J dp 1:2 , 3 F s * M (p 1 .. 3 )^, M (p 1 ..3)«+^ 2 +^ 3 -i?), (6.16) 
(*s| H$ A |* 3 ) = 2(|^3 E S / dpi - 6Fs *i^^ (pi - 3)jFls2 ^' t2(p4 - 6) 

colour spin 

x H^^^j tti 2 ,j 2 ,k 2 {L\ + L 2 + L 3 + L 4 + L 5 + Le), (6.17) 



(*3|^|*3) = ^-r fcc ^m B m c X: £ / 



dPl..6 



1 ' spin colour J V Pi P4 L "P2 L "p 5 u 'p3 U 'p6 / 

x ^7 1 ,ri,ti(Pl-3)-P , S2,r 2 ,t2(P4..6) ^iiji.fci ^i2,hM 

x <5(pi + p 2 + p 3 P4 P 5 - Pe) (Ti + T 2 + T 3 + T 4 + T 5 + T 6 ) , (6.18) 



where, the L-type terms are 



Li = S(pi +p 2 P4 Ps)<$(P3 - Pe) ' ' " 



(«« 5 )V2 

(p^ -pi) 2 



x S(pi, S i)7^(P4,52)S(p2,r 1 ) 7 ^(p 5 ,r 2 )^ il2 T« j2 ^ , (6.19) 



mAm B 



L 2 = 6( P1 + P2 P4 pb) 5( P3 - P6 ) ( A A B B)1/2 

V Pi ^p 4 P2 P5 / 

^fcifc 2 $tit 2 
(pi - pi) 



x u(p 2 ,r 1 )^u(p 5 ,r 2 )u(p 1 ,s 1 )ru(p4,s 2 )T^ ^ (6.20) 



L 3 = S(pi +p 3 - P4 - Pe)(5(P2 - Ps) ■ 



(««) 1/2 



X S( Pl , Sl )7^(p4, S 2)S(p3,tl)7^(P6,t2)T^ 2 r fc a ife2 ^| J ipC)2 ' ( 6 ' 21 ) 
£4 = S(pi + p 3 - P4 - Pe) S(p 2 - P5) ' ' 



(««) 1/2 



x S( P 3,t 1 ) 7 ^(p 6 ,t 2 )S(p 1 , Sl )7^(P4, S 2)T fe a ifc2 Tf iJ2 V^ r ' , (6.22) 

IPi — Pi ) 

L 5 = S(p 2 +p 3 - Ps - Pe)^(Pi -P4) ■ 



(iO B UJ B UJ C UJ C )V2 
rp2 P5 P3 Pt' 



8 S 

x S(p 2 ,r 1 ) 7AlU (p5,r 2 )S(p3,ti) 7 ^(P6,t2)r^ 2 T fe a ifc2 ( ^' 2 _^ 2 , (6.23) 
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Le = S(p 2 +P3 -P5 -Pe)<5(pi - P4) 



(iO B L0 B UJ C L0 C )V2 
V P2 P5 P3 P6/ 



x «(P3^i)7^(P6^2)S(p 2 ,r 1 )7^(P5,r 2 )T fe a ifc2 T^ 2 (6.24) 

IP 2 P5 J 



The colour factors are calculated, for instance, as follows: 



a '*i,ji,fei s 'i2,i2,fe2 - L i 1 i 2 - L jij2 " k i k 2 ~ p. tjl^'lfel t k 1 j 2 i 2 J -i 1 i 2 ± j 1 j 



a ($3132 ^»i»2 »2 ^'132) ^hio^jih 



= l(Tt 2i2 T^ 2 -T^ i2 Tt 2j2 ) = ~tr [T a T a ] = -| (6.25) 
where the fact that the generators are traceless and the following identity have been used: 

Qj/c €imn — &jm $hn $jn $km- (6.26) 

The T-type terms are 



Ti = u(pi, si)7„u(p 4 , s 2 )u(p 2 , ri)7 M u(p5, r 2 )u(p3, ti)7"u(p6, t 2 ) 

v (Pi -P4 )m /Tia ^Tl6 ^THC /(• nry\ 

A /-„/! „AN2^„B „B\2^„C „CN2 J «i«2 .71.72 1 k 1 k 2 i 



r 2 = u(pi, Si)7i/«(P4, s 2 )u(p 3 , t!)7 M w(p 6 , t 2 )u(p2, r 1 )7 ly u(p 5 , r 2 ) 

T? ii2 Tl k2 T* ij2 , (6.28) 



(^-^(p^-P^Of -pf) 2 4142 klk2 

T 3 = u(p 2 , ri)7^u(p 5 , r 2 )w(pi, si)7 M w(p 4 , s 2 )u(p 3 , ti)j"u(p 6 , t 2 ) 

K3A 2 n ik2 , (6-29) 



(Pi-Pi) 2 (pt-Pi) 2 (p§-P^) 2 ±nn Ml12 ±klk2 ' 



T 4 = u(p 2 , r 1 )7 ly u(p 5 , r 2 )w(p 3 , t 1 )7^u(p 6 , t 2 )u(pi, s 1 )7 1/ u(p 4 , s 2 ) 
(p?-P?)V 



^.,2^2^ (6-30) 



(pf-pf) 2 ^-^) 2 ^-^) 2 nn ^ 
?5 = w(p 3 , h)j v u(p 6 , t 2 )u(pi, si)7 Al u(p 4 , s 2 )w(p 2 , ri)7 l/ u(p 5 , r 2 ) 



(P3 P6 )m 71a rp6 ™ la 01 \ 



C _ „C\ 2 /'„A _ ^A^^B _ „B\ 2 ^^1^2 ii«2 .71.72' 



(p^-p^) 2 (pf-pt) 2 (p§-pf) 
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T 6 = u(p 3 , ii)7„u(p6, t 2 )u(p 2 , ri)7 M w(p 5 , r 2 )u(pi, Si)7"u(p4, s 2 ) 



G>3 -Pe% 



(P3 C -P^) 2 fe B -pf) 2 (pf-ff) 2 felfc2 3132 



(6.32) 



6.3 Products of SU(2) and SU(3) Representations 

The Young tableaux method provides a way to find the dimensions of products of group repre- 
sentations. The general rules of the Young tableaux method for particle physicists are described 
in the book by Georgi [51]. In this section, the Young tableaux for the product of SU(2) and 
SU(3) representations along with the corresponding dimensions of the resulting representations 
are shown. The products of the SU(2) representations refer to the addition of spin while those of 
SU(3) refer to colour. 

The Young tableaux and the dimensions of the representations corresponding to the product 
of the fundamental and conjugate SU (2) representations are 



□ ® □ 



2 ® 2 = 3 © 1. 



(6.33) 



Adding a fundamental SU(2) representation to the above yields: 



□ 



2®2®2 = 2©2©4. 



(6.34) 



Lastly, adding a conjugate representation to the the above yields: 



® □ 
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2<g)2<g>2<g)2 = lffilffi3ffi3ffi3ffi5. 



(6.35) 



The Young tableaux and the dimensions of the res-presentations corresponding to the product 
of the fundamental and conjugate SU (3) representations are 



a 






a | 


]L 




b 




3 


= i 


i 6 


3 1 



(6.36) 



Adding a fundamental SU(3) representation to the above yields 



3 = 3 



15. 



(6.37) 



Lastly, adding a conjugate SU(3) representation to the above yields: 









a 




b 











a 




b_ 









a 






b 





\a\ 


b 









a 




b 





3®3®3®3 = lel 



10 © 10 © 27. 



(6.38) 



The labels a and b are included because the representations associated with multiple boxes, in 
this case 3, must be incorporated with other representations with the assistance of these labels. 
This is described in detail in the book by Georgi [51]. 

The Young tableaux and the dimensions of the representations corresponding to the product 
of two fundamental SU (3) representations are 



□ ® □ 



3 <g> 3 = 3 © 6. 
Ill 



(6.39) 



(6.40) 



6.4 Multi-Component Quark-Antiquark Trial State 

The matrix element in the Schrodinger picture of the trial state (3.84) with the reformulated 
Hamiltonian operator (3.63) is comprised of the following contributions: 

(*t| H R -E |tf t ) - (* t | H^-E |* t ) + (* 2 | H$ A |tf 2 ) + (* 2 | H$ A |*4> + (* 2 | H* A |tf 4 ) 
+ (tf 4 | H* A |M> 4 ) + (M> 4 | H $ A |*2) + (* 4 | H£ A |*2) + (tf 4 | H* A |*4) + (* 4 | |* 4 ). (6.41) 

The free and the linear contributions to the matrix element are 

(* 2 |Jfy-i5|* 2 > = \C F \ 2 J2 / dPi,2CAi(Pi.2)^i,Ai(Pi,2)«+<--B) 
\C F \ 2 g 2 m A 



KlXl 

'- ,,2 ™2 

(* 2 |^|*2> 



2(2tt) e 



X u(pi, Ki)7m u (P3j k 2) U(P4, A 2 )7 A MP2, Al)- 



(P4 -P2) 2 

Tr(T a )Tr(T a ) 



- u(pi, Kl)7^(P2, Ai) U(P4, A 2 )7 A 'm(P3, k 2 ) 

- u(pi, «i)7^(P2, Ai) u(p4, A 2 )7 Ai u(p 3 , k 2 ) 

+ u(P4, A 2 )7M ,y (P2 ; Ai) u(pi, Ki)7 Ai u(p 3 , ^2)7^—^ ) , (6.42) 



(Pi + P2) 2 
Tr(T a )Tr(T a ) 

(P3+P4) 2 
Tr(T a T a ) 

(P3 - Pi) 2 

where the two virtual annihilation terms are set to zero by the colour factors (i.e. the factor 

coming from the summation of the colour indices) since the T matrices are traceless. The gluon 

4 

exchange terms have the colour factor of Tr (T a T a ) = —. 

O 
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The contribution to the matrix element from the linear cross term is 



x 

Kl,2 ^1,2 



(Pl-P3) 2 (P5+P6Y 



+ 



2(2tt) s 

^ y ^Pl..6 ^.Ai (Pi, 2)G K2 ,A 2 , W2 , I a 2 (P3..6) 

^1,2 M2,^2 

£(Pl ~ P3 ~ P5 ~ Pe) ^(P2 - P4) ft\iA 2 

fw^ w- 4 w s cv B ) 1/2 

\ Pi P3 PB P6/ 

x u(pi, Ki)7 I/ u(p3, k 2 ) v(p 6 , iy 2 )j' J u(p 5 ,^2) 

S(P2 - P4 - Ps - Pe) S(Pi - P3) S K1K2 
(uj a lo a lo b uj b ) 1/2 

V P2 P4 ps pel 

XW(P4 ) A 2 )7^(P2,A 1 )W(P 6 ,^) 7 ^(P5,M2) ( foj^p - ^5+^)2 ) }' ( 6 ' 43 ) 

where the colour factor is calculated as before. 

The contribution to the matrix element from the cubic cross term is 

,. T; 1 „ i.t. . i C* F C G g 4 m 2 A m B 
(y 2 \H 3A |* 4 ) - 

/" dPi..6 ^(Pi +P2 -P3-P4- P5 -Pe) / \n I \ 

X 2^ / { A A A A ~R RNl/2 i, Kl,A 1 (Pl,2)tx K2 ,A 2lAl2 ^ 2 (P3.. 6 ) 

«i,2 A ll2 J («««) 

M2^2 

x I (ps - Pi V «(pi, M)7^(P3, k 2 ) w(P4, A 2 )7^(P2, Ai) w(pe, ^)7"«(P5, M2) / afcc Tr(T a T c T b ) 

+ (Ps - Pi)„ «(Pi, «i)7^«(P3, k 2 ) «(P6, ^h"«(P5, /fc) «(p 4 , A 2 ) 7 ^(p2, Ai) / a6c Tr(T a T b T c ) 
+ (P4 - K)^ w(p 4 , A 2 )7^(p 2 , Ai) w(p6, ^h^Ps, /x 2 ) «(pi, Ki)7"«(P3, K2) / abc Tr(T a T c T 6 ) 
+ (P4 - P2) P w(p 4 , A 2 )7^(p 2 , Ai) u(pi, Ki)7"u(p 3 , /S2) w(Pe, ^)7"«(P5, Afc) / abc Tr(T a T b T c ) 
+ (PB + Pe)^ w(pe, ^2)7^(p 5 , M2) u(pi, Ki)y , «(p 3) k 2 ) w(p4, A 2 )7^(p 2 , A 2 ) / abc Tr(T a T c T 6 ) 

+ (PB + Pe) M w(P6, ^)7^«(P5, M2) w(P4, A 2 ) 7 M «(p 2 , Ai) u(pi, Ki)7^(p 3 , k 2 ) / o6c Tr(T a T b T c )| 

1 1 1 



(Pi - P3) 2 (P2 - P4) 2 (Ps + Pe) 2 



(6.44) 
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The colour factor, for instance, can be evaluated as follows: 



N c = f abc Tr(T a T b T c ) 






(6.45) 



where the anti-symmetric property of f abc and the commutator definition of the generating ma- 
trices (3.30) has been used. 

Similarly, one can obtain the other index combination 



Notice that a factor of i in the matrix element (6.44) and the i in Nc multiply out to give a real 

result as is expected in an energy calculation. The product f abc f abc is a positive definite quantity 

C F 

whose value can be absorbed into the variational parameter R = ——. Therefore, it is actually 
not required to calculate R. 

The contributions (vt^l H^ A and H^ A \^f 2 ) to the matrix element are just the complex 
conjugates of equations (6.43) and (6.44) respectively. Without actually making any effort one 
can see immediately that (vf^l H^ A Y&t) vanishes because of the colour factor as in equations 
(3.82). An extra delta function arising from the "spectator" quark would be of no consequence 
for the colour factor calculation in (vt^l H^ A The remaining contributions (vf^l H^ A 1^4) and 
(^4 1 H^ A 1^4) are irrelevant and therefore are not included here. 

The calculation leading up to equation (3.107) involves the non-relativistic versions of kernels 
(3.94) and (3.96) and the spin index function of equation (3.103). In particular, there are 
extensive summations over the spin and spinor indices. In order to perform them it is prudent 
to use Maple rather to do these summations by hand. The results are obtained in Maple and 




(6.46) 
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presented below where the ~ symbol refers to the fact that only factor carrying the spin and 
spinor indices have been included. 

The calculation involving the linear kernel 3^2,2 is 

@Kl,Al (3^,2)^1 (Pl..4) @K 2 ,A 2 

~ 0«!,Ai |u(pi,ki) j,,u(p 3 ,K 2 ) v(p 4 ,A 2 ) 7 /i w(p 2 ,Ai) |e K2 ,A 2 

= 0«i,Ai |u t (pi,Kl)u(p 3 ,'«2)w t (p4, A 2 )w(p 2 ,Ai) 

- u t (pi,Ki)7° 7 i «(P3,«2) w t (p4, A 2 )7° 7 J w(p2, Ai)| 6 K2 ,a 2 

- ~r~T- XL (<Ti <r • p 3 + <r • pi CTf) x K2 ??| 2 (er • p 4 ct, + er t <t • p 2 ) t] Xl \ <d K2M 
^m A ) 

, £«i Ai £k 2 A 2 + I , \ t I \ 

= 1 g — 2 Xki (ct, <T • p 3 + <T • pi CTi) Xk 2 ?? A2 (er • p 4 cr 4 + cr, er • p 2 ) T] Xl 

= 1 - 7-V f 3 (Pi • P2 + P3 • P4) - (Pi ' P4 + P2 ' P3) ) • (6-47) 
4m A V / 

The contribution involving summations in the kernel 3^,4 has two terms each distinguished by 
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a Roman numeral subscript: 

~ ©ki.Ai |s(pi,Kl)7 I ,u(p3,K2)S(p 6 ,^ 2 )7 ,yu (P5,A i 2)'5A 1 ,A 2 j =k 2 ,A 2iAI2 , i > 2 

= 6^,^! |u f (pi, ki) u(P3, K2) v f (P6, v 2 ) w(p5, a*2) 

- M t (Pi,'«i)7 7 I '«(p3,'«2)w t (p6,^2)7 7 4M (P5,A*2)| <5 Al A2 S K2A2iAl2 ^ 2 

©I ^Ai A 2 <^«i re 2 t / \ 
«i,Ai < vl 2 (o- • Pe + er • p 5 ) x M2 

<5Ai A 2 + , s 

- Xkr T • p 3 + <X ■ pi CTi) X K2 

t / q- ■ p 6 <jj cr ■ p 5 



X r lv 2 y^i + 4 m 2 J A 2 ,^ 2 ,iy 2 

<^A A 

= ~ 9ki ' Ai a X«! cr ' P3 + cr • pi ct 2 ) Xk2 

X ^ 2 (cr • p 6 CTi cr • p 5 ) 

X/i 2 ,=, K2.A 2 ,/J 2 ,^ 2 

= - 96m | mA (P3 ' (P6 x P5) + Pi ' (P6 x Ps)) + 0(m- A ) 7 (6.48) 

(@ki,Ai {y^AT^X^ 2 ^ 2 (Pl..6)S K2l A 2)tt2 ,i* 2 ) 

~ 6«i,Ai w(p 4 ,A 2 )7 l ,u(p2,Ai)w(p 6 ,^ 2 )7' / w(P5,M2)'5 Kl ^ 2 S K2 ,a 2 , M2!I , 2 

~ - ©ki.Ai o A 2 lA2 (cr • p 4 CT 2 + CT 4 cr • p 2 ) 7] Xl 77^ (cr • p 6 CT, cr • p 5 ) X M2 ^K 2 ,\ 2 ,p, 2 ,u 2 

o 771 ^ 777 4 

(p 4 • (pe x p 5 ) + p 2 • (pa x p 5 )) + C(m- 4 ). (6.49) 



96 

where similar steps have been used to obtained equation (6.49) as for (6.48). 

The contribution involving summations in the kernel C24 has six terms each distinguished by 
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a Roman numeral index: 



O^m)^^ 2 '" 2 (pi..e) S K2) A 2 ,^ 2lI / 2 ) 

~ 9«i,Ai (P3 -pi) M w(Pl,Kl)7i/W(P3,K2) 



x u(p 4 , A 2 ) 7^ v(p 2 , Ai) u(p 6 , i> 2 ) 7" w(p 5 , A«2) - K2 ,A 2 



^2^2 



~ ©ki.Ai ^(P3 -Pl)o^Ai,A 2 - ^~( P3 "POi^As (<T • p 4 <7j + (7; <7 • P2)??Aij 
< (<x • p 6 + er • p 5 ) Xa12 
1 t / 

- 2^ Xki °" ' P3 + ^ ' Pi Oi) X*2 

t ( _i_ _ ' Pe _ ' Ps ^ I w 

^2 ^ <T « ' 4m 2 / J = ' k 2,a 2 ,/j 2 ,i/ 2 

= 2^i^ (P3_pi) '( (P5+P6)x(P4_P2) ) 

+ 7^2-(P3 - Pi) ' ((Ps - Pi) x (p 4 - P2) - 2pi x (p 4 - p 2 )) + ©(m" 3 ). (6.50) 
The remaining terms can be determined using analogous steps: 

(©Ki.Ai ( C 2,4)^ A ,2 ' AI2 '" 2 (Pl..e) S K2) A 2)At2 , 1 ^) 

3 i i / 

= 75 2" (P3 - Pi) ' (P5 x Pe) + 7; 2" (P3 - Pl) • ( 3 (p 3 X p 4 ) + (p 2 X p 3 ) 

-(Pixp 2 )-(pixp 4 ))+0(m- 3 ), (6.51) 



(@«i,Ai O^M)^ 2 ^ 2 ^ 2 '" 2 (P1..6) S K2,A 2 , ( ii 2 , 1 / 2 ) iii 

3 i i / 

= 7, T (P4 " P2) ' (P5 X p 6 ) + 7) 2" (P4 - P2) • ( 3 (p 3 X p 4 ) + (p 2 X p 3 ) 

Zi TTlg Z TJX v 

- (Pi x pa) - (pi x p 4 )) + 0(m" 3 ), (6.52) 
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= 7J " (P4 - P2> ( (P6 + P5) X (Pi - P3) J 

+ tt-T (P4 - P2) • (-(P4-P2) x (p 3 -p!) + 2p 2 x (p 3 -pi)J +0(m- 3 ), 
2 m A V / 

(6.53) 



(^m)^ 2 ''* 2 ' 1 ' 2 (pi..e) S K2 ,a 2 ,^ 2 , 
3 i 

= 7i (P2 + P4) • (P5 X p 6 ) 

(Pi - Pa) • (2 p 2 x (p 5 + pa) + (p 2 - p 4 ) x (p 5 + pa)) + ©(m^ 3 ), (6.54) 



2m\ 



(®Kl,\l ( C 2a)Ti,\ 2,I12 ' V2 (Pl-S) 5 K2,A2,;*2,^) 

3i 

(Pi +P3) • (Ps x p 6 ) 



2 ITlA TUB 

+ 7T~~T (P2 - P4) • (- 2 pi x (p 5 + p 6 ) - ( Pl - p 3 ) x (p 5 + p 6 )) + 0(m" 3 ). 

(6.55) 

6.5 The Determination of the Optimal Variational Parameters 

This section shows the details of the calculation in equation (3.120). Recall that the function 
/ as given by ansatz (3.110), which is a properly normalized function. The calculation on the 
right hand side of (3.120) is straightforward: 

dp£.r(p).f(p) =£, (6.56) 

and, the kinetic energy contribution, calculated in spherical coordinates, is 

/ dpf p / p = / dpp 4 

J m A m a tt J (p 2 a 2 + V 4 



32 

m A 



Y 

■ 2 r a P 4 32 f n\_ m A a i 8 (Rvr . 

Va^Jo dP (pi + 1)4 - -m~A~^ V32 ) ~ ~A^> (b - b{) 

118 



where the substitution for a dimensionless variable of integration has been made and the con- 



sequent radial integral can be calculated using software like Maple or Mathematica, and the 

9 2 

constants arc A = a vha ol s and a s = - — . 

4-7T 

The contribution to the energy from the linear interaction term H^ A is divided into three 
terms 7i, T2 and %. The calculation of the first term is as follows: 



7;= / dPi,3.r(Pi){j 7 2, 2 (Pi, 3 )} i /(p 3 ) 

32 aV f 1 1 1 

57T 2 (2tt) 3 J dPl ' 3 (p2 a 2 + l)2 (p2 fl 2 + 1) 2 (p 3 _ pi )2 



3 

256 a 3 a 
T(2 



■ OL s f , , , 1 

7T 3 J (pf a 2 



1 



+ 1) 2 (p 2 a 2 + l) 



+1 



1 



{pI+P 2 i - 2pip 2 M) 



. (6.58) 



The evaluation of the \i integral can be done with the substitution v = p\ + p\ — 2pip 3 fi, where 
the limits of integration become: v\ = (ps + pi) 2 corresponds to /1 = —1 and v-i = (ps — p{) 2 
corresponds to fi — +1. The expression, upon the integration over the variable \i and the angular 
coordinates of the vector pi, becomes: 



Ti 



1024 a 3 a 



3(2 

1024 a s 
3a(2 7r) 2 

128 a s 
3 air 



1 a s f 



1 



1 



' PlP3 (p 2 a 2 + l) 2 (p 2 a 2 + l) 2 



In 



P3+P1 



P3 - Pi 



dpi dp. 



PiPi 



In 



P3+P1 



P3 - Pi 



dpi 



3 (p\ + l) 2 (p 2 + l) 2 

p\ 4a s 4 rriACt 2 

3p? + 3pf +P6 + 1 ~ "3T ~ 3 A ' 



(6.59) 



where the substitution for dimensionless variables of integration have been made and the radial 
integrals over the variables pi and p 3 have been found in Maple. 
The second terms T2 is calculated as follows: 



7~2 = J dpi, 3 /*(Pi) |3 ; 2,2(pi,3)} ii /(P3) 

as r 

a 3 it 2 J 



4 a 3 a 
3m 2 A 



P 2 3 



{pW + i? 



p\ 



(pi + 1) 2 J ^ 6 (pI + 1) 2 
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3 m 2 A a 3 



4 m A aj 
3 A 3 ' 



(6.60) 



where, once again, one can easily evaluate the radial integrals over the variables p\ and ps in 
Maple. 

Finally, the calculation of the third term % is as follows: 



% = /dpi,3/*(Pi) {y2,2(pi,3)} m /(P3) 



la 3 a s 



3m 2 A TT A 



p\ 



Pi 



aPl ' 3 (p 2 a 2 + l) 2 (p§a 2 + l) 2 (p 3 -pi) 2 



64 a 3 a 



3 771^ 



7T 2 J 



dpi,3PlP3 

4 16 Ci 4 



(pI+pI) 



3 a 3 7T 2 3 



niA a 



(p? + l) 2 (p| + l) 2 
4 16 d 



In 



P3 + Pl 



P3 - Pi 



3.61) 



where similar steps as those leading up to equation (6.59) has been taken to obtain this result, 
and the evaluation of the double-quadrature in Maple yields the constant 

pq(p 2 + q 2 ) 



Ci 



dpdq 



In 



p + q 



1.85055. 



(6.62) 



(p 2 + l) 2 (q 2 + l) 2 p-q 

The contribution to the energy from the cubic interaction term is discussed in the course 
of the dissertation. 
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